A scalable parallel algorithm for incomplete factor preconditioning

被引:83
作者
Hysom, D [1 ]
Pothen, A
机构
[1] Old Dominion Univ, Norfolk, VA 23529 USA
[2] NASA, ICASE, Langley Res Ctr, Hampton, VA 23681 USA
关键词
incomplete factorization; ILU; preconditioning; parallel preconditioning;
D O I
10.1137/S1064827500376193
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
We describe a parallel algorithm for computing incomplete factor (ILU) preconditioners. The algorithm attains a high degree of parallelism through graph partitioning and a two-level ordering strategy. Both the subdomains and the nodes within each subdomain are ordered to preserve concurrency. We show through an algorithmic analysis and through computational results that this algorithm is scalable. Experimental results include timings on three parallel platforms for problems with up to 20 million unknowns running on up to 216 processors. The resulting preconditioned Krylov solvers have the desirable property that the number of iterations required for convergence is insensitive to the number of processors.
引用
收藏
页码:2194 / 2215
页数:22
相关论文
共 37 条
[1]  
Axelsson O., 1994, ITERATIVE SOLUTION M
[2]  
BALAY S, 1999, PETSC HOME PAGE
[3]   PARALLELIZATION OF ROBUST MULTIGRID METHODS - ILU FACTORIZATION AND FREQUENCY DECOMPOSITION METHOD [J].
BASTIAN, P ;
HORTON, G .
SIAM JOURNAL ON SCIENTIFIC AND STATISTICAL COMPUTING, 1991, 12 (06) :1457-1470
[4]  
Benzi M., 1999, ELECTRON T NUMER ANA, V8, P88
[5]   Object-oriented design of preconditioned iterative methods in diffpack [J].
Bruaset, AM ;
Langtangen, HP .
ACM TRANSACTIONS ON MATHEMATICAL SOFTWARE, 1997, 23 (01) :50-80
[6]  
CHAN TF, 1997, ICASE LARC INTERDISC, V4, P167
[7]   Experimental study of ILU preconditioners for indefinite matrices [J].
Chow, E ;
Saad, Y .
JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS, 1997, 86 (02) :387-414
[8]   An object-oriented framework for block preconditioning [J].
Chow, E ;
Heroux, MA .
ACM TRANSACTIONS ON MATHEMATICAL SOFTWARE, 1998, 24 (02) :159-183
[9]   ORDERING METHODS FOR PRECONDITIONED CONJUGATE-GRADIENT METHODS APPLIED TO UNSTRUCTURED GRID PROBLEMS [J].
DAZEVEDO, EF ;
FORSYTH, PA ;
TANG, WP .
SIAM JOURNAL ON MATRIX ANALYSIS AND APPLICATIONS, 1992, 13 (03) :944-961
[10]   Ordering strategies and related techniques to overcome the trade-off between parallelism and convergence in incomplete factorizations [J].
Doi, S ;
Washio, T .
PARALLEL COMPUTING, 1999, 25 (13-14) :1995-2014