I'm trying to solve a symmetric indefinite system with PARDISO but the returned solution is somewhat strange, it's made of only large values (like 3.53e+18). I'm attaching a reproducer, I would really appreciate some insight into this issue.
Thank you in advance.
That is right, but it has a really low rank deficiency of 1. For other test cases with similar spectral properties, I have an actual solution, so I was wondering what is happening here.