Fast CG-Based Methods for Tikhonov--Phillips Regularization

  • Authors:
  • Andreas Frommer;Peter Maass

  • Affiliations:
  • -;-

  • Venue:
  • SIAM Journal on Scientific Computing
  • Year:
  • 1999

Quantified Score

Hi-index 0.01

Visualization

Abstract

Tikhonov--Phillips regularization is one of the best-known regularization methods for inverse problems. A posteriori criteria for determining the regularization parameter $\alpha$ require solving $$(*) (A^*A+\alpha I) x =A^* y^{\delta}$$ for different values of $\alpha$.We investigate two methods for accelerating the standard cg-algorithm for solving the family of systems (*). The first one utilizes a stopping criterion for the cg-iterations which depends on $\alpha$ and $\delta$. The second method exploits the shifted structure of the linear systems (*), which allows us to solve (*) simultaneously for different values of $\alpha$. We present numerical experiments for three test problems which illustrate the practical efficiency of the new methods. The experiments as well as theoretical considerations show that run times are accelerated by a factor of at least 3.