pardisoinit

Initialize Intel MKL PARDISO with default parameters in accordance with the matrix type.

Syntax

call pardisoinit (pt, mtype, iparm)

Include Files

  • mkl.fi, mkl_pardiso.f90

Description

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.

Note

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.

Note

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.

Input Parameters

perm

INTEGER

Ignored.

Note

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.

mtype

INTEGER

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

Output Parameters

pt

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.

Note

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.

iparm

INTEGER

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.

For more complete information about compiler optimizations, see our Optimization Notice.