Intel® oneAPI Math Kernel Library
Ask questions and share information with other developers who use Intel® Math Kernel Library.
Announcements
Welcome to the Intel Community. If you get an answer you like, please mark it as an Accepted Solution to help others. Thank you!
6433 Discussions

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

Sergey_S_1
Beginner
91 Views

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

0 Kudos
4 Replies
Gennady_F_Intel
Moderator
91 Views

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

Sergey_S_1
Beginner
91 Views

Hi Gennady,

Thank you for quick response!

Sergey Solovev

Gennady_F_Intel
Moderator
91 Views

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

Sergey_S_1
Beginner
91 Views

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

Reply