pardisoinit

Syntax

Fortran:

call pardisoinit (pt, mtype, iparm)

C:

pardisoinit (pt, &mtype, iparm);

Include Files

The FORTRAN 77 interface is specified in the mkl_pardiso.f77 include file, the Fortran 90 interface is specified in the mkl_pardiso.f90 include file, and the C interface is specified in the mkl_pardiso.h include file.

Description

This function initializes 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 distributed by the University of Basel.

Note

An alternative way to set default PARDISO iparm values is to call pardiso 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 PARDISO. Switching on the out-of core version of 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

Note

Parameters types in this section are specified in FORTRAN 77 notation. See PARDISO parameters in tabular form section for detailed description of types of PARDISO parameters in C/Fortran 90 notations.

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 PARDISO parameters in tabular form for more details about the default values of PARDISO.

Output Parameters

pt

INTEGER for 32-bit architectures

INTEGER*8 for 64-bit architectures

Array, DIMENSION (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 PARDISO. After that first call you should never modify the pointer, as a serious memory leak can occur.

iparm

INTEGER

Array, dimension (64). This array is used to pass various parameters to 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 PARDISO parameters in tabular form for more details about the default values of PARDISO.