Palindromic 3-stage splitting integrators, a roadmap

Abstract

The implementation of multi-stage splitting integrators is essentially the same as the implementation of the familiar Strang/Verlet method. Therefore multi-stage formulas may be easily incorporated into software that now uses the Strang/Verlet integrator. We study in detail the two-parameter family of palindromic, three-stage splitting formulas and identify choices of parameters that may outperform the Strang/Verlet method. One of these choices leads to a method of effective order four suitable to integrate in time some partial differential equations. Other choices may be seen as perturbations of the Strang method that increase efficiency in molecular dynamics simulations and in Hybrid Monte Carlo sampling.

20 pages, 8 figures, 2 tables

Cited by in corpus (9)

• Geometric integrators and the Hamiltonian Monte Carlo method
• Multi-stage splitting integrators for sampling with modified Hamiltonian Monte Carlo methods
• Geometric Methods for Sampling, Optimisation, Inference and Adaptive Agents
• Adaptive multi-stage integration schemes for Hamiltonian Monte Carlo
• HMC: avoiding rejections by not using leapfrog and some results on the acceptance rate
• The shifted ODE method for underdamped Langevin MCMC
• Cayley Splitting for Second-Order Langevin Stochastic Partial Differential Equations
• Symmetrically processed splitting integrators for enhanced Hamiltonian Monte Carlo sampling
• Is the NUTS algorithm correct?