I was looking into how to call pardiso and I'm a bit puzzled by the argument pt.
The documentation says:
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.
The example provided with OneAPI 2024.2 does
INCLUDE 'mkl_pardiso.f90'
PROGRAM pardiso_sym_f90
USE mkl_pardiso
TYPE(MKL_PARDISO_HANDLE), ALLOCATABLE :: pt(:)
ALLOCATE (pt(64))
DO i = 1, 64
pt(i)%DUMMY = 0
END DO
This is coherent with the content of mkl_pardiso.f90, where MKL_PARDISO_HANDLE is defined as
TYPE MKL_PARDISO_HANDLE; INTEGER(KIND=8) DUMMY; END TYPE
The interface to pardiso is
INTERFACE PARDISO
SUBROUTINE PARDISO_D( PT, MAXFCT, MNUM, MTYPE, PHASE, N, A, IA, JA, PERM, NRHS, IPARM, MSGLVL, B, X, ERROR )
USE MKL_PARDISO_PRIVATE
TYPE(MKL_PARDISO_HANDLE), INTENT(INOUT) :: PT(*)
Which is coherent with the example, but not with the documentation.
Can you please clarify this aspect?
P.S.: since 32 bit architectures aren't supported anymore, you should probably remove
INTEGER for 32-bit or 64-bit architectures
