Algebraic multigrid preconditioning of the Hessian in PDE-constrained optimization

01/30/2020
by   Andrew T. Barker, et al.
0

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.

READ FULL TEXT

page 1

page 11

page 14

research
11/13/2019

A new preconditioner for elliptic PDE-constrained optimization problems

We propose a preconditioner to accelerate the convergence of the GMRES i...
research
04/19/2023

Operations for D-Algebraic Functions

A function is differentially algebraic (or simply D-algebraic) if there ...
research
02/13/2023

Technical Note: PDE-constrained Optimization Formulation for Tumor Growth Model Calibration

We discuss solution algorithms for calibrating a tumor growth model usin...
research
08/25/2020

Optimization with learning-informed differential equation constraints and its applications

Inspired by applications in optimal control of semilinear elliptic parti...
research
04/17/2022

Effective numerical computation of p(x)-Laplace equations in 2D

In this article we implement a method for the computation of a nonlinear...
research
09/01/2019

A least-squares Galerkin approach to gradient and Hessian recovery for nondivergence-form elliptic equations

We propose a least-squares method involving the recovery of the gradient...
research
06/20/2021

A PDE approach to centroidal tessellations of domains

We introduce a class of systems of Hamilton-Jacobi equations that charac...

Please sign up or login with your details

Forgot password? Click here to reset