The following technical report is available from
http://aib.informatik.rwth-aachen.de:
Symbolic vs. Algorithmic Differentiation of GSL Integration Routines
Niloofar Safiran and Uwe Naumann
AIB 2015-14
Forward and reverse modes of algorithmic differentiation (AD) transform
implementations of multivariate vector functions as computer programs
into tangent and adjoint code, respectively. The adjoint mode is of
particular interest in large-scale functions due to the independence of
its computational cost on the number of free variables. The additional
memory requirement for the computation of derivatives of the output with
respect to parameters by a fully algorithmic method (derived by AD) can
quickly become prohibitive for large values of n. This can be reduced
significantly by the symbolic approach to differentiation of the
underlying integration routine. Vectorizing gsl routines for integration
and applying symbolic adjoint on them has considerably less memory
requirement with nearly the same runtime overhead and in most cases
faster convergence in comparison with algorithmic adjoint.
