Algebraic multigrid preconditioning of the Hessian in optimization constrained by a partial differential equation

We construct an algebraic multigrid (AMG) based preconditioner for the reduced Hessian of a linear-quadratic optimization problem constrained by an elliptic partial differential equation. While the preconditioner generalizes a geometric multigrid preconditioner introduced in earlier works, its construction relies entirely on a standard AMG infrastructure built for solving the forward elliptic equation, thus allowing for it to be implemented using a variety of AMG methods and standard packages. Our analysis establishes a clear connection between the quality of the preconditioner and the AMG method used. The proposed strategy has a broad and robust applicability to problems with unstructured grids, complex geometry, and varying coefficients. The method is implemented using the Hypre package and several numerical examples are presented.

[1]  L. R. Scott,et al.  The Mathematical Theory of Finite Element Methods , 1994 .

[2]  G. Biros,et al.  A multilevel algorithm for inverse problems with elliptic PDE constraints , 2008 .

[3]  Robert D. Falgout,et al.  Scaling Hypre's Multigrid Solvers to 100, 000 Cores , 2011, High-Performance Scientific Computing.

[4]  Scott P. MacLachlan,et al.  Theoretical bounds for algebraic multigrid performance: review and analysis , 2014, Numer. Linear Algebra Appl..

[5]  Todd F. Dupont,et al.  Optimal order multilevel preconditioners for regularized ill-posed problems , 2008, Math. Comput..

[6]  Ludmil T. Zikatanov,et al.  Algebraic multigrid methods * , 2016, Acta Numerica.

[7]  Panayot S. Vassilevski,et al.  An improved convergence analysis of smoothed aggregation algebraic multigrid , 2012, Numer. Linear Algebra Appl..

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

[9]  Marian Brezina,et al.  Algebraic multigrid by smoothed aggregation for second and fourth order elliptic problems , 2005, Computing.

[10]  Xiaozhe Hu,et al.  Modifying AMG coarse spaces with weak approximation property to exhibit approximation in energy norm , 2019, SIAM J. Matrix Anal. Appl..

[11]  M. Pachter,et al.  Optimal control of partial differential equations , 1980 .

[12]  Jyoti Saraswat,et al.  Optimal-Order Preconditioners for Linear Systems Arising in the Semismooth Newton Solution of a Class of Control-Constrained Problems , 2014, SIAM J. Matrix Anal. Appl..

[13]  V. E. Henson,et al.  BoomerAMG: a parallel algebraic multigrid solver and preconditioner , 2002 .

[14]  Hehu Xie,et al.  A Multilevel Correction Method for Optimal Controls of Elliptic Equations , 2014, SIAM J. Sci. Comput..

[15]  Cosmin G. Petra,et al.  Multigrid Preconditioning of Linear Systems for Interior Point Methods Applied to a Class of Box-constrained Optimal Control Problems , 2010, SIAM J. Numer. Anal..

[16]  Johann Sienz,et al.  GE Jet Engine Bracket Challenge: A Case Study in Sustainable Design , 2014 .