zbMATH — the first resource for mathematics

Fast higher-order derivative tensors with Rapsodia. (English) Zbl 1166.65009
Summary: 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.

65D25 Numerical differentiation
68W30 Symbolic computation and algebraic computation
65Y20 Complexity and performance of numerical algorithms
65Y15 Packaged methods for numerical algorithms
AD01; Rapsodia
Full Text: DOI