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
- Software Development SDKs and Libraries
- Intel® oneAPI Math Kernel Library
- Sparse Matrix Structure leads Zero Pivot Error

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

diogoneto

Beginner

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

09-03-2013
10:01 AM

119 Views

Sparse Matrix Structure leads Zero Pivot Error

Hi all,

I am using the DSS from MKL to solve the linear system of equations Ax=b, for which the resulting matrix A is sparse and not symmetric, but its structure is symmetric.

Concerning the options, I am using MKL_DSS_SYMMETRIC_STRUCTURE in the definition of the structure and MKL_DSS_INDEFINITE in the factorization.

I am using the following example (15*15 matrix in attachment) to try and solve a system considering that the matrix A is sparse. Unfortunately, the DSS solver gives an error of “MKL-DSS-DSS-Error, Zero Pivot detected”, but the matrix is non-singular. The determinant of the matrix is about -3.1e+62 (calculated in excel) and I get the solution both in excel and matlab.

When I switch the definition of the matrix structure for MKL_DSS_NON_SYMMETRIC, I get the correct solution. I don’t understand why!

In order to understand this problem I create a new subroutine. Using the initial definition of matrix structure, I give all elements of A as non-zero (dense matrix), I am not take advantage of the sparse matrix. With this structure, DSS solver gives the correct solution.

Someone can help me why I cannot use the MKL_DSS_SYMMETRIC_STRUCTURE definition in the sparse matrix structure.

Both subroutines are in attachment.

Thank you.

Diogo

Link Copied

5 Replies

Gennady_F_Intel

Moderator

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

09-04-2013
01:38 AM

119 Views

Gennady_F_Intel

Moderator

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

09-04-2013
09:51 PM

119 Views

Hi Diogo,

Yes, we confirmed that this is the bug and the cause of the problem has been found. How important this issue is for you? Can you use the MKL_DSS_NON_SYMMETRIC case as a temporarily workaroud? we expect no performance degradation in that case.

-Gennady

diogoneto

Beginner

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

09-05-2013
02:37 AM

119 Views

Hi Gennady,

Thank you for your quick response. I have not problem to use the option MKL_DSS_NON_SYMMETRIC and I check yesterday that the performance is not affected with a simple problem.

Diogo

Gennady_F_Intel

Moderator

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

09-05-2013
08:27 PM

119 Views

Gennady_F_Intel

Moderator

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

10-27-2013
11:44 PM

119 Views

Diogo,

please check the latest update 1 ( MKL 11.1 Update 1) where the fix of the problem has been provided. Pls let us know if the problem is still there.

Topic Options

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

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