Automatic code generator for higher order integrators

Abstract

Some explicit algorithms for higher order symplectic integration of a large class of Hamilton's equations have recently been discussed by Mushtaq \emph{et. al}. Here we present a Python program for automatic numerical implementation of these algorithms for a given Hamiltonian, both for double precision and multiprecision computations. We provide examples of how to use this program, and illustrate behaviour of both the code generator and the generated solver module(s).

23 pages, 12 figures

References in corpus (2)

• The NumPy array: a structure for efficient numerical computation
• Higher order splitting methods with modified integrators for a class of Hamiltonian systems

Cited by in corpus (3)

• How to Differentiate Collective Variables in Free Energy Codes: Computer-Algebra Code Generation and Automatic Differentiation
• On the Benefits of Surrogate Lagrangians in Optimal Control and Planning Algorithms
• Surrogate Lagrangians for Variational Integrators: High Order Convergence with Low Order Schemes