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.

Free knot linear interpolation and the Milstein scheme for stochastic differential equations

A strong order $3/4$ method for SDEs with discontinuous drift coefficient

An Analysis of the Milstein Scheme for SPDEs without a Commutative Noise Condition

Sharp lower error bounds for strong approximation of SDEs with piecewise Lipschitz continuous drift coefficient

The Optimal Rate for Linear KB-splines and LKB-splines Approximation of High Dimensional Continuous Functions and its Application

Estimation of Knots in Linear Spline Models

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

An Antithetic Multilevel Monte Carlo-Milstein Scheme for Stochastic Partial Differential Equations

Higher order time discretization method for a class of semilinear stochastic partial differential equations with multiplicative noise

A tamed-adaptive Milstein scheme for stochastic differential equations with low regularity coefficients

Regular polynomial interpolation and approximation of global solutions of linear partial differential equations

Approximation of SDEs: a stochastic sewing approach

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

Asymptotic Error for the Milstein Scheme for SDEs with Stochastic Evaluation Times

Strong convergence of an adaptive time-stepping Milstein method for SDEs with monotone coefficients

KNOT PLACEMENT FOR B-SPLINE CURVE APPROXIMATION VIA L∞,1-Norm AND DIFFERENTIAL EVOLUTION ALGORITHM

On the convergence of adaptive approximations for stochastic differential equations

Milstein schemes and antithetic multilevel Monte Carlo sampling for delay McKean-Vlasov equations and interacting particle systems

Higher order numerical methods for SDEs without globally monotone coefficients

On the complexity of strong approximation of stochastic differential equations with a non-Lipschitz drift coefficient