English

LSRN: A Parallel Iterative Solver for Strongly Over- or Under-Determined Systems

Data Structures and Algorithms 2012-02-21 v2 Mathematical Software Numerical Analysis

Abstract

We describe a parallel iterative least squares solver named \texttt{LSRN} that is based on random normal projection. \texttt{LSRN} computes the min-length solution to minxRnAxb2\min_{x \in \mathbb{R}^n} \|A x - b\|_2, where ARm×nA \in \mathbb{R}^{m \times n} with mnm \gg n or mnm \ll n, and where AA may be rank-deficient. Tikhonov regularization may also be included. Since AA is only involved in matrix-matrix and matrix-vector multiplications, it can be a dense or sparse matrix or a linear operator, and \texttt{LSRN} automatically speeds up when AA is sparse or a fast linear operator. The preconditioning phase consists of a random normal projection, which is embarrassingly parallel, and a singular value decomposition of size γmin(m,n)×min(m,n)\lceil \gamma \min(m,n) \rceil \times \min(m,n), where γ\gamma is moderately larger than 1, e.g., γ=2\gamma = 2. We prove that the preconditioned system is well-conditioned, with a strong concentration result on the extreme singular values, and hence that the number of iterations is fully predictable when we apply LSQR or the Chebyshev semi-iterative method. As we demonstrate, the Chebyshev method is particularly efficient for solving large problems on clusters with high communication cost. Numerical results demonstrate that on a shared-memory machine, \texttt{LSRN} outperforms LAPACK's DGELSD on large dense problems, and MATLAB's backslash (SuiteSparseQR) on sparse problems. Further experiments demonstrate that \texttt{LSRN} scales well on an Amazon Elastic Compute Cloud cluster.

Keywords

Cite

@article{arxiv.1109.5981,
  title  = {LSRN: A Parallel Iterative Solver for Strongly Over- or Under-Determined Systems},
  author = {Xiangrui Meng and Michael A. Saunders and Michael W. Mahoney},
  journal= {arXiv preprint arXiv:1109.5981},
  year   = {2012}
}

Comments

19 pages

R2 v1 2026-06-21T19:11:13.533Z