Block-ADI preconditioners for solving sparse nonsymmetric linear systems of equations

There is currently a regain of interest in the Alternating Direction Implicit (ADI) algorithms as preconditioners for iterative methods for solving large sparse linear systems, because of their suitability for parallel computing. However, the classical ADI iteration is not directly applicable to Finite Element (FE) matrices, tends to converge too slowly for 3-D problems, and the selection of adequate acceleration parameters, remains a diicult task. In this paper we propose a Block-ADI approach, which overcomes some of these problems. In particular we derive a simple inexpensive heuristic for selecting the acceleration parameters. The new approach can be viewed as a combination of the classical ADI method and a domain decomposition approach.

[1]  William J. Layton,et al.  Peaceman - Rachford procedure and domain decomposition for finite element problems , 1995, Numer. Linear Algebra Appl..

[2]  Yousef Saad,et al.  ILUT: A dual threshold incomplete LU factorization , 1994, Numer. Linear Algebra Appl..

[3]  Youcef Saad,et al.  Highly Parallel Preconditioners for General Sparse Matrices , 1994 .

[4]  Eugene L. Wachspress,et al.  Alternating direction implicit iteration for systems with complex spectra , 1991 .

[5]  Y. Saad,et al.  Krylov Subspace Methods on Supercomputers , 1989 .

[6]  S. Eisenstat,et al.  An experimental study of methods for parallel preconditioned Krylov methods , 1989, C3P.

[7]  S. Lennart Johnsson,et al.  Solving narrow banded systems on ensemble architectures , 1985, TOMS.

[8]  Owe Axelsson,et al.  A survey of preconditioned iterative methods for linear systems of algebraic equations , 1985 .

[9]  H. Elman.  Iterative methods for large, sparse, nonsymmetric systems of linear equations , 1982 .

[10]  O. Axelsson.  Conjugate gradient type methods for unsymmetric and inconsistent systems of linear equations , 1980 .

[11]  J. Meijerink,et al.  An iterative solution method for linear systems of which the coefficient matrix is a symmetric -matrix , 1977 .

[12]  Louis A. Hageman,et al.  Iterative Solution of Large Linear Systems. , 1971 .

[13]  J. J. Douglas.  Alternating direction methods for three space variables , 1962 .

[14]  David Young,et al.  Alternating Direction Implicit Methods , 1962, Adv. Comput..

[15]  J. Gillis,et al.  Matrix Iterative Analysis , 1961 .