topic Hi Kirill, in Intel® oneAPI Math Kernel Library

Pardiso Hermitian matrix input

mullervki — Mon, 18 May 2020 18:56:38 GMT

Hello,

If I have a Hermitian matrix (mtype=4 or -4), does it matter whether I enter the upper or lower matrices? For example, would both

int ia[4] = {1,3,5,6};

int ja[5] = {

1, 3,

2, 3,

3};

and

int ia[4] = {1,2,3,6};

int ja[5] = {

1, 2, 3};

be equally valid? I ask because in the second case I'm getting an error:

*** Error in PARDISO (incorrect input matrix ) error_num= 24
*** Input check: j=3, ja(j)=3, ja(j+1)=1 are incompatible

But I can't understand what's wrong with the second input.

Thanks.

Kirill_V_Intel — Mon, 18 May 2020 19:06:33 GMT

Hi,

You need to pass the upper triangular part only.

Best,
Kirill

mullervki — Mon, 18 May 2020 21:49:35 GMT

Hi Kirill,

Thanks! I looked throughout the documentation but missed the page you pointed me to. Much appreciated.

-Arthur