Quantcast
Channel: Intel® Software - Intel® oneAPI Math Kernel Library & Intel® Math Kernel Library
Viewing all articles
Browse latest Browse all 3005

Pardiso unable to solve a symmetric matrix

$
0
0

Hello MKL tech support,

I have a positive definite matrix and using Pardiso to solve the problem. I kept receiving the error -4 for the solver. It is interative refinement problem in phase 33.

My params for pardiso are:

        iparm[7] = 15;
        iparm[10] = 0;
        iparm[12] = 0;
        iparm[0] = 1;
    iparm[1] = 0;
    iparm[3] = 0;
    iparm[4] = 0;
    iparm[5] = 0;
    iparm[6] = 0;
    iparm[9] = 13;
    
    iparm[11] = 0;
    iparm[13] = 0;

    iparm[17] = -1;
    iparm[18] = 0;
    iparm[19] = 0;
    iparm[20] = 0;
 
    iparm[26] = 1;
    iparm[27] = 0;
    iparm[34] = 1;
    iparm[59] = 0;

I am trying to solve Ax = y but there was no error in the factorization phase until the phase 33.

Thanks,

AttachmentSize
Downloadpile-ill-cond.zip24.08 MB

Viewing all articles
Browse latest Browse all 3005

Trending Articles



<script src="https://jsc.adskeeper.com/r/s/rssing.com.1596347.js" async> </script>