topic Re: pardiso documentation inconsistent with example in Intel® oneAPI Math Kernel Library

pardiso documentation inconsistent with example

Riccardo — Thu, 05 Dec 2024 13:47:12 GMT

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

Re:pardiso documentation inconsistent with example

Aleksandra_K — Tue, 10 Dec 2024 10:55:06 GMT

Dear Riccardo,

Thank you for noting the problem in the documentation. I have escalated the issue to have it corrected.

Best regards,

Aleksandra

Re: pardiso documentation inconsistent with example

Riccardo — Thu, 12 Dec 2024 08:07:19 GMT

Thanks a lot Aleksandra! Please write the reply here when you get it.

Best regards,

Riccardo

Re:pardiso documentation inconsistent with example

Aleksandra_K — Mon, 16 Dec 2024 16:26:04 GMT

Dear Riccardo,

I wanted to let you know that we will track the documentation issue in our internal systems. For now, I will be closing this ticket.

If you need anything else, please feel free to contact us.

Best regards,

Aleksandra