Abstract:In this paper, we study numerical approximations for stochastic differential equations (SDEs) that use adaptive step sizes. In particular, we consider a general setting where decisions to reduce step sizes are allowed to depend on the future trajectory of the underlying Brownian motion. Since these adaptive step sizes may not be previsible, the standard mean squared error analysis cannot be directly applied to show that the numerical method converges to the solution of the SDE. Building upon the pioneering work of Gaines and Lyons, we instead use rough path theory to establish pathwise convergence for a wide class of adaptive numerical methods on general Stratonovich SDEs (with sufficiently smooth vector fields). To our knowledge, this is the first convergence guarantee that applies to standard solvers, such as the Milstein and Heun methods, with non-previsible step sizes. In our analysis, we require the sequence of adaptive step sizes to be nested and the SDE solver to have unbiased "Lévy area" terms in its Taylor expansion. We conjecture that for adaptive SDE solvers more generally, convergence is still possible provided the method does not introduce "Lévy area bias". We present a simple example where the step size control can skip over previously considered times, resulting in the numerical method converging to an incorrect limit (i.e. not the Stratonovich SDE). Finally, we conclude with an experiment demonstrating the accuracy of Heun's method and a newly introduced Splitting Path-based Runge-Kutta scheme (SPaRK) when used with adaptive step sizes.

Optimal sampling design for global approximation of jump diffusion SDEs

Adaptive step-size control for global approximation of SDEs driven by countably dimensional Wiener process

Randomized Milstein algorithm for approximation of solutions of jump-diffusion SDEs

A higher order approximation method for jump-diffusion SDEs with discontinuous drift coefficient

Lower error bounds and optimality of approximation for jump-diffusion SDEs with discontinuous drift

On the $\varepsilon$-Euler-Maruyama scheme for time inhomogeneous jump-driven SDEs

Strong Approximation of Monotone Stochastic Partial Differential Equations Driven by Multiplicative Noise

An adaptive strong order 1 method for SDEs with discontinuous drift coefficient

Strong convergence of a class of adaptive numerical methods for SDEs with jumps

On the Rate of Convergence of Weak Euler Approximation for Nondegenerate Itô Diffusion and Jump Processes

New algorithms for sampling and diffusion models

A multilevel Monte Carlo algorithm for SDEs driven by countably dimensional Wiener process and Poisson random measure

Global Optimization Via Schrödinger-Föllmer Diffusion.

Optimizing the diffusion coefficient of overdamped Langevin dynamics

Convergence of a Metropolized Integrator for Stochastic Differential Equations with Variable Diffusion Coefficient

Efficient estimation for ergodic diffusion processes sampled at high frequency

Adaptive Weak Approximation of Diffusions with Jumps

On the convergence of adaptive approximations for stochastic differential equations

Importance Sampling for the Empirical Measure of Weakly Interacting Diffusions

A randomized and fully discrete Galerkin finite element method for semilinear stochastic evolution equations

An Adaptive Importance Sampling for Locally Stable Point Processes