Abstract:SIAM Journal on Numerical Analysis, Volume 62, Issue 1, Page 500-532, February 2024. In this paper, we introduce a new simple approach to developing and establishing the convergence of splitting methods for a large class of stochastic differential equations (SDEs), including additive, diagonal, and scalar noise types. The central idea is to view the splitting method as a replacement of the driving signal of an SDE, namely, Brownian motion and time, with a piecewise linear path that yields a sequence of ODEs—which can be discretized to produce a numerical scheme. This new way of understanding splitting methods is inspired by, but does not use, rough path theory. We show that when the driving piecewise linear path matches certain iterated stochastic integrals of Brownian motion, then a high order splitting method can be obtained. We propose a general proof methodology for establishing the strong convergence of these approximations that is akin to the general framework of Milstein and Tretyakov. That is, once local error estimates are obtained for the splitting method, then a global rate of convergence follows. This approach can then be readily applied in future research on SDE splitting methods. By incorporating recently developed approximations for iterated integrals of Brownian motion into these piecewise linear paths, we propose several high order splitting methods for SDEs satisfying a certain commutativity condition. In our experiments, which include the Cox–Ingersoll–Ross model and additive noise SDEs (noisy anharmonic oscillator, stochastic FitzHugh–Nagumo model, underdamped Langevin dynamics), the new splitting methods exhibit convergence rates of [math] and outperform schemes previously proposed in the literature.

ADAPTIVE TIME SPLITTING METHOD FOR MULTI-SCALE EVOLUTIONARY PARTIAL DIFFERENTIAL EQUATIONS

Adaptive Partitioned Methods for the Time-Accurate Approximation of the Evolutionary Stokes–Darcy System

Time multiscale modeling of sorption kinetics I: uniformly accurate schemes for highly oscillatory advection-diffusion equation

An Adaptive Linear Time Stepping Algorithm For Second-Order Linear Evolution Problems

Parallel in time partially explicit splitting scheme for high contrast multiscale problems

Splitting methods with complex coefficients for linear and nonlinear evolution equations

Practical splitting methods for the adaptive integration of nonlinear evolution equations. Part I: Construction of optimized schemes and pairs of schemes

Error estimation and adaptivity for differential equations with multiple scales in time

A New Splitting Method for Time-dependent Convection-dominated Diffusion Problems

Adaptive Third Order Adams-Bashforth Time Stepping for Extended Boussinesq Equations

A numerical strategy to discretize and solve the Poisson equation on dynamically adapted multiresolution grids for time-dependent streamer discharge simulations

Stability-Optimized High Order Methods and Stiffness Detection for Pathwise Stiff Stochastic Differential Equations

Adaptive Space-Time Domain Decomposition for Multiphase Flow in Porous Media with Bound Constraints

High Order Splitting Methods for SDEs Satisfying a Commutativity Condition

A second-order adaptive time filter algorithm with different subdomain variable time steps for the evolutionary Stokes/Darcy model

A splitting method for SDEs with locally Lipschitz drift: Illustration on the FitzHugh-Nagumo model

On the convergence of adaptive approximations for stochastic differential equations

A parallel splitting wavelet method for 2D conservation laws

Local time stepping methods and discontinuous Galerkin methods applied to diffusion advection reaction equations

Stochastic methods for solving high-dimensional partial differential equations

Adaptive timestepping strategies for nonlinear stochastic systems