Implicitly restarting the LSQR algorithm

Authors

James Baglama, University of Rhode Island
Daniel J. Richmond, University of Rhode Island

Document Type

Article

Date of Original Version

1-1-2014

Abstract

The LSQR algorithm is a popular method for solving least-squares problems. For some matrices, LSQR may require a prohibitively large number of iterations to determine an approximate solution within a desired accuracy. This paper develops a strategy that couples the LSQR algorithm with an implicitly restarted Golub-Kahan bidiagonalization method to improve the convergence rate. The restart is carried out by first applying the largest harmonic Ritz values as shifts and then using LSQR to compute the solution to the least-squares problem. Theoretical results show how this method is connected to the augmented LSQR method of Baglama, Reichel, and Richmond [Numer. Algorithms, 64 (2013), pp. 263-293] in which the Krylov subspaces are augmented with the harmonic Ritz vectors corresponding to the smallest harmonic Ritz values. Computed examples show the proposed method to be competitive with other methods.

Publication Title, e.g., Journal

Electronic Transactions on Numerical Analysis

Volume

42

Citation/Publisher Attribution

Baglama, James, and Daniel J. Richmond. "Implicitly restarting the LSQR algorithm." Electronic Transactions on Numerical Analysis 42, (2014): 85-105. https://digitalcommons.uri.edu/math_facpubs/71