Solving complex-valued linear systems via equivalent real formulations

被引:95
作者
Day, D [1 ]
Heroux, MA [1 ]
机构
[1] Sandia Natl Labs, Dept Appl Math, Massively Parallel Comp Res Lab, Albuquerque, NM 87185 USA
关键词
complex linear systems; iterative methods; sparse matrices; preconditioning;
D O I
10.1137/S1064827500372262
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
Most preconditioned iterative methods apply to both real- and complex-valued linear systems. At the same time, most iterative linear solver packages available today focus exclusively on real- valued systems or deal with complex-valued systems as an afterthought. By recasting the complex problem in a real formulation, a real- valued solver can be applied to the equivalent real system. On one hand, real formulations have been dismissed due to their unfavorable spectral properties. On the other hand, using an equivalent preconditioned real formulation can be very effective. We give theoretical and experimental evidence that an equivalent real formulation is useful in a number of practical situations. Furthermore, we show how to use the advanced features of modern solver packages to formulate equivalent real preconditioners that are computationally efficient and mathematically identical to their complex counterparts. The effectiveness of equivalent real formulations is demonstrated by solving ill-conditioned complex-valued linear systems for a variety of large scale applications. Moreover, the circumstances under which certain equivalent real formulations are competitive is more clearly delineated.
引用
收藏
页码:480 / 498
页数:19
相关论文
共 16 条
[1]  
[Anonymous], 1999, OFFICIAL AZTEC USERS
[2]  
Balay S, 1997, MODERN SOFTWARE TOOLS FOR SCIENTIFIC COMPUTING, P163
[3]  
Balay S., 1998, ANL9511
[4]  
Barrett R., 1994, Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods, V2nd ed.
[5]  
Fischer B, 1996, Tech. rep.
[6]  
FRAYSEE V, 1997, TRPA9749 CERFACS
[7]   QMRPACK: A package of QMR algorithms [J].
Freund, RW ;
Nachtigal, NM .
ACM TRANSACTIONS ON MATHEMATICAL SOFTWARE, 1996, 22 (01) :46-77
[8]   CONJUGATE GRADIENT-TYPE METHODS FOR LINEAR-SYSTEMS WITH COMPLEX SYMMETRICAL COEFFICIENT MATRICES [J].
FREUND, RW .
SIAM JOURNAL ON SCIENTIFIC AND STATISTICAL COMPUTING, 1992, 13 (01) :425-448
[9]  
Golub G.H., 2013, MATRIX COMPUTATIONS
[10]  
HEROUX MA, 2000, KOMPLEX SOLVER PACKA