| Full text | |
| Author(s): |
Borges, Leonardo S.
;
Viloche Bazan, Fermin S.
;
Cunha, Maria C. C.
Total Authors: 3
|
| Document type: | Journal article |
| Source: | COMPUTATIONAL & APPLIED MATHEMATICS; v. 34, n. 3, p. 23-pg., 2015-10-01. |
| Abstract | |
The numerical treatment of large-scale discrete ill-posed problems is often accomplished iteratively by projecting the original problem onto a -dimensional subspace with acting as regularization parameter. Hence, to filter out the contribution of noise in the computed solution, the iterative process must be stopped early. In this paper, we analyze in detail a stopping rule for LSQR proposed recently by the authors, and show how to extend it to Krylov subspace methods such as GMRES, MINRES, etc. Like the original rule, the extended version works well without requiring a priori knowledge about the error norm and stops automatically (k) over bar + 1 after steps where (k) over bar is the computed regularization parameter. The performance of the stopping rule on several test problems is illustrated numerically. (AU) | |