A multilevel algorithm for inverse problems with elliptic PDE constraints

We present a multilevel algorithm for the solution of a source identification problem in which the forward problem is an elliptic partial differential equation on the 2D unit box. The Hessian corresponds to a Tikhonov-regularized first-kind Fredholm equation. Our method uses an approximate Hessian operator for which, first, the spectral decomposition is known, and second, there exists a fast algorithm that can perform the spectral transform. Based on this decomposition we propose a conjugate gradients solver which we precondition with a multilevel subspace projection scheme. The coarse-level preconditioner is an exact solve and the finer-levels preconditioner is one step of the scaled Richardson iteration. As a model problem, we consider the 2D-Neumann Poisson problem with variable coefficients and partial observations. The approximate Hessian for this case is the Hessian related to a problem with constant coefficients and full observations. We can use a fast cosine transform to compute the spectral transforms. We examine the effect of using Galerkin or level-discretized Hessian operators and we provide results from numerical experiments that indicate the effectiveness of the method for full and partial observations.

[1]  J. Craggs.  Applied Mathematical Sciences , 1973 .

[2]  D. Brandt,et al.  Multi-level adaptive solutions to boundary-value problems math comptr , 1977 .

[3]  Yau Shu Wong,et al.  Spectral multigrid methods for elliptic equations II , 1982 .

[4]  Wolfgang Hackbusch,et al.  Multi-grid methods and applications , 1985, Springer series in computational mathematics.

[5]  William L. Briggs,et al.  A multigrid tutorial , 1987 .

[6]  Wilhelm Heinrichs.  Line relaxation for spectral multigrid methods , 1988 .

[7]  Josef Stoer,et al.  Numerische Mathematik 1 , 1989 .

[8]  J. King.  On the construction of preconditioners by subspace decomposition , 1990 .

[9]  A 3D spectral multigrid method , 1991 .

[10]  J. King,et al.  Multilevel algorithms for ill-posed problems , 1992 .

[11]  R. Kress,et al.  Inverse Acoustic and Electromagnetic Scattering Theory , 1992 .

[12]  Eyal Arian,et al.  MULTIGRID ONE SHOT METHODS FOR OPTIMAL CONTROL PROBLEMS: INFINITE DIMENSIONAL CONTROL , 1994 .

[13]  O. Axelsson.  Iterative solution methods , 1995 .

[14]  Peter Monk,et al.  A new algorithm in electromagnetic inverse scattering theory with an application to medical imaging , 1997 .

[15]  Martin Hanke,et al.  Two-level preconditioners for regularized inverse problems I: Theory , 1999, Numerische Mathematik.

[16]  William L. Briggs,et al.  A multigrid tutorial, Second Edition , 2000 .

[17]  V. Schulz,et al.  Multigrid optimization in applications , 2000 .

[18]  S. McCormick,et al.  A multigrid tutorial (2nd ed.) , 2000 .

[19]  Y. Y. Belov,et al.  Inverse Problems for Partial Differential Equations , 2002 .

[20]  Barbara Kaltenbacher.  V-cycle convergence of some multigrid methods for ill-posed problems , 2003, Math. Comput..

[21]  A. Borzì.  MULTIGRID METHODS FOR OPTIMALITY SYSTEMS , 2004 .

[22]  George Biros,et al.  DYNAMIC DATA-DRIVEN INVERSION FOR TERASCALE SIMULATIONS: REAL-TIME IDENTIFICATION OF AIRBORNE CONTAMINANTS , 2005, ACM/IEEE SC 2005 Conference (SC'05).

[23]  Alfio Borzì,et al.  High-order discretization and multigrid solution of elliptic nonlinear constrained optimal control problems , 2007 .

[24]  George Biros,et al.  Multigrid Algorithms for Inverse Problems with Linear Parabolic PDE Constraints , 2008, SIAM J. Sci. Comput..

[25]  Aaas News,et al.  Book Reviews , 1893, Buffalo Medical and Surgical Journal.