I am using (INTEL_MKL) PARDISO solver in a visco-elasto-plastic Stokes flow code. It is a finite difference code which uses marker-in-cell technique and has a mixed formulation with linear velocity and constant pressure.
During the reordering and symbolic factorization (METIS) I keep getting SEG-faults (and occasionally a SIGBUS) at random moments, where it will probably runs through the next time. This occurs both when using METIS = 3 or 2, however it runs fine without re-ordering.
I have checked the physics, matrix conditioning etc. but the fact that solver only has a problem with longer runs, whereas it runs fine with the restart file, suggests that there is some algorithmic/compiler specific issue. I have traced the memory utilization and it is far less than the RAM.
This seems to be the problem with the INTEL_MKL version of Pardiso only.