Iterative refinement enhances the stability of QR factorization methods for solving linear equations

Higham, Nicholas J. (1991) Iterative refinement enhances the stability of QR factorization methods for solving linear equations. BIT Numerical Mathematics, 31. pp. 447-468. ISSN 1572-9125

[thumbnail of g75558h0wv05713u.pdf] PDF
g75558h0wv05713u.pdf
Restricted to Repository staff only

Download (1MB)

Abstract

Iterative refinement is a well-known technique for improving the quality of an approximate solution to a linear system. In the traditional usage residuals are computed in extended precision, but more recent work has shown that fixed precision is sufficient to yield benefits for stability. We extend existing results to show that fixed precision iterative refinement renders anarbitrary linear equations solver backward stable in a strong, componentwise sense, under suitable assumptions. Two particular applications involving theQR factorization are discussed in detail: solution of square linear systems and solution of least squares problems. In the former case we show that one step of iterative refinement suffices to produce a small componentwise relative backward error. Our results are weaker for the least squares problem, but again we find that iterative refinement improves a componentwise measure of backward stability. In particular, iterative refinement mitigates the effect of poor row scaling of the coefficient matrix, and so provides an alternative to the use of row interchanges in the HouseholderQR factorization. A further application of the results is described to fast methods for solving Vandermonde-like systems.

Item Type: Article
Uncontrolled Keywords: Iterative refinement - linear system - least squares problem - QR factorization - Gaussian elimination - partial pivoting - rounding error analysis - backward error - componentwise error bounds - Householder transformations - Givens transformations - confluent Vandermonde-like matrices
Subjects: MSC 2010, the AMS's Mathematics Subject Classification > 15 Linear and multilinear algebra; matrix theory
MSC 2010, the AMS's Mathematics Subject Classification > 65 Numerical analysis
Depositing User: Ms Lucy van Russelt
Date Deposited: 28 Jun 2006
Last Modified: 20 Oct 2017 14:12
URI: https://eprints.maths.manchester.ac.uk/id/eprint/334

Actions (login required)

View Item View Item