Iterative refinement for linear systems and LAPACK
- 1 October 1997
- journal article
- research article
- Published by Oxford University Press (OUP) in IMA Journal of Numerical Analysis
- Vol. 17 (4) , 495-509
- https://doi.org/10.1093/imanum/17.4.495
Abstract
The technique of iterative refinement for improving the computed solution to a linear system was used on desk calculators and computers in the 1940s and has remained popular. In the 1990s iterative refinement is well supported in software libraries, notably in LAPACK. Although the behaviour of iterative refinement in floating point arithmetic is reasonably well understood, the existing theory is not sufficient to justify the use of fixed precision iterative refinement in all the LAPACK routines in which it is implemented. We present analysis that provides the theoretical support needed for LAPACK. The analysis covers both mixed and fixed precision iterative refinement with an arbitrary number of iterations, makes only a general assumption on the underlying solver, and is relatively short. We identify some remaining open problems.This publication has 3 references indexed in Scilit:
- Stability of the Diagonal Pivoting Method with Partial PivotingSIAM Journal on Matrix Analysis and Applications, 1997
- Iterative Refinement in Floating PointJournal of the ACM, 1967
- Compatibility of approximate solution of linear equations with given error bounds for coefficients and right-hand sidesNumerische Mathematik, 1964