Showing results for

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

xian-zhong_guous_cd-

Beginner

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

10-06-2014
12:31 PM

56 Views

mkl metis SIGSEGV

We have a symmetric matrix (attached). MKL crashes on metis(see the attached callstack). Can you take a look?

symmetric_matrix.txt:

first line: number of eqautions (n)

second line: index based (0 means 0 based)

third line: number of nonzeros (nz)

rest: non zeros in triplet format (r c a)

Thanks,

Link Copied

3 Replies

xian-zhong_guous_cd-

Beginner

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

10-06-2014
12:34 PM

56 Views

Major version: 11

Minor version: 1

Update version: 0

Product status: Product

Build: n20130711

mecej4

Black Belt

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

10-06-2014
02:40 PM

56 Views

The outcome is not surprising, because the matrix that you have given does not appear to be symmetric as claimed. The convention is that you pass only the upper triangle of a symmetric matrix to MKL/Pardiso, in the CSR format. The MKL manual says "The Intel MKL direct sparse solvers use a row-major upper triangular storage format: the matrix is compressed row-by-row and for symmetric matrices only non-zero elements in the upper triangular half of the matrix are stored".

Your data file is in the COO format, and line 5282 contains '739 733 -0.0202558', which implies that row 740 of the matrix starts in column 734, which would be located in the lower triangle. The Pardiso routine is, probably, unable to address this discrepancy. There are many other entries (over 50) in the lower triangle. You have to resolve this inconsistency before trying to reorder/factorize the matrix using MKL/Pardiso. I have done a bit of speculation with the comments here since you did not state how you did the COO to CSR conversion nor did you provide information on the arguments that you called Pardiso with.

xian-zhong_guous_cd-

Beginner

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

10-06-2014
02:45 PM

56 Views

Thanks for the quick response.

- 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.