Sparse Schur Complement is Computed Wrong

[Havlu]

Hi,

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

[ArpitaP_Intel]

Hi Abdullah,

Thanks for reaching out to us.

This issue is a duplicate of the below

https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Sparse-Schur-Complement-is-Computed-Wrong/m-p/1297389#M31723

We are working on your issue and you can track your issue using the above link .

We are closing this thread and would no longer monitor it.

Thanks!