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.

Effectiveness of Implicit Methods for Stiff Stochastic Differential Equations

Effectiveness of implicit methods for stiff stochastic differential equations

A family of fully implicit Milstein methods for stiff stochastic differential equations with multiplicative noise

Mean-square convergence rates of implicit Milstein type methods for SDEs with non-Lipschitz coefficients

Stability of the Drift-Implicit and Double-Implicit Milstein Schemes for Nonlinear SDEs

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

Explicit Stabilized Multirate Method for Stiff Stochastic Differential Equations

Stability of the Semi-Implicit Milstein Numerical Method for the Stochastic Differential Delay Equations

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

S-Rock Methods for Stiff Ito Sdes

Analysis of multiscale methods for stochastic differential equations

A Family of Fully Implicit Strong Itô-Taylor Numerical Methods for Stochastic Differential Equations

Implicit-explicit Relaxation Runge-Kutta Methods: Construction, Analysis and Applications to PDEs

Delay-Dependent Stability Analysis of Numerical Methods for Stochastic Delay Differential Equations

Analysis of Numerical Techniques for Multiscale Stochastic Dynamical Systems

Performance Evaluation of Single-step Explicit Exponential Integration Methods on Stiff Ordinary Differential Equations

Convergence and Stability of the Semi-Implicit Euler Method for Linear Stochastic Delay Integro-Differential Equations.

EFFICIENT STOCHASTIC ASYMPTOTIC-PRESERVING IMPLICIT-EXPLICIT METHODS FOR TRANSPORT EQUATIONS WITH DIFFUSIVE SCALINGS AND RANDOM INPUTS

A fast implicit difference scheme for solving the generalized time-space fractional diffusion equations with variable coefficients

Exponential Stability of Implicit Numerical Solution for Nonlinear Neutral Stochastic Differential Equations with Time-Varying Delay and Poisson Jumps