A sparse approximate inverse preconditioner for the conjugate gradient method

被引:279
作者
Benzi, M
Meyer, CD
Tuma, M
机构
[1] CERFACS, F-31057 TOULOUSE, FRANCE
[2] N CAROLINA STATE UNIV, DEPT MATH, RALEIGH, NC 27695 USA
[3] ACAD SCI CZECH REPUBL, INST COMP SCI, PRAGUE 18207 8, CZECH REPUBLIC
关键词
sparse approximate inverses; preconditioned conjugate gradient method; H-matrices; incomplete factorizations;
D O I
10.1137/S1064827594271421
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
A method for computing a sparse incomplete factorization of the inverse of a symmetric positive definite matrix A is developed, and the resulting factorized sparse approximate inverse is used as an explicit preconditioner for conjugate gradient calculations. It is proved that in exact arithmetic the preconditioner is well defined if A is an H-matrix. The results of numerical experiments are presented.
引用
收藏
页码:1135 / 1149
页数:15
相关论文
共 29 条
[1]  
ALVARADO FL, 1994, GRAPH THEORY SPA IMA, V56, P141
[2]  
Axelsson O., 1994, ITERATIVE SOLUTION M
[3]  
Benzi, 1993, THESIS N CAROLINA ST
[4]   A DIRECT PROJECTION METHOD FOR SPARSE LINEAR-SYSTEMS [J].
BENZI, M ;
MEYER, CD .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 1995, 16 (05) :1159-1176
[5]  
BENZI M, 1994, SIAM PROC S, P294
[6]  
BENZI M, 1995, IN PRESS P 2 IMACS I
[7]  
Chow E, 1994, P COL C IT METH APR
[8]   BLOCK PRECONDITIONING FOR THE CONJUGATE-GRADIENT METHOD [J].
CONCUS, P ;
GOLUB, GH ;
MEURANT, G .
SIAM JOURNAL ON SCIENTIFIC AND STATISTICAL COMPUTING, 1985, 6 (01) :220-252
[9]   APPROXIMATE INVERSE PRECONDITIONINGS FOR SPARSE LINEAR-SYSTEMS [J].
COSGROVE, JDF ;
DIAZ, JC ;
GRIEWANK, A .
INTERNATIONAL JOURNAL OF COMPUTER MATHEMATICS, 1992, 44 (1-4) :91-110
[10]   DECAY-RATES FOR INVERSES OF BAND MATRICES [J].
DEMKO, S ;
MOSS, WF ;
SMITH, PW .
MATHEMATICS OF COMPUTATION, 1984, 43 (168) :491-499