The Discrete Adjoint Method for Exponential Integration

The implementation of the discrete adjoint method for exponential time differencing (ETD) schemes is considered. This is important for parameter estimation problems that are constrained by stiff time-dependent PDEs when the discretized PDE system is solved using an exponential integrator. We also discuss the closely related topic of computing the action of the sensitivity matrix on a vector, which is required when performing a sensitivity analysis. The PDE system is assumed to be semi-linear and can be the result of a linearization of a nonlinear PDE, leading to exponential Rosenbrock-type methods. We discuss the computation of the derivatives of the $\varphi$-functions that are used by ETD schemes and find that the derivatives strongly depend on the way the $\varphi$-functions are evaluated numerically. A general adjoint exponential integration method, required when computing the gradients, is developed and its implementation is illustrated by applying it to the Krogstad scheme. The applicability of the methods developed here to pattern formation problems is demonstrated using the Swift-Hohenberg model.