2008.111: Updating the QR factorization and the least squares problem
2008.111: Sven Hammarling and Craig Lucas (2008) Updating the QR factorization and the least squares problem.
Full text available as:
| PDF - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader 316 Kb |
Abstract
In this paper we treat the problem of updating the QR factorization, with applications to the least squares problem. Algorithms are presented that compute the factorization A1 = Q1 R1, where A1 is the matrix A = QR after it has had a number of rows or columns added or deleted. This is achieved by updating the factors Q and R, and we show this can be much faster than computing the factorization of A1 from scratch. We consider algorithms that exploit the Level 3 BLAS where possible and place no restriction on the dimensions of A or the number of rows and columns added or deleted. For some of our algorithms we present Fortran 77 LAPACK-style code and show the backward error of our updated factors is comparable to the error bounds of the QR factorization of A1.
| Item Type: | MIMS Preprint |
|---|---|
| Subjects: | MSC 2000 > 15 Linear and multilinear algebra; matrix theory MSC 2000 > 65 Numerical analysis |
| MIMS number: | 2008.111 |
| Deposited By: | Sven Hammarling |
| Deposited On: | 14 November 2008 |
Download Statistics: last 4 weeks
Repository Staff Only: edit this item