I am comparing several solvers, and while some ones give expected results, PARDISO keeps giving very different results each time I run the program, or even different results if I call it two consecutive times.
It happens mainly with large matrices (80000 rows, 1000000 nonzeros) while other tests I made with very small matrices (eg: 5x5) give correct results.
With some settings it even crashes.
I mainly use default options for iparm.
I also tried DSS interface, and I get almost the same behavior.
Any idea ?