Abstract
We present a randomized iterative algorithm that exponentially converges in the mean square to the minimum ℓ2-norm least squares solution of a given linear system of equations. The expected number of arithmetic operations required to obtain an estimate of given accuracy is proportional to the squared condition number of the system multiplied by the number of nonzero entries of the input matrix. The proposed algorithm is an extension of the randomized Kaczmarz method that was analyzed by Strohmer and Vershynin.
Original language | English (US) |
---|---|
Pages (from-to) | 773-793 |
Number of pages | 21 |
Journal | SIAM Journal on Matrix Analysis and Applications |
Volume | 34 |
Issue number | 2 |
DOIs | |
State | Published - 2013 |
Keywords
- Iterative method
- LAPACK
- Linear least squares
- Minimum-length solution
- Overdetermined system
- Random sampling
- Randomized algorithms
- Sparse matrix
- Underdetermined system
ASJC Scopus subject areas
- Analysis