topic have you tried the latest 11 in Intel® oneAPI Math Kernel Library

PARDISO error for structurally symmetric matrix

asd__asdqwe — Tue, 23 Feb 2016 17:20:41 GMT

Hello,

I'm trying to solve a structurally symmetric problem, but for some reason, with the attached program, everything fails except if mtype = 11. Could you help me understand why please ? I've exported the matrix to MATLAB and max(sort(I) - sort(J)) = 0, so the matrix should be indeed structurally symmetric.

Thank you in advance.

have you tried the latest 11

Gennady_F_Intel — Wed, 24 Feb 2016 04:18:18 GMT

have you tried the latest 11.3 update 2 which we received the last week?

I see the problem on the

Gennady_F_Intel — Wed, 24 Feb 2016 06:27:15 GMT

I see the problem on the reordering phase with the latest 11.3 u2 version of MKL. We will investigate the cause of the issue and keep you informed with the status.

Thank you.

asd__asdqwe — Wed, 24 Feb 2016 13:43:22 GMT

Thank you.

MKL doesn’t check what the

Gennady_F_Intel — Thu, 25 Feb 2016 08:25:57 GMT

MKL doesn’t check what the exact type of input matrix is. It might be the good feature request for the future implementation, but we don’t do that right now. For this partially case, it seems your input matrix is not real and structurally symmetric but this is nonsymmetric. You may try to check it by himself if the input matrix is nonsymmetric. You may try to use mkl_dcsradd ( ….. ) routine with the request == 1 and when the second matrix will be the same but transposed. Then you need to check the number of nonzeroes into resulted matrix C. In the case if this nnz will be different compare with original matrix, that’s mean the input matrix is nonsymmetric.

Hello,

MariaZh — Fri, 26 Feb 2016 17:02:44 GMT

Hello,

During investigation of your case I found out that the matrix you've provided doesn't contain all diagonal elements.
According to documentation all diagonal elements must be stored for structurally symmetric matrix (please find description of ja array in https://software.intel.com/ru-ru/node/521678). You can set them to zero and put in values array to avoid problem on reordering phase.

Best regards,
Maria