Deflated and augmented Krylov subspace methods: Basic Facts and a Breakdown-free deflated MINRES

In this paper we consider deflation and augmentation techniques for accelerating the convergence of Krylov subspace methods for the solution of nonsingular linear algebraic systems. The two techniques are conceptually different from preconditioning. Deflation "removes" certain parts from the operator, while augmentation adds a subspace to the Krylov subspace. Both approaches have been used in a variety of methods and settings. For Krylov subspace methods that satisfy a (Petrov-) Galerkin condition we show that augmentation can in general be achieved implicitly by projecting the residuals appropriately and correcting the approximate solutions in a final step. In this context, we analyze known methods to deflate CG, GMRes and MinRes. Our analysis reveals that the recently proposed RMinRes method can break down. We show how such breakdowns can be avoided by choosing a special initial guess, and we derive a breakdown-free deflated MinRes method. In numerical experiments we study the properties of different variants of MinRes analyzed in this paper.

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

[2]  Reinhard Nabben,et al.  Deflation and Balancing Preconditioners for Krylov Subspace Methods Applied to Nonsymmetric Matrices , 2008, SIAM J. Matrix Anal. Appl..

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

[4]  Cornelis Vuik,et al.  A Comparison of Deflation and Coarse Grid Correction Applied to Porous Media Flow , 2004, SIAM J. Numer. Anal..

[5]  Cornelis Vuik,et al.  A Comparison of Two-Level Preconditioners Based on Multigrid and Deflation , 2010, SIAM J. Matrix Anal. Appl..

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

[7]  Cornelis Vuik,et al.  On the Construction of Deflation-Based Preconditioners , 2001, SIAM J. Sci. Comput..

[8]  Cornelis Vuik,et al.  A comparison of abstract versions of deflation, balancing and additive coarse grid correction preconditioners , 2008, Numer. Linear Algebra Appl..

[9]  Cornelis Vuik,et al.  A Comparison of Deflation and the Balancing Preconditioner , 2005, SIAM J. Sci. Comput..

[10]  Oliver G. Ernst,et al.  Analysis of acceleration strategies for restarted minimal residual methods , 2000 .

[11]  Z. Dostál Conjugate gradient method with preconditioning by projector , 1988 .

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

[13]  E. Sturler,et al.  Large‐scale topology optimization using preconditioned Krylov subspace methods with recycling , 2007 .

[14]  Cornelis Vuik,et al.  Comparison of Two-Level Preconditioners Derived from Deflation, Domain Decomposition and Multigrid Methods , 2009, J. Sci. Comput..

[15]  L. Yu. Kolotilina,et al.  Twofold deflation preconditioning of linear algebraic systems. I. Theory , 1998 .

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

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

[18]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .

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

[20]  H. Walker,et al.  GMRES On (Nearly) Singular Systems , 1997, SIAM J. Matrix Anal. Appl..

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

[22]  Frédéric Guyomarc'h,et al.  A Deflated Version of the Conjugate Gradient Algorithm , 1999, SIAM J. Sci. Comput..