A Framework for Deflated and Augmented Krylov Subspace Methods

We consider deflation and augmentation techniques for accelerating the convergence of Krylov subspace methods for the solution of nonsingular linear algebraic systems. Despite some formal similarity, the two techniques are conceptually different from preconditioning. Deflation (in the sense the term...

Full description

Saved in:
Bibliographic Details
Published in:SIAM journal on matrix analysis and applications 2013-01, Vol.34 (2), p.495-518
Main Authors: Gaul, Andre, Gutknecht, Martin H, Liesen, Jorg, Nabben, Reinhard
Format: Article
Language:English
Subjects:
Applied mathematics
Approximation
Augmentation
Breakdown
Eigenvalues
Eigenvectors
Mathematical models
Methods
Operators
Preconditioning
Subspace methods
Subspaces
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:We consider deflation and augmentation techniques for accelerating the convergence of Krylov subspace methods for the solution of nonsingular linear algebraic systems. Despite some formal similarity, the two techniques are conceptually different from preconditioning. Deflation (in the sense the term is used here) "removes" certain parts from the operator making it singular, while augmentation adds a subspace to the Krylov subspace (often the one that is generated by the singular operator); in contrast, preconditioning changes the spectrum of the operator without making it singular. Deflation and augmentation have been used in a variety of methods and settings. Typically, deflation is combined with augmentation to compensate for the singularity of the operator, but both techniques can be applied separately. We introduce a framework of Krylov subspace methods that satisfy a Galerkin condition. It includes the families of orthogonal residual and minimal residual methods. We show that in this framework augmentation can be achieved either explicitly or, equivalently, implicitly by projecting the residuals appropriately and correcting the approximate solutions in a final step. We study conditions for a breakdown of the deflated methods, and we show several possibilities to avoid such breakdowns for the deflated minimum residual (MinRes) method. Numerical experiments illustrate properties of different variants of deflated MinRes analyzed in this paper. [PUBLICATION ABSTRACT]
ISSN:0895-4798
1095-7162
DOI:10.1137/110820713