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.
Thanks for reaching out to us.
This issue is a duplicate of the below
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.