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!
6427 Discussions

PARDISO error for structurally symmetric matrix

asd__asdqwe
Beginner
79 Views

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.

0 Kudos
5 Replies
Gennady_F_Intel
Moderator
79 Views

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

Gennady_F_Intel
Moderator
79 Views

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.

asd__asdqwe
Beginner
79 Views

Thank you.

Gennady_F_Intel
Moderator
79 Views

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.

MariaZh
Employee
79 Views

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

Reply