Postprocessed integrators for the high order integration of ergodic SDEs

Abstract

The concept of effective order is a popular methodology in the deterministic literature for the construction of efficient and accurate integrators for differential equations over long times. The idea is to enhance the accuracy of a numerical method by using an appropriate change of variables called the processor. We show that this technique can be extended to the stochastic context for the construction of new high order integrators for the sampling of the invariant measure of ergodic systems. The approach is illustrated with modifications of the stochastic $\theta$-method applied to Brownian dynamics, where postprocessors achieving order two are introduced. Numerical experiments, including stiff ergodic systems, illustrate the efficiency and versatility of the approach.

21 pages, to appear in SIAM J. Sci. Comput

Cited by in corpus (7)

• High-order integrator for sampling the invariant distribution of a class of parabolic SPDEs with additive space-time noise
• Optimal explicit stabilized integrator of weak order one for stiff and ergodic stochastic differential equations
• Order conditions for sampling the invariant measure of ergodic stochastic differential equations on manifolds
• Time correlation functions of equilibrium and nonequilibrium Langevin dynamics: Derivations and numerics using random numbers
• A uniformly accurate scheme for the numerical integration of penalized Langevin dynamics
• On the Numerical Stationary Distribution of Overdamped Langevin Equation in Harmonic System
• Hopf algebra structures for the backward error analysis of ergodic stochastic differential equations