Fast higher-order derivative tensors with Rapsodia. 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.
Keywords for this software
References in zbMATH (referenced in 7 articles , 1 standard article )
Showing results 1 to 7 of 7.
- Goldsztejn, Alexandre; Cruz, Jorge; Carvalho, Elsa: Convergence analysis and adaptive strategy for the certified quadrature over a set defined by inequalities (2014)
- Griewank, Andreas; Lehmann, Lutz; Leovey, Hernan; Zilberman, Marat: Automatic evaluations of cross-derivatives (2014)
- Charpentier, I.: On higher-order differentiation in nonlinear mechanics (2012)
- Reed, James A.; Utke, Jean; Abdel-Khalik, Hany S.: Combining automatic differentiation methods for high-dimensional nonlinear models (2012)
- Lampoh, Komlanvi; Charpentier, Isabelle; Daya, El Mostafa: A generic approach for the solution of nonlinear residual equations. III: Sensitivity computations (2011)
- Charpentier, I.; Utke, J.: Fast higher-order derivative tensors with Rapsodia (2009)
- Charpentier, Isabelle; dal Cappello, Claude; Utke, Jean: Efficient higher-order derivatives of the hypergeometric function (2008)