LSMR: An iterative algorithm for sparse least-squares problems

Details LSMR: An iterative algorithm for sparse least-squares problems LSMR: An iterative algorithm for sparse least-squares problems

2010-06-04

arxiv.org/abs/1006.0758v2

Computer Science

Mathematical Software

21 pages

Scientific paper

10.1016/j.ijheatmasstransfer.201

An iterative method LSMR is presented for solving linear systems $Ax=b$ and least-squares problem $\min \norm{Ax-b}_2$, with $A$ being sparse or a fast linear operator. LSMR is based on the Golub-Kahan bidiagonalization process. It is analytically equivalent to the MINRES method applied to the normal equation $A\T Ax = A\T b$, so that the quantities $\norm{A\T r_k}$ are monotonically decreasing (where $r_k = b - Ax_k$ is the residual for the current iterate $x_k$). In practice we observe that $\norm{r_k}$ also decreases monotonically. Compared to LSQR, for which only $\norm{r_k}$ is monotonic, it is safer to terminate LSMR early. Improvements for the new iterative method in the presence of extra available memory are also explored.

Affiliated with

Also associated with

No associations

Rating

LSMR: An iterative algorithm for sparse least-squares problems does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.

If you have personal experience with LSMR: An iterative algorithm for sparse least-squares problems, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and LSMR: An iterative algorithm for sparse least-squares problems will most certainly appreciate the feedback.

Rate now

Comments { 0 }

Profile ID: LFWR-SCP-O-722407