topic I can now extract and view in Intel® oneAPI Math Kernel Library

choice of iparms in Pardiso

Antoine__A_ — Thu, 29 Nov 2012 14:57:11 GMT

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

Thanks for any help you can give me
Antoine

In producing the attachment,

mecej4 — Fri, 30 Nov 2012 14:55:00 GMT

In producing the attachment, you seem to have used a version of tar whose output my Cygwin version cannot process:

tar: Ignoring unknown extended header keyword `SCHILY.dev' tar: Ignoring unknown extended header keyword `SCHILY.ino' tar: Ignoring unknown extended header keyword `SCHILY.nlink'

ok, you can try this

Antoine__A_ — Fri, 30 Nov 2012 15:28:19 GMT

ok, you can try this

I can now extract and view

mecej4 — Fri, 30 Nov 2012 17:59:09 GMT

I can now extract and view the files, thanks. However, it is not clear what the files contain. Does 36_A.mat contain the upper triangle of a symmetric matrix? Or does it contain a triangular matrix? Do you expect the matrix to be positive-definite? What is the reason behind giving values for the vector x? Your matrix appears to be 252 X 252, but the files for x and b have 1296 lines each. How do you explain this?

Sorry, you are right, I

Antoine__A_ — Mon, 03 Dec 2012 10:38:00 GMT

Sorry, you are right, I should have specified my problem. The matrix I sent you was a 1296 x 1296 matrix in COO format. I send it to you in CSR format, it will be easier. The new file for A is composed of 4 lines : the size of data (row column val), the ranks of the first element of the row, the columns, and the values. This matrix has no particularity : non-symmetrical, not positive-definite... I give 2 vectors x here : one obtained with Matlab (very good acuracy), and one calculated by Pardiso. Maybe you need it to compare, or just seeing how wrong is Pardiso with my parameters.