Turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- Intel Community
- Software Development SDKs and Libraries
- Intel® oneAPI Math Kernel Library & Intel® Math Kernel Library
- PARDISO error for structurally symmetric matrix

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Mute
- Printer Friendly Page

Highlighted
##

asd__asdqwe

Beginner

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

02-23-2016
09:20 AM

18 Views

PARDISO error for structurally symmetric matrix

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.

5 Replies

Highlighted
##

Gennady_F_Intel

Moderator

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

02-23-2016
08:18 PM

18 Views

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

Highlighted
##

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.

Gennady_F_Intel

Moderator

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

02-23-2016
10:27 PM

18 Views

Highlighted
##

asd__asdqwe

Beginner

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

02-24-2016
05:43 AM

18 Views

Thank you.

Highlighted
##

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.

Gennady_F_Intel

Moderator

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

02-25-2016
12:25 AM

18 Views

Highlighted
##

MariaZh

Employee

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

02-26-2016
09:02 AM

18 Views

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

For more complete information about compiler optimizations, see our Optimization Notice.