I have implemented pardiso solver in our model for flow simulation. Pardiso can solve most of our problem well and produce correct results. But recently I have a case that pardiso generates quite different results with our origional solver (ws209), which has been used for over ten years for our model.I guess if there is something wrong with the setting.
Please find the test codes (pardiso_unsym_f.f90), sparse matrices exported from our model (a_i.txt, b_i.txt, ia_i.txt, ja_i.txt), result generated by pardiso (x_out_i.txt) and result generated by ws209 solver (x_ws209_i.txt). The first value in each file is the number of values.
Can anybody help to have a check?
Thanks and regards,