Splitting methods

RI McLachlan, GRW Quispel - Acta Numerica, 2002 - cambridge.org
I thought that instead of the great number of precepts of which logic is composed, I would
have enough with the four following ones, provided that I made a firm and unalterable …

Geometric integrators for ODEs

RI McLachlan, GRW Quispel - Journal of Physics A: Mathematical …, 2006 - iopscience.iop.org
Geometric integration is the numerical integration of a differential equation, while preserving
one or more of its' geometric'properties exactly, ie to within round-off error. Many of these …

Symplectic fourth-order maps for the collisional N-body problem

W Dehnen, DM Hernandez - Monthly Notices of the Royal …, 2016 - academic.oup.com
We study analytically and experimentally certain symplectic and time-reversible N-body
integrators which employ a Kepler solver for each pair-wise interaction, including the …

Switching integrators reversibly in the astrophysical N-body problem

DM Hernandez, W Dehnen - Monthly Notices of the Royal …, 2023 - academic.oup.com
We present a simple algorithm to switch between N-body time integrators in a reversible
way. We apply it to planetary systems undergoing arbitrarily close encounters and highly …

Symplectic integration for the collisional gravitational N-body problem

DM Hernandez, E Bertschinger - Monthly Notices of the Royal …, 2015 - academic.oup.com
We present a new symplectic integrator designed for collisional gravitational N-body
problems which makes use of Kepler solvers. The integrator is also reversible and …

Accurate orbit propagation in the presence of planetary close encounters

D Amato, G Bau, C Bombardelli - Monthly Notices of the Royal …, 2017 - academic.oup.com
We present an efficient strategy for the numerical propagation of small Solar system objects
undergoing close encounters with massive bodies. The trajectory is split into several …

Fast and reliable symplectic integration for planetary system N-body problems

DM Hernandez - Monthly Notices of the Royal Astronomical …, 2016 - academic.oup.com
We apply one of the exactly symplectic integrators, which we call HB15, of Hernandez &
Bertschinger, along with the Kepler problem solver of Wisdom & Hernandez, to solve …

Should N-body integrators be symplectic everywhere in phase space?

DM Hernandez - Monthly Notices of the Royal Astronomical …, 2019 - academic.oup.com
Symplectic integrators are the preferred method of solving conservative N-body problems in
cosmological, stellar cluster, and planetary system simulations because of their superior …

[HTML][HTML] Tensor-product adaptive grids based on coordinate transformations

PA Zegeling - Journal of computational and applied mathematics, 2004 - Elsevier
In this paper we discuss a two-dimensional adaptive grid method that is based on a tensor-
product approach. Adaptive grids are a commonly used tool for increasing the accuracy and …

An efficient multiple time-scale reversible integrator for the gravitational N-body problem

B Leimkuhler - Applied numerical mathematics, 2002 - Elsevier
A large gravitational (or classical atomic) N-body simulation typically includes fast binary
stars, planet-moon systems, or other tightly bound objects, demanding a small timestep and …