I am writing to you in order to inform you about an unusual behavior that pardiso solver exhibit in the latest version of mkl library. More specifically, we use this solver in our finite element code in order to solve the linear system of equations. Up to version ..., we could solve systems that arise when periodic boundary conditions are imposed. Unfortunately, the latest version can not solve them and a NAN appears during back- substitution. We have used both intel and gnu compilers as well as multithreaded and sequential versions without any improvement. Finally, these can be solved with WATSON and HSL solvers.
I modified the example pardiso_unsym_f.f from the MKL 11.0.1 distribution to read your data files, running under Windows 7 X64 using the 13.0.1.119 compilers (32-bit and 64-bit). Your system was solved successfully and there were no NaNs in the solution vector.
Please provide code that reads your data files and produces the behavior that you described (NaNs in the solution).
PARDISO PROBLEM
Dear developer:
I am writing to you in order to inform you about an unusual behavior that pardiso solver exhibit in the latest version of mkl library. More specifically, we use this solver in our finite element code in order to solve the linear system of equations. Up to version ..., we could solve systems that arise when periodic boundary conditions are imposed. Unfortunately, the latest version can not solve them and a NAN appears during back- substitution. We have used both intel and gnu compilers as well as multithreaded and sequential versions without any improvement. Finally, these can be solved with WATSON and HSL solvers.
I am looking forward for your response.
- Dimitrios Fraggedakis