Hello,
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.