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

I have a 3x4 matrix

0 1 0 0

0 0 2 0

0 0 0 3

the corresponding CSR format is:

ia= [0,1,2,3]

ja= [1,2,3]

va= [1,2,3]

However, when I check it with sparse_matrix_checker, it says:

*** Input check: i=2, ja[i]=3, neqns=3 are incompatible

Matrix check details: (2, 3, 3)

Matrix check result: MKL_SPARSE_CHECKER_OUT_OF_RANGE

Is there any special request for MKL's sparse matrix that are not listed in the reference ?

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

The sparse matrix checker is part of the sparse solvers system APIs. the matrix handle there only has a single dimension size meaning it is assuming that the matrix is square (sparse solvers in pardiso depend on matrix being square even though QR can be rectangular, it is a separate set of algorithms and APIs), so the sparse_matrix_checker assumes your matrix is square.

In your case, you have provided neqns = 3 so it is assuming you have a 3x3 matrix , but you have a rectangular matrix input with a column index which is 3 (base 0) which is out of range, hence the MKL_SPARSE_CHECKER_OUT_OF_RANGE return value.

Spencer

Link Copied

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

The sparse matrix checker is part of the sparse solvers system APIs. the matrix handle there only has a single dimension size meaning it is assuming that the matrix is square (sparse solvers in pardiso depend on matrix being square even though QR can be rectangular, it is a separate set of algorithms and APIs), so the sparse_matrix_checker assumes your matrix is square.

In your case, you have provided neqns = 3 so it is assuming you have a 3x3 matrix , but you have a rectangular matrix input with a column index which is 3 (base 0) which is out of range, hence the MKL_SPARSE_CHECKER_OUT_OF_RANGE return value.

Spencer

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