topic Hi Gennady, in Intel® oneAPI Math Kernel Library

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

Sergey_S_1 — Fri, 18 Jan 2013 12:27:52 GMT

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

Hi Sergey,

Gennady_F_Intel — Fri, 18 Jan 2013 13:04:12 GMT

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,

Sergey_S_1 — Mon, 21 Jan 2013 06:45:23 GMT

Hi Gennady,

Thank you for quick response!

Sergey Solovev

Sergey, pls check the updated

Gennady_F_Intel — Sat, 06 Apr 2013 20:44:03 GMT

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

Gennady, I have checked the

Sergey_S_1 — Mon, 15 Apr 2013 06:06:05 GMT

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