ias15: a fast, adaptive, high-order integrator for gravitational dynamics, accurate to machine precision over a billion orbits

We present IAS15, a 15th-order integrator to simulate gravitational dynamics. The integrator is based on a Gau-Radau quadrature and can handle conservative as well as non-conservative forces. We develop a step-size control that can automatically choose an optimal timestep. The algorithm can handle close encounters and high-eccentricity orbits. The systematic errors are kept well below machine precision, and long-term orbit integrations over 10 9 orbits show that IAS15 is optimal in the sense that it follows Brouwer's law, i.e. the energy error behaves like a random walk. Our tests show that IAS15 is superior to a mixed-variable symplectic integrator and other popular integrators, including high-order ones, in…

• NS
  National Science Foundation

  Awards: AST-0807444, 0807444