On the eigenvalue distribution of a class of preconditioning methods
Numerische Mathematik
On the rate of convergence of the preconditioned conjugate gradient method
Numerische Mathematik
Mixed and hybrid finite element methods
Mixed and hybrid finite element methods
A preconditioned iterative method for saddlepoint problems
SIAM Journal on Matrix Analysis and Applications
Fast iterative solution of stabilised Stokes systems part II: using general block preconditioners
SIAM Journal on Numerical Analysis
Iterative solution methods
Mathematical physiology
A domain embedding preconditioner for the Lagrange multiplier systems
Mathematics of Computation
Multigrid
The Numerical Solution of the Steady State Solid Fuel Ignition Model and Its Optimal Control
SIAM Journal on Scientific Computing
SIAM Journal on Control and Optimization
Preconditioners for saddle point problems arising in computational fluid dynamics
Applied Numerical Mathematics
Uniform preconditioners for the time dependent Stokes problem
Numerische Mathematik
A Multigrid Scheme for Elliptic Constrained Optimal Control Problems
Computational Optimization and Applications
SIAM Journal on Scientific Computing
SIAM Journal on Matrix Analysis and Applications
Multigrid Algorithms for Inverse Problems with Linear Parabolic PDE Constraints
SIAM Journal on Scientific Computing
Multigrid Methods for PDE Optimization
SIAM Review
Hi-index | 0.00 |
This paper is devoted to the numerical treatment of linear optimality systems (OS) that arise in connection with inverse problems for partial differential equations. If such inverse problems are regularized by Tikhonov regularization, then it follows from standard theory that the associated OS is well-posed, provided that the regularization parameter $\alpha$ is positive and that the involved state equation satisfies suitable assumptions. We explain and analyze how certain mapping properties of the operators appearing in the OS can be employed to define efficient preconditioners for finite element (FE) approximations of such systems. The key feature of the scheme is that the number of iterations needed to solve the preconditioned problem by the minimal residual method is bounded independently of the mesh parameter $h$, used in the FE discretization, and increases only moderately as $\alpha\rightarrow0$. More specifically, if the stopping criterion for the iteration process is defined in terms of the associated energy norm, then the number of iterations required (in the severely ill-posed case) cannot grow faster than $O((\ln(\alpha))^2)$. Our analysis is based on a careful study of the operators involved, which yields the distribution of the eigenvalues of the preconditioned OS. Finally, the theoretical results are illuminated by a number of numerical experiments addressing both a model problem studied by Borzi, Kunisch, and Kwak [SIAM J. Control Optim., 41 (2003), pp. 1477-1497] and an inverse problem arising in connection with electrocardiography [Nielsen, Cai, and Lysaker, Math. Biosci., 210 (2007), pp. 523-553].