Computer Science Technical Reports
CS at VT

On Consistency Properties of Discrete Adjoint Linear Multistep Methods

Sandu, Adrian (2007) On Consistency Properties of Discrete Adjoint Linear Multistep Methods. Technical Report TR-07-40, Computer Science, Virginia Tech.

Full text available as:
PDF - Requires Adobe Acrobat Reader or other PDF viewer.
dadj_lmm_report.pdf (286249)

Abstract

In this paper we analyze the consistency properties of discrete adjoints of linear multistep methods. Discrete adjoints are very popular in optimization and control since they can be constructed automatically by reverse mode automatic differentiation. The consistency analysis reveals that the discrete linear multistep adjoints are, in general, inconsistent approximations of the adjoint ODE solution along the trajectory. However, the discrete adjoints at the initial time (and therefore the discrete adjoint gradients) converge to the adjoint ODE solution with the same order as the original linear multistep method. Discrete adjoints inherit the zero-stability properties of the forward method. Numerical results confirm the theoretical findings.

Item Type:Departmental Technical Report
Keywords:Linear multistep methods, discrete adjoints
Subjects:Computer Science > Numerical Analysis
ID Code:1007
Deposited By:Constantinescu, Mr. Emil M.
Deposited On:28 November 2007