Introducing an Explicit Symplectic Integration Scheme for Riemannian
  Manifold Hamiltonian Monte Carlo
  We introduce a recent symplectic integration scheme derived for solving
physically motivated systems with non-separable Hamiltonians. We show its
relevance to Riemannian manifold Hamiltonian Monte Carlo (RMHMC) and provide an
alternative to the currently used generalised leapfrog symplectic integrator,
which relies on solving multiple fixed point iterations to convergence. Via
this approach, we are able to reduce the number of higher-order derivative
calculations per leapfrog step. We explore the implications of this integrator
and demonstrate its efficacy in reducing the computational burden of RMHMC. Our
code is provided in a new open-source Python package, hamiltorch.
