Abstract:Stochastic differential equations (SDE) often exhibit large random transitions. This property, which we denote as pathwise stiffness, causes transient bursts of stiffness which limit the allowed step size for common fixed time step explicit and drift-implicit integrators. We present four separate methods to efficiently handle this stiffness. First, we utilize a computational technique to derive stability-optimized adaptive methods of strong order 1.5 for SDEs. The resulting explicit methods are shown to exhibit substantially enlarged stability regions which allows for them to solve pathwise stiff biological models orders of magnitude more efficiently than previous methods like SRIW1 and Euler-Maruyama. Secondly, these integrators include a stiffness estimator which allows for automatically switching between implicit and explicit schemes based on the current stiffness. In addition, adaptive L-stable strong order 1.5 implicit integrators for SDEs and stochastic differential algebraic equations (SDAEs) in mass-matrix form with additive noise are derived and are demonstrated as more efficient than the explicit methods on stiff chemical reaction networks by nearly 8x. Lastly, we developed an adaptive implicit-explicit (IMEX) integration method based off of a common method for diffusion-reaction-convection PDEs and show numerically that it can achieve strong order 1.5. These methods are benchmarked on a range of problems varying from non-stiff to extreme pathwise stiff and demonstrate speedups between 5x-6000x while showing computationally infeasibility of fixed time step integrators on many of these test equations.

Weak Convergence of the Split-Step Backward Euler Method for Stochastic Delay Integro-Differential Equations

On one-step numerical schemes of weak convergence for SDEs with super-linear coefficients

Weak error analysis for strong approximation schemes of SDEs with super-linear coefficients II: finite moments and higher-order schemes

Strong and weak convergence rates of finite element method for stochastic partial differential equation with non-sided Lipschitz coefficient

The tamed Milstein method for commutative stochastic differential equations with non-globally Lipschitz continuous coefficients

An adaptive Euler-Maruyama scheme for SDDEs

Weak Backward Error Analysis for SDEs

Stochastic fractional integro-differential equations with weakly singular kernels: Well-posedness and Euler–Maruyama approximation

Stochastic C-stability and B-consistency of explicit and implicit Euler-type schemes

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

Strong convergence rates for long-time approximations of SDEs with non-globally Lipschitz continuous coefficients

Long-time weak convergence analysis of a semi-discrete scheme for stochastic Maxwell equations

Weak approximation schemes for SDEs with super-linearly growing coefficients

Strong convergence of multiscale truncated Euler-Maruyama method for super-linear slow-fast stochastic differential equations

Strong Predictor-Corrector Euler-Maruyama Methods for Stochastic Differential Equations with Markovian Switching

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

Weak Convergence Of Tamed Exponential Integrators for Stochastic Differential Equations

Semi-implicit Euler–Maruyama method for non-linear time-changed stochastic differential equations

Convergence and stability of the semi-tamed Milstein method for commutative stochastic differential equations with non-globally Lipschitz continuous coefficients

Strong convergence in the infinite horizon of numerical methods for stochastic differential equations

Strong convergence of the truncated Euler–Maruyama method for stochastic functional differential equations