Discrete first- and second-order adjoints and automatic differentiation for the sensitivity analysis of dynamic models

We describe the use of first- and second-order tangent-linear and adjoint models of the residual of linear-implicit autonomous differential algebraic systems in the context of an extrapolated Euler scheme. The derivative code compiler dcc is applied to a C-implementation of the residual to get first derivative code. Second-(and higher-)order derivative models are obtained by reapplication of dcc to its own output. The resulting solver serves as a first proof of concept of a new platform for source-level manipulation of mathematical models that is currently under development at RWTH Aachen University.