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.
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]);
}
