We propose a hybrid method for solving large sparse linear least squares problems. The method is iterative in nature, as it is based on preconditioned LSQR. However, the preconditioner comes from an orthogonal factorization of a submatrix of the original matrix associated with the least squares problem; the construction of the preconditioner is based on well-known techniques in sparse direct methods. In this talk, we will discuss the quality of the preconditioner and the choices of the submatrix.
Additional Information
Esmond G. Ng (Computational Research Division, Lawrence Berkeley National Laboratory)