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

## Low Accuracy on dgeevx results, how to increase accuracy?

Beginner
379 Views

Dear all,

I am trying to solve a non-symmetric tridiagonal eigenvalue problem with dgeevx from MKL, however the accuracy of the eigenvalues are not that good. I am also using dsyevr for symmetric tridiagonal matrices and for that routine there is a tolerance setting. I was wondering if there is a way to increase the accuracy of dgeevx. At the moment, I am solving without balancing and scaling.

Best regards,

Umut

2 Replies
Honored Contributor III
379 Views

On what basis do you conclude that DGEEVX is of low accuracy? Is this experienced with any tridiagonal matrix, or can you give a specific example? What is the size of the matrix?

Have you considered using specialized software, e.g., see the link to Dario Bini's Eigensolve at the end of the page http://fibonacci.dm.unipi.it/~bini/Ricerca/ric.html ?

Beginner
379 Views

mecej4 wrote:

On what basis do you conclude that DGEEVX is of low accuracy? Is this experienced with any tridiagonal matrix, or can you give a specific example? What is the size of the matrix?

Have you considered using specialized software, e.g., see the link to Dario Bini's Eigensolve at the end of the page http://fibonacci.dm.unipi.it/~bini/Ricerca/ric.html ?

I am comparing the results with the results of MATLAB, eig , for accuracy.

These are the reduced problems of a Lanczos solver and sizes range from 40-50 to a value around 500(if the block  sizes are large )

Yes I can give an example, but need more time, I will soon send a matrix where the accuracy seems to be low.

Thx.

Umut