Fast higher-order derivative tensors with Rapsodia

Authors:
I. Charpentier;J. Utke
Affiliations:
Laboratoire de Physique et Mecanique des Materiaux, Metz Cedex 1, France;Mathematics and Computer Science Division, Argonne National Laboratory, Argonne, IL, USA
Venue:
Optimization Methods & Software
Year:
2009

Citing 4
Cited 2

Algorithm 755: ADOL-C: a package for the automatic differentiation of algorithms written in C/C++

ACM Transactions on Mathematical Software (TOMS)
Evaluating derivatives: principles and techniques of algorithmic differentiation

Evaluating derivatives: principles and techniques of algorithmic differentiation
Evaluating higher derivative tensors by forward propagation of univariate Taylor series

Mathematics of Computation
A simple automatic derivative evaluation program

Communications of the ACM

Using Multicomplex Variables for Automatic Computation of High-Order Derivatives

ACM Transactions on Mathematical Software (TOMS)
Convergence analysis and adaptive strategy for the certified quadrature over a set defined by inequalities

Journal of Computational and Applied Mathematics

Quantified Score

Hi-index	0.00

Visualization

Abstract

A number of practical problems in physics can be solved by using accurate higher-order derivatives. Such derivatives can be obtained with automatic differentiation. However, one has to be concerned with the complexity of computing higher-order derivative tensors even for a modest order and number of independents. Initial experiments using univariate Taylor polynomials with interpolation and operator overloading with unrolled loops showed better runtimes than using other automatic differentiation tools. Motivated by these results, we developed the Rapsodia code generator that produces Fortran and C++libraries for the most common intrinsics. Here we explain the algorithmic approach, implementation, and present test results on a select set of applications. Further details on the Rapsodia tool, and an example for user extensions are given in the Appendix.