Developer Reference

  • 2020.2
  • 07/15/2020
  • Public Content
Contents

pardiso

Calculates the solution of a set of sparse linear equations with single or multiple right-hand sides.

Syntax

call pardiso
(
pt
,
maxfct
,
mnum
,
mtype
,
phase
,
n
,
a
,
ia
,
ja
,
perm
,
nrhs
,
iparm
,
msglvl
,
b
,
x
,
error
)
Include Files
  • mkl.fi
    ,
    mkl_pardiso.f90
Description
The routine
pardiso
calculates the solution of a set of sparse linear equations
A
*
X
=
B
with single or multiple right-hand sides, using a parallel
LU
,
LDL
, or
LL
T
factorization, where
A
is an
n
-by-
n
matrix, and
X
and
B
are
n
-by-
nrhs
vectors or matrices.
This routine supports the Progress Routine feature. See Progress Function for details. The case of
iparm
(24)
=10 does not support this feature.
Optimization Notice
Intel's compilers may or may not optimize to the same degree for non-Intel microprocessors for optimizations that are not unique to Intel microprocessors. These optimizations include SSE2, SSE3, and SSSE3 instruction sets and other optimizations. Intel does not guarantee the availability, functionality, or effectiveness of any optimization on microprocessors not manufactured by Intel. Microprocessor-dependent optimizations in this product are intended for use with Intel microprocessors. Certain optimizations not specific to Intel microarchitecture are reserved for Intel microprocessors. Please refer to the applicable product User and Reference Guides for more information regarding the specific instruction sets covered by this notice.
Notice revision #20110804
This notice covers the following instruction sets: SSE2, SSE4.2, AVX2, AVX-512.
Input Parameters
The types given for parameters in this section are specified in FORTRAN 77 notation. See Intel MKL PARDISO Parameters in Tabular Form for detailed description of types of
Intel® MKL
PARDISO parameters in Fortran 90 notation.
pt
INTEGER
for 32-bit or 64-bit architectures
INTEGER*8
for 64-bit architectures
Array with size of 64.
Handle to internal data structure. The entries must be set to zero prior to the first call to
pardiso
. Unique for factorization.
After the first call to
pardiso
do not directly modify
pt
, as that could cause a serious memory leak.
Use the pardiso_handle_store or
pardiso_handle_store_64
routine to store the content of
pt
to a file. Restore the contents of
pt
from the file using pardiso_handle_restore or
pardiso_handle_restore_64
. Use
pardiso_handle_store
and
pardiso_handle_restore
with
pardiso
, and
pardiso_handle_store_64
and
pardiso_handle_restore_64
with
pardiso_64
.
maxfct
INTEGER
Maximum number of factors with identical sparsity structure that must be kept in memory at the same time. In most applications this value is equal to 1. It is possible to store several different factorizations with the same nonzero structure at the same time in the internal data structure management of the solver.
pardiso
can process several matrices with an identical matrix sparsity pattern and it can store the factors of these matrices at the same time. Matrices with a different sparsity structure can be kept in memory with different memory address pointers
pt
.
mnum
INTEGER
Indicates the actual matrix for the solution phase. With this scalar you can define which matrix to factorize. The value must be:
1
mnum
maxfct
.
In most applications this value is 1.
mtype
INTEGER
Defines the matrix type, which influences the pivoting method. The
Intel® MKL
PARDISO solver supports the following matrices:
1
real and structurally symmetric
2