Regularization by inexact Krylov methods with applications to blind deblurring

05/16/2021
by   Silvia Gazzola, et al.
0

This paper is concerned with the regularization of large-scale discrete inverse problems by means of inexact Krylov methods. Specifically, we derive two new inexact Krylov methods that can be efficiently applied to unregularized or Tikhonov-regularized least squares problems, and we study their theoretical properties, including links with their exact counterparts and strategies to monitor the amount of inexactness. We then apply the new methods to separable nonlinear inverse problems arising in blind deblurring. In this setting inexactness stems from the uncertainty in the parameters defining the blur, which may be recovered using a variable projection method leading to an inner-outer iteration scheme (i.e., one cycle of inner iterations is performed to solve one linear deblurring subproblem for any intermediate values of the blurring parameters computed by a nonlinear least squares solver). The new inexact solvers can naturally handle varying inexact blurring parameters while solving the linear deblurring subproblems, allowing for a much reduced number of total iterations and substantial computational savings with respect to their exact counterparts.

READ FULL TEXT

page 15

page 19

page 21

page 23

research
05/29/2021

An ℓ_p Variable Projection Method for Large-Scale Separable Nonlinear Inverse Problems

The variable projection (VarPro) method is an efficient method to solve ...
research
10/23/2019

Krylov Methods for Low-Rank Regularization

This paper introduces new solvers for the computation of low-rank approx...
research
07/12/2019

Adaptive Regularization Parameter Choice Rules for Large-Scale Problems

This paper derives a new class of adaptive regularization parameter choi...
research
08/11/2020

On inner iterations of the joint bidiagonalization based algorithms for solving large scale linear discrete ill-posed problems

The joint bidiagonalization process of a large matrix pair A,L can be us...
research
01/10/2023

First-projection-then-regularization hybrid algorithms for large-scale general-form regularization

The paper presents first-projection-then-regularization hybrid algorithm...
research
06/26/2023

Subspace Recycling for Sequences of Shifted Systems with Applications in Image Recovery

For many applications involving a sequence of linear systems with slowly...
research
07/08/2020

A stochastic approach to mixed linear and nonlinear inverse problems with applications to seismology

We derive an efficient stochastic algorithm for computational inverse pr...

Please sign up or login with your details

Forgot password? Click here to reset