PARDISO reordering: incorrect description perm() array in the MKL manual

PARDISO reordering: incorrect description perm() array in the MKL manual

Hi! I want to use own reordering in PARDISO. From manual I see following description of perm() array:  “Let A be the original matrix and C = P*A*PT be the permuted matrix. Row (column) i of A is the perm(i) row (column) of C.” Unfortunately, I got incorrect results, i.e. incorrect number of nonzero elements in L factors. After investigation, I found out that description of perm() array is incorrect. User should to set inversion of permutation instead permutation. I.e. “Row (column) i of C is the perm(i) row (column) of A.” Could you clarify, I am right or not?  

P.S. Test case has been attached.

Sergey Solovev

AttachmentSize
Downloadtext/x-c++src pds-reor-bug.cpp8.81 KB
5 posts / 0 new
Last post
For more complete information about compiler optimizations, see our Optimization Notice.

Hi Sergey, 

the problem is reproduced and at the first glance it looks like the documentation's problem - the description should be slightly changed accordingly. we will check the behavior again and let you know.

--Gennady

Hi Gennady,

Thank you for quick response!

Sergey Solovev

Sergey, pls check the updated doc in the latest 11.0 update3. regards, Gennady

Gennady, I have checked the fix. All is OK, thank you and MKL support team!

Leave a Comment

Please sign in to add a comment. Not a member? Join today