topic Thank you Gennady, I will in Intel® oneAPI Math Kernel Library

Linear Systems and Unbalanced Matrix

Gianluca_G_ — Tue, 01 Mar 2016 07:31:46 GMT

Hello,

We are using this LAPACKE_dgesv (MKL function) in order to solve a linear system with dense coefficient matrix.

Sometimes the matrix is unbalanced (too different values), so bad conditioned! in this case results can be wrong.

Is there any mkl function that can alert about this situation, before running the solver?

Thank you

Gianluca

Guinluca,

Gennady_F_Intel — Wed, 02 Mar 2016 05:54:51 GMT

Guinluca,

there are a number of lapack's outines for estimating the condition number of a matrix. please have a look at the gecon() as an example. You may check the cond number of your case in advance solving the system.

or you may try to play with gerfs() routine which performs an iterative refinement of the solution to a system of linear equations

-- Gennady

Thank you Gennady, I will

Gianluca_G_ — Thu, 03 Mar 2016 07:08:13 GMT

Thank you Gennady, I will check it

Gianluca