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.

[1]  Christian Bischof,et al.  Parallel calculation of sensitivity derivatives for aircraft design using automatic differentiation , 1994 .

[2]  Boghos B. Joulakian,et al.  Double ionization of helium by fast electrons: use of correlated two electron wavefunctions , 1992 .

[3]  Mohamed Masmoudi,et al.  Computation of high order derivatives in optimal shape design , 1994 .

[4]  Arthur B. Baggeroer,et al.  An overview of matched field methods in ocean acoustics , 1993 .

[5]  S. Jones,et al.  Single and Double Ionization of Atoms by Photons, Electrons, and Ions , 2003 .

[6]  I. Charpentier.  Adjoint modelling experiments on eruptive columns , 2007 .

[7]  I. Charpentier.  Variational coupling of Plinian column models and data: Application to El Chichón volcano , 2008 .

[8]  Philippe Roux,et al.  MODE AND WAVENUMBER INVERSION IN SHALLOW WATER USING AN ADJOINT METHOD , 2004 .

[9]  Andreas Griewank,et al.  Evaluating higher derivative tensors by forward propagation of univariate Taylor series , 2000, Math. Comput..

[10]  Andreas Griewank,et al.  Evaluating derivatives - principles and techniques of algorithmic differentiation, Second Edition , 2000, Frontiers in applied mathematics.

[11]  Jamal Berakdar,et al.  Comparative theoretical study of (e, 3e) on helium: Coulomb-waves versus close-coupling approach , 2002 .

[12]  J. Briggs,et al.  Triply-differential cross sections for ionisation of hydrogen atoms by electrons and positrons , 1989 .

[13]  Toichiro Kinoshita,et al.  GROUND STATE OF THE HELIUM ATOM , 1957 .

[14]  R. E. Wengert,et al.  A simple automatic derivative evaluation program , 1964, Commun. ACM.

[15]  W. Kuperman,et al.  Computational Ocean Acoustics , 1994 .

[16]  Andreas Griewank,et al.  Algorithm 755: ADOL-C: a package for the automatic differentiation of algorithms written in C/C++ , 1996, TOMS.

[17]  Isabelle Charpentier,et al.  High order cross derivative computation for the differential cross section of double ionization of helium by electron impact , 2005 .