topic Simple code gets error=-1 inconsistent input in Pardiso in Intel® oneAPI Math Kernel Library

Simple code gets error=-1 inconsistent input in Pardiso

RichardY1 — Thu, 13 Apr 2023 06:02:47 GMT

Hello,

I'm new on using intel MKL Pardiso solver (2023.1.0 version) to solve for a small matrix as a demo on CentOS7. On calling pardiso(), output error value shows to be -1 for both phase 12 and phase 33, which said to be inconsistent input in documentation. That's the only info I can find throughout the Internet. I have no way to debug this issue. Could you help with this? Thanks a lot. Please check the code below.

#include "mkl.h"

int main (void)

{

MKL_INT n = 5;

MKL_INT ia[ ] = {1, 2, 4, 5, 6, 7};

MKL_INT ja[ ] = {1, 2, 4, 3, 4, 5};

double a[ ] = {1, 1, 2, 1, 1, 1};

MKL_INT mtype = 11;

void *pt[64];

double b[5], x[5];

MKL_INT nrhs = 1;

MKL_INT iparm[64];

MKL_INT maxfct, mnum, phase, error, msglvl;

MKL_INT i;

double ddum; /* Double dummy */

MKL_INT idum; /* Integer dummy. */

for ( i = 0; i < 64; i++ ){

iparm[i] = 0;

pt[i] = 0;

}

iparm[0] = 1; /* No solver default */

iparm[1] = 2; /* Fill-in reordering from METIS */

iparm[7] = 2; /* Max numbers of iterative refinement steps */

iparm[9] = 13; /* Perturb the pivot elements with 1E-13 */

iparm[10] = 1; /* Use nonsymmetric permutation and scaling MPS */

iparm[11]= 0; /* Transposed Matrix */

iparm[17] = -1; /* Output: Number of nonzeros in the factor LU */

iparm[18] = -1; /* Output: Mflops for LU factorization */

iparm[34] = 0; /* zero-based indexing */

iparm[36] = 0; /* CSR format */

maxfct = 1; /* Maximum number of numerical factorizations. */

mnum = 1; /* Which factorization to use. */

msglvl = 1; /* Print statistical information in file */

error = 0; /* Initialize error flag */

phase = 12;

PARDISO (pt, &maxfct, &mnum, &mtype, &phase,

&n, a, ia, ja, &idum, &nrhs, iparm, &msglvl, &ddum, &ddum, &error);

phase = 33;

b[0] = 1;

b[1] = 10;

b[2] = 3;

b[3] = 4;

b[4] = 5;

PARDISO (pt, &maxfct, &mnum, &mtype, &phase,

&n, a, ia, ja, &idum, &nrhs, iparm, &msglvl, b, x, &error);

printf ("\nThe solution of the system is: ");

for ( i = 0; i < n; i++ ){

printf ("\n x [%d] = % f", i, x[i]);

}

Re: Simple code gets error=-1 inconsistent input in Pardiso

RichardY1 — Fri, 14 Apr 2023 06:22:12 GMT

I found a fix.

There are several things that are related to this issue. First of all, compiling MKL package looks complicated to me. I'm using cmake with ninja to compile my program with gnu compiler. At first I manually specified the library and include path of MKL and referred to this Link method for compiling flags. With these setup, the above simple code snippet reports error=-1.

Then I found this guide that shows example cmake code. Noticing that I should set MKL_INTERFACE to lp64 in CMakeList.txt since I'm using gnu compiler, otherwise the compiler complains that the definition of MKL_INT is long long, which can not be converted to int. Another thing is that I can switch between sequential and parallel with MKL_THREADING variable.

Re: Simple code gets error=-1 inconsistent input in Pardiso

Hsbsh — Thu, 13 Apr 2023 10:16:11 GMT

LLP64 is windows. LP64 is linux. So it is simple as that. You cannot use LLP64 on linux, it is just not how it works.

Re: Simple code gets error=-1 inconsistent input in Pardiso

RichardY1 — Fri, 14 Apr 2023 01:51:49 GMT

I don't think that's correct. In MKLConfig.cmake, it said "Values: lp64, ilp64, GNU or INTEL interface will be selected based on Compiler"

Re: Simple code gets error=-1 inconsistent input in Pardiso

Hsbsh — Fri, 14 Apr 2023 01:55:00 GMT

What is not correct? ILP64 is even more rare, was used on Solaris.

LLP64 will not be there, as I said, that is only Windows.

Re: Simple code gets error=-1 inconsistent input in Pardiso

VarshaS_Intel — Fri, 14 Apr 2023 18:04:20 GMT

Hi,

Thanks for posting in Intel Communities.

We have tried at our end using the g++ and icpx compiler and we are able to get the expected results. Please find the below screenshot:

Could you please let us know the GNU compiler version you are using? Also, could you please provide us with the CMakeLists you are using to investigate more at our end?

Thanks & Regards,

Varsha

Re: Simple code gets error=-1 inconsistent input in Pardiso

RichardY1 — Mon, 17 Apr 2023 08:31:31 GMT

Hi Varsha,

Thank you for your reply. The problem is actually due to my incorrect usage and I put the solution in the first reply. I don't have further questions and you can close this topic.

Regards,

Re:Simple code gets error=-1 inconsistent input in Pardiso

VarshaS_Intel — Tue, 18 Apr 2023 06:17:23 GMT

Hi Richard,

Thanks for the confirmation. Glad to know that your issue is resolved. This thread will no longer be monitored by Intel. If you need additional information, please start a new thread.

Have a Good day!

Thanks & Regards,

Varsha