Intel® oneAPI Math Kernel Library
Ask questions and share information with other developers who use Intel® Math Kernel Library.
7287 Discussions

Comparison of MKL DSS with PARDISO and MATLAB for indefinite unsymmetric matrix

mchoi
Novice
635 Views

Hello, 

I have compared the solutions from using MKL DSS, MKL PARDISO, and MATLAB. This was motivated by the fact that DSS gives seemingly inaccurate solution when the matrix is badly conditioned, compared with PARDISO and MATLAB.  In the attached files, I have two examples of unsymmetric indefinite (system) matrices:

 

1. input_matrix_ex1: small matrix

 

2. input_matrix_ex2: larger matrix (badly conditioned)

 

In the second example, the DSS solution is significantly different from those using MKL PARDISO and MATLAB, see the figures below.

 

My question is what makes the difference of solutions between DSS and PARDISO? 

 

comparison_dss_pardiso_matlab.pngcomparison_dss_pardiso_matlab_magnified.png

 

Thank you.

 

0 Kudos
0 Replies
Reply