- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content

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

Link Copied

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content

**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'

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page