I am trying to
change a finite element code that currently uses IMSL routines
(specifically LFTQS and LFSQS) to the corresponding MKL routines (PBTRF
The code will run fine for a few timesteps, but eventually, the outputs become unstable and the code stops running. Everything works fine when I run it with the IMSL library.
Everything I have read on the MKL vs IMSL threads state that the MKL has higher performance. Is there a reason they are less stable? Are there any steps I can try to get the code to run stably with the MKL routines?