Computational Optimization and Applications
Algorithm 873: LSTRS: MATLAB software for large-scale trust-region subproblems and regularization
ACM Transactions on Mathematical Software (TOMS)
Computers & Mathematics with Applications
An iterative Lagrange method for the regularization of discrete ill-posed inverse problems
Computers & Mathematics with Applications
A Subspace Minimization Method for the Trust-Region Step
SIAM Journal on Optimization
An Efficient Iterative Approach for Large-Scale Separable Nonlinear Inverse Problems
SIAM Journal on Scientific Computing
Accelerating the LSTRS Algorithm
SIAM Journal on Scientific Computing
A Regularized Gauss-Newton Trust Region Approach to Imaging in Diffuse Optical Tomography
SIAM Journal on Scientific Computing
Hi-index | 0.00 |
We consider large-scale least squares problems where the coefficient matrix comes from the discretization of an operator in an ill-posed problem, and the right-hand side contains noise. Special techniques known as regularization methods are needed to treat these problems in order to control the effect of the noise on the solution. We pose the regularization problem as a quadratically constrained least squares problem. This formulation is equivalent to Tikhonov regularization, and we note that it is also a special case of the trust-region subproblem from optimization. We analyze the trust-region subproblem in the regularization case and we consider the nontrivial extensions of a recently developed method for general large-scale subproblems that will allow us to handle this case. The method relies on matrix-vector products only, has low and fixed storage requirements, and can handle the singularities arising in ill-posed problems. We present numerical results on test problems, on an inverse interpolation problem with field data, and on a model seismic inversion problem with field data.