On the Construction of Deflation-Based Preconditioners

In this article we introduce new bounds on the effective condition number of deflated and preconditioned-deflated symmetric positive definite linear systems. For the case of a subdomain deflation such as that of Nicolaides [SIAM J. Numer. Anal., 24 (1987), pp. 355--365], these theorems can provide direction in choosing a proper decomposition into subdomains. If grid refinement is performed, keeping the subdomain grid resolution fixed, the condition number is insensitive to the grid size. Subdomain deflation is very easy to implement and has been parallelized on a distributed memory system with only a small amount of additional communication. Numerical experiments for a steady-state convection-diffusion problem are included.

[1]  J. Meijerink,et al.  An Efficient Preconditioned CG Method for the Solution of a Class of Layered Problems with Extreme Contrasts in the Coefficients , 1999 .

[2]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[3]  Alan B. Williams,et al.  A Deflation Technique for Linear Systems of Equations , 1998, SIAM J. Sci. Comput..

[4]  Robert J. Plemmons,et al.  Nonnegative Matrices in the Mathematical Sciences , 1979, Classics in Applied Mathematics.

[5]  Cornelis Vuik,et al.  A parallel block-preconditioned GCR method for incompressible flow problems , 2001, Future Gener. Comput. Syst..

[6]  E. F. Kaasschieter,et al.  Preconditioned conjugate gradients for solving singular systems , 1988 .

[7]  Ronald B. Morgan,et al.  A Restarted GMRES Method Augmented with Eigenvectors , 1995, SIAM J. Matrix Anal. Appl..

[8]  Cornelis Vuik,et al.  A Parallel Block Preconditioner Accelerated by Coarse Grid Correction , 2000, HPCN Europe.

[9]  Lois Mansfield.  On the conjugate gradient solution of the Schur complement system obtained from domain decomposition , 1990 .

[10]  S. A. Kharchenko,et al.  Eigenvalue translation based preconditioners for the GMRES(k) method , 1995, Numer. Linear Algebra Appl..

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

[12]  Per Weinerfelt,et al.  Coarse grid correction scheme for implicit multiblock Euler calculations , 1995 .

[13]  Y. Saad.  Analysis of Augmented Krylov Subspace Methods , 1997, SIAM J. Matrix Anal. Appl..

[14]  I. Gustafsson.  A class of first order factorization methods , 1978 .

[15]  S. Eisenstat,et al.  Variational Iterative Methods for Nonsymmetric Systems of Linear Equations , 1983 .

[16]  W. Hackbusch.  Iterative Solution of Large Sparse Systems of Equations , 1993 .

[17]  Cornelis Vuik,et al.  Parallel implementation of a multiblock method with approximate subdomain solution , 1999 .

[18]  Lois Mansfield.  Damped Jacobi Preconditioning and Coarse Grid Deflation for Conjugate Gradient Iteration on Parallel Computers , 1991, SIAM J. Sci. Comput..

[19]  Cornelis Vuik,et al.  GMRESR: a family of nested GMRES methods , 1994, Numer. Linear Algebra Appl..

[20]  Gene H. Golub,et al.  Adaptively Preconditioned GMRES Algorithms , 1998, SIAM J. Sci. Comput..

[21]  Gene H. Golub,et al.  Matrix computations , 1983 .

[22]  Cornelis Vuik,et al.  How fast the Laplace equation was solved in 1995 , 1997 .

[23]  K. Burrage,et al.  Restarted GMRES preconditioned by deflation , 1996 .

[24]  Gene H. Golub,et al.  Matrix computations (3rd ed.) , 1996 .

[25]  H. V. D. Vorst,et al.  The rate of convergence of Conjugate Gradients , 1986 .

[26]  O. Axelsson,et al.  On the eigenvalue distribution of a class of preconditioning methods , 1986 .

[27]  Y. Saad,et al.  Deflated and Augmented Krylov Subspace Techniques , 1997 .

[28]  J. H. Wilkinson.  The algebraic eigenvalue problem , 1966 .

[29]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[30]  R. Nicolaides.  Deflation of conjugate gradients with applications to boundary value problems , 1987 .

[31]  Yousef Saad,et al.  Deflated and Augmented Krylov Subspace Techniques , 1997, Numer. Linear Algebra Appl..