Community
cancel
Showing results for 
Search instead for 
Did you mean: 
Iosif_M_
Beginner
72 Views

Sparse eigensolver

Jump to solution

Dear colleagues, is it possible to find eigenvalues and eigenvectors by means of MKL routines in the following case: Ax = kx, where A is a sparse, general, complex matrix. I found the routines for sparse symmetric/hermitian matrices but not for general sparse matrices.

0 Kudos
1 Solution
Gennady_F_Intel
Moderator
72 Views

Iosif, that's not possible for general sparse cases for now. one of the existing solution is feast 3.0 which is on our radar but without some timeline at this moment. Do you have some specific timeline when you need to have such functionality available?

View solution in original post

5 Replies
Ying_H_Intel
Employee
72 Views

Hi Losif, 

For general sparse matrix, The Pardiso or DSS solver may be you wanted. 

Best Regards,

Ying 

 

Iosif_M_
Beginner
72 Views

Dear Ying,

Thank you for the response. As far as I understand PARDISO can solve sparse SLAE, but it is not eigensolver.

Best regards,

Iosif

Ying_H_Intel
Employee
72 Views

Hi Losif, 

then other option may be : 

The Extended Eigensolver functionality offers:
• Real/Complex and Single/Double precisions: double precision is recommended to provide better accuracy
of eigenpairs.
• Reverse communication interfaces (RCI) provide maximum flexibility for specific applications. RCI are
independent of matrix format and inner system solvers, so you must provide your own linear system
solvers (direct or iterative) and matrix-matrix multiply routines.
• Predefined driver interfaces for dense, LAPACK banded, and sparse (CSR) formats are less flexible but are
optimized and easy to use:
• The Extended Eigensolver interfaces for dense matrices are likely to be slower than the comparable
LAPACK routines because the FEAST algorithm has a higher computational cost.
7 Intel® Math Kernel Library Reference Manual
1628
• The Extended Eigensolver interfaces for banded matrices support banded LAPACK-type storage.
• The Extended Eigensolver sparse interfaces support compressed sparse row format and use the Intel
MKL PARDISO solver.

Best wishes, 

Ying

Gennady_F_Intel
Moderator
73 Views

Iosif, that's not possible for general sparse cases for now. one of the existing solution is feast 3.0 which is on our radar but without some timeline at this moment. Do you have some specific timeline when you need to have such functionality available?

View solution in original post

Iosif_M_
Beginner
72 Views

Dear Gennady,

I have also found that functionality in the FEAST solver. Just now it is enough for us to use FEAST and switch to MKL/FEAST later.

Iosif

Reply