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

AdjuntoTamaño
Descargar pds-reor-bug.cpp8.81 KB
publicaciones de 5 / 0 nuevos
Último envío
Para obtener más información sobre las optimizaciones del compilador, consulte el aviso sobre la optimización.

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!

Inicie sesión para dejar un comentario.