Estimating the Backward Error in LSQR

We propose practical stopping criteria for the iterative solution of sparse linear least squares (LS) problems. Although we focus our discussion on the algorithm LSQR of Paige and Saunders, the ideas discussed here may also be applicable to other algorithms. We review why the 2-norm of the projection of the residual vector onto the range of A is a useful measure of convergence, and show how this projection can be estimated efficiently at every iteration of LSQR. We also give practical and cheaply-computable estimates of the backward error for the LS problem. Key words. linear least squares, iterative methods, sparse matrix problems, stopping criteria, backward perturbation analysis, backward error AMS subject classifications. 65F10, 65F20, 65F50, 65G50

• Type: Journal
• Year: 2010
• Where: SIAMMAX
• Authors: Pavel Jiránek, David Titley-Péloquin
• Keywords: Algorithm Lsqr, Backward Error, Propose Practical Stopping, SIAMMAX 2010