Towards Matrix‐Free AD‐Based Preconditioning of KKT Systems in PDE‐Constrained Optimization

Griesse, Roland; Walther, Andrea

doi:10.1002/pamm.200510013

Cited by 2 publications

(2 citation statements)

References 7 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…Discrete second order adjoints (obtained by linearizing the numerical approximations of the model) have been obtained by automatic differentiation [7,9,10].…”

Section: Introductionmentioning

confidence: 99%

Second-order adjoints for solving PDE-constrained optimization problems

Cioaca

Alexe

Sandu

2012

Optimization Methods and Software

View full text Add to dashboard Cite

Inverse problems are of utmost importance in many fields of science and engineering. In the variational approach inverse problems are formulated as PDE-constrained optimization problems, where the optimal estimate of the uncertain parameters is the minimizer of a certain cost functional subject to the constraints posed by the model equations. The numerical solution of such optimization problems requires the computation of derivatives of the model output with respect to model parameters. The first order derivatives of a cost functional (defined on the model output) with respect to a large number of model parameters can be calculated efficiently through first order adjoint sensitivity analysis. Second order adjoint models give second derivative information in the form of matrix-vector products between the Hessian of the cost functional and user defined vectors. Traditionally, the construction of second order derivatives for large scale models has been considered too costly. Consequently, data assimilation applications employ optimization algorithms that use only first order derivative information, like nonlinear conjugate gradients and quasi-Newton methods.In this paper we discuss the mathematical foundations of second order adjoint sensitivity analysis and show that it provides an efficient approach to obtain Hessian-vector products. We study the benefits of using of second order information in the numerical optimization process for data assimilation applications. The numerical studies are performed in a twin experiment setting with a two-dimensional shallow water model. Different scenarios are considered with different discretization approaches, observation sets, and noise levels. Optimization algorithms that employ second order derivatives are tested against widely used methods that require only first order derivatives. Conclusions are drawn regarding the potential benefits and the limitations of using high-order information in large scale data assimilation problems.

show abstract

“…Discrete second order adjoints (obtained by linearizing the numerical approximations of the model) have been obtained by automatic differentiation [7,9,10].…”

Section: Introductionmentioning

confidence: 99%

Second-order adjoints for solving PDE-constrained optimization problems

Cioaca

Alexe

Sandu

2012

Optimization Methods and Software

View full text Add to dashboard Cite

show abstract

“…Discrete second order adjoints (obtained by linearizing the numerical approximations of the model) have been obtained by automatic differentiation [19,40,7].…”

mentioning

confidence: 99%

Discrete second order adjoints in atmospheric chemical transport modeling

Sandu

Zhang

2008

Journal of Computational Physics

View full text Add to dashboard Cite

Abstract. Atmospheric chemical transport models (CTMs) are essential tools for the study of air pollution, for environmental policy decisions, for the interpretation of observational data, and for producing air quality forecasts. Many air quality studies require sensitivity analyses, i.e., the computation of derivatives of the model output with respect to model parameters. The derivatives of a cost functional (defined on the model output) with respect to a large number of model parameters can be calculated efficiently through adjoint sensitivity analysis. While the traditional (first order) adjoint models give the gradient of the cost functional with respect to parameters, second order adjoint models give second derivative information in the form of products between the Hessian of the cost functional and a user defined vector.In this paper we discuss the mathematical foundations of the discrete second order adjoint sensitivity method and present a complete set of computational tools for performing second order sensitivity studies in three-dimensional atmospheric CTMs. The tools include discrete second order adjoints of Runge Kutta and of Rosenbrock time stepping methods for stiff equations together with efficient implementation strategies. Numerical examples illustrate the use of these computational tools in important applications like sensitivity analysis, optimization, uncertainty quantification, and the calculation of directions of maximal error growth in three-dimensional atmospheric CTMs.

show abstract

Towards Matrix‐Free AD‐Based Preconditioning of KKT Systems in PDE‐Constrained Optimization

Cited by 2 publications

References 7 publications

Second-order adjoints for solving PDE-constrained optimization problems

Second-order adjoints for solving PDE-constrained optimization problems

Discrete second order adjoints in atmospheric chemical transport modeling

Contact Info

Product

Resources

About