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.
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.