Initialize Intel MKL PARDISO with default parameters in accordance with the matrix type.
call pardisoinit (pt, mtype, iparm)
- mkl.fi, mkl_pardiso.f90
This function initializes Intel MKL PARDISO internal address pointer pt with zero values (as needed for the very first call of pardiso) and sets default iparm values in accordance with the matrix type. Intel MKL supplies the pardisoinit routine to be compatible with PARDISO 3.2 or lower.
An alternative way to set default iparm values is to call pardiso in the analysis phase with iparm(1)=0. In this case you must initialize the internal address pointer pt with zero values manually.
The pardisoinit routine initializes only the in-core version of Intel MKL PARDISO. Switching on the out-of core version of Intel MKL PARDISO as well as changing default iparm values can be done after the call to pardisoinit but before the first call to pardiso.
Parameters types in this section are specified in FORTRAN 77 notation. See Intel MKL PARDISO Parameters in Tabular Form section for detailed description of types of Intel MKL PARDISO parameters in Fortran 90 notation.
This scalar value defines the matrix type. Based on this value pardisoinit sets default values for the iparm array. Refer to the section Intel MKL PARDISO Parameters in Tabular Form for more details about the default values of Intel MKL PARDISO
INTEGER for 32-bit or 64-bit architectures
INTEGER*8 for 64-bit architectures
Array with a dimension of 64. Solver internal data address pointer. These addresses are passed to the solver, and all related internal memory management is organized through this array. The pardisoinit routine nullifies the array pt.
It is very important that the pointer pt is initialized with zero before the first call of Intel MKL PARDISO. After that first call you should never modify the pointer, because it could cause a serious memory leak or a crash.
Array with a dimension of 64. This array is used to pass various parameters to Intel MKL PARDISO and to return some useful information after execution of the solver. The pardisoinit routine fills-in the iparm array with the default values. Refer to the section Intel MKL PARDISO Parameters in Tabular Form for more details about the default values of Intel MKL PARDISO.