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

pardiso_getdiag

Alemdar__Bulent
Beginner
390 Views

Hello

I am calling "pardiso_getdiag( pt, df, da, mnum, error)" after PARDISO call (phase = 22) and it returns diagonal terms of initial and final (factorized) matrix.

As noticed, these results returned are not in the same order of the original matrix (they are permuted).

 

Is there a way to get these results in the same order of original matrix.

 

Bulent

 

 

 

0 Kudos
1 Reply
Alexander_K_Intel2
390 Views

Hi,

Currently PARDISO doesn't support ability to reorder return vector to achieve initial one. That's come from the fact that for non-symmetric matrix we use non-symmetric permutation and it is not clear what is order of diagonal elements in such case.

Thanks,

Alex

0 Kudos
Reply