A sparse approximate inverse preconditioner for nonsymmetric linear systems

被引:226
作者
Benzi, M
Tuma, M
机构
[1] Los Alamos Natl Lab, Sci Comp Grp CIC19, Los Alamos, NM 87545 USA
[2] Acad Sci Czech Republic, Inst Comp Sci, Prague 18207 8, Czech Republic
关键词
preconditioning; approximate inverses; sparse linear systems; sparse matrices; incomplete factorizations; conjugate gradient-type methods;
D O I
10.1137/S1064827595294691
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
This paper is concerned with a new approach to preconditioning for large, sparse linear systems. A procedure for computing an incomplete factorization of the inverse of a nonsymmetric matrix is developed, and the resulting factorized sparse approximate inverse is used as an explicit preconditioner for conjugate gradient-type methods. Some theoretical properties of the preconditioner are discussed, and numerical experiments on test matrices from the Harwell-Boeing collection and from Tim Davis's collection are presented. Our results indicate that the new preconditioner is cheaper to construct than other approximate inverse preconditioners. Furthermore, the new technique insures convergence rates of the preconditioned iteration which are comparable with those obtained with standard implicit preconditioners.
引用
收藏
页码:968 / 994
页数:27
相关论文
共 57 条
[1]  
ANDERSON EC, 1988, 805 CSRD U ILL URB C
[2]  
Axelsson O., 1994, ITERATIVE SOLUTION M
[3]  
Barrett R., 1994, Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods, V2nd ed.
[4]   PARALLEL ALGORITHMS FOR THE SOLUTION OF CERTAIN LARGE SPARSE LINEAR-SYSTEMS [J].
BENSON, M ;
KRETTMANN, J ;
WRIGHT, M .
INTERNATIONAL JOURNAL OF COMPUTER MATHEMATICS, 1984, 16 (3-4) :245-260
[5]  
Benzi, 1993, THESIS N CAROLINA ST
[6]   A DIRECT PROJECTION METHOD FOR SPARSE LINEAR-SYSTEMS [J].
BENZI, M ;
MEYER, CD .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 1995, 16 (05) :1159-1176
[7]   A sparse approximate inverse preconditioner for the conjugate gradient method [J].
Benzi, M ;
Meyer, CD ;
Tuma, M .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 1996, 17 (05) :1135-1149
[8]  
BENZI M, 1994, SIAM PROC S, P294
[9]  
BENZI M, 1995, 653 CZECH AC SCI I C
[10]   KRYLOV METHODS PRECONDITIONED WITH INCOMPLETELY FACTORED MATRICES ON THE CM-2 [J].
BERRYMAN, H ;
SALTZ, J ;
GROPP, W ;
MIRCHANDANEY, R .
JOURNAL OF PARALLEL AND DISTRIBUTED COMPUTING, 1990, 8 (02) :186-190