Minimum residual methods for augmented systems

B. Fischer*, A. Ramage, D. J. Silvester, A. J. Wathen

*Corresponding author for this work
145 Citations (Scopus)

Abstract

For large systems of linear equations, iterative methods provide attractive solution techniques. We describe the applicability and convergence of iterative methods of Krylov subspace type for an important class of symmetric and indefinite matrix problems, namely augmented (or KKT) systems. Specifically, we consider preconditioned minimum residual methods and discuss indefinite versus positive definite preconditioning. For a natural choice of starting vector we prove that when the definite and indefinite preconditioners are related in the obvious way, MINRES (which is applicable in the case of positive definite preconditioning) and full GMRES (which is applicable in the case of indefinite preconditioning) give residual vectors with identical Euclidean norm at each iteration. Moreover, we show that the convergence of both methods is related to a system of normal equations for which the LSQR algorithm can be employed. As a side result, we give a rare example of a non-trivial normal(1) matrix where the corresponding inner product is explicitly known: a conjugate gradient method therefore exists and can be employed in this case.

Original languageEnglish
JournalBIT Numerical Mathematics
Volume38
Issue number3
Pages (from-to)527-543
Number of pages17
ISSN0006-3835
DOIs
Publication statusPublished - 01.01.1998

Fingerprint

Dive into the research topics of 'Minimum residual methods for augmented systems'. Together they form a unique fingerprint.

Cite this