Intel® oneAPI Math Kernel Library
Ask questions and share information with other developers who use Intel® Math Kernel Library.
Announcements
This community is designed for sharing of public information. Please do not share Intel or third-party confidential information here.
6559 Discussions

Sparse Schur Complement is Computed Wrong

Havlu
Novice
1,173 Views

Hi,

[Apologies for cross-posting, I don't know which board would be proper to use]

I am trying to use pardiso to compute the Schur complement of a matrix. I had some stability issues for some large matrices, so I decided to put in a tiny matrix and see if I can figure out the issue. However, it looks like for the small matrix I have, the Schur complement is computed wrong.

This matrix is real symmetric and 8-by-8. The code is attached, and it is commented for ease of reading. Any help is appreciated.

Thank you,
Abdullah

Labels (3)
0 Kudos
8 Replies
MRajesh_intel
Moderator
1,123 Views

Hi,

 

Thanks for posting your query. I have executed your sample code and the results are attached below. If these are not the right values, can you please provide us the right results? Can you please also share the MKL version used, OS version used?

 

Thanks

Rajesh.

 

Havlu
Novice
1,116 Views

Hi Rajesh,

 

Thank you for the reply. The matrix outputted is wrong, and to my knowledge, it is not even a valid CSR3 matrix.

Here is the correct Schur complement to 3 digits after the point:

// [ 1.667 0.167 -0.204]
// [ 0.167 1.667 0.204]
// [-0.204 0.204 1.500]

You can also find this in the source code, it is provided there for convenience. Since it is symmetric, according to the documentation, Pardiso should return the lower triangular part of this matrix in the CSR3 format.

I am on Fedora Linux, and I am using MKL 2021.3.0.

Best,
Abdullah

Khang_N_Intel
Employee
1,083 Views

Hi Abdullah,

Just want to let you know that I am discussing with the engineering team right now and will let you know what the next step is.

Khang


Havlu
Novice
1,055 Views

Hi Khang,

Thanks for the update. I am looking forward to hearing from you.

Best,
Abdullah

Havlu
Novice
1,028 Views

Hi Khang,

Just checking in. Are there any updates, or news on the next steps?

Best,
Abdullah

Khang_N_Intel
Employee
971 Views

Hi Abdullah,


The issue has been confirmed. We have already escalated it to the engineering team. We will let you know when we receive any information from the engineering team.


Best,

Khang



Havlu
Novice
377 Views

Hey Khang,

 

Hope all is well. Are there any updates on this matter? I was hoping to see a bugfix in the releases since, but I haven't seen anything regarding this issue and it seems to persist in the latest version.

 

Best,

Abdullah

Gennady_F_Intel
Moderator
269 Views

Abdullah,

there are no updates on this issue at this moment. we will keep this thread informed as soon as any updates.

-Gennady


Reply