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

Sparse Iteartive solvers with Complex numbers

Dogan__Hakan
Beginner
542 Views
Hi all,

I am trying to use iterative solvers for sparse systems, such as subroutines ILUT for preconditioning ang DFGMRES for the iteration but apparently these work with DOUBLE PRECISION but not COMPLEX numbers.

Do you have any suggestions how can proceed further?
Thank you very much in advance.

Hakan
0 Kudos
4 Replies
Gennady_F_Intel
Moderator
542 Views
Hi Hakan,
Currently MKL doesnt support RSI solvers for complex data. Our suggestion - you may try to use direct solver for complex data types.
--Gennady
0 Kudos
Dogan__Hakan
Beginner
542 Views
Hi Gennady,

Well actually my program works with PARDISO, but I need to solve larger scale problems and I thought iterative solvers can handle more. Am I true with this thought? And finally would you recommend 'SPARSKIT' library to import into the Intel Fortran Complier?

Thank you very much.

Hakan
0 Kudos
Victor_K_Intel1
Employee
542 Views
Hakan,

Formally speaking, support of complex data is not crucial. For example, if you need to solve a system of linear equations Ax=f with complex matrix A=B+i*C and complex right hand side vector f=g+i*h assuming the solutionhas the form x=y+i*z you can get a system with real coefficient matrix and right hand side:
(B -C) (y) =(g)
(C B) (z) (h)
The system is nothing esle but equations you obtain when equate separately real and imaginary part of the original system. The coefficient matrix orderand the rhs vector length are doubled. This approach might benotthe most efficient from the performance point of view and requires more memory but this is another story.

With best regards,
Victor

0 Kudos
Dogan__Hakan
Beginner
542 Views
Hello again,

Thank you for the answer, this will help I think. But before trying that I want to ask something. I have got the latest version of IMSL library which has some improved mathematical functions. Do you reckon usage of IMSL would make it possible to use GMRES with complex numbers?

Thank you in advance,
Hakan
0 Kudos
Reply