When I do a
call pardiso_getdiag( pt, df, da, mnum, error)
the arrays DF and DA are filled with 'factored' and 'original' diagonal pivots.
The pivots are in the same 'order' of the 'original' matrix or they are 'permuted'?
Both arrays returned for permuted matrix.
For more complete information about compiler optimizations, see our Optimization Notice.