Intel® oneAPI Math Kernel Library
Ask questions and share information with other developers who use Intel® Math Kernel Library.

## choice of iparms in Pardiso

Beginner
266 Views

hello,

I would like to solve a linear system of type Ax = b with A,B and x in the file joined.

The matrix A is well-conditioned, so I must be wrong with the iparms :

357         iparm[0] = 1; /* No solver default */
358         iparm[1] = 2; /* Fill-in reordering from METIS */
359         iparm[2] = 8;
360         iparm[3] = 0; /* CGS */
361         iparm[4] = 0; /* No user fill-in reducing permutation */
362         iparm[5] = 0; /* Write solution into x */
363         iparm[6] = 0; /* Not in use */
364         iparm[7] = 0; /* Max numbers of iterative refinement steps */
365         iparm[8] = 0; /* Not in use */
366         iparm[9] = 13; /* Perturb the pivot elements with 1E-13 */
367         iparm[10] = 1; /* Use nonsymmetric permutation and scaling MPS */
368         iparm[11] = 0; /* Not in use */
369         iparm[12] = 0; /* Not in use */
370         iparm[13] = 0; /* Output: Number of perturbed pivots */
371         iparm[14] = 0; /* Not in use */
372         iparm[15] = 0; /* Not in use */
373         iparm[16] = 0; /* Not in use */
374         iparm[17] = -1; /* Output: Number of nonzeros in the factor LU */
375         iparm[18] = -1; /* Output: Mflops for LU factorization */
376         iparm[19] = 0; /* Output: Numbers of CG Iterations */
377         iparm[27] = 1; /* check the data structure */
378         iparm[31] = 1; /* iterative solver*/
380         maxfct = 1; /* Maximum number of numerical factorizations. */
381         mnum = 1; /* Which factorization to use. */

The Package ID of mkl is : l_mkl_p_10.0.011

Antoine

4 Replies
Black Belt
266 Views
In producing the attachment, you seem to have used a version of tar whose output my Cygwin version cannot process: