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.

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

Explicit Stabilized Multirate Method for Stiff Stochastic Differential Equations

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

Balanced implicit Patankar–Euler methods for positive solutions of stochastic differential equations of biological regulatory systems

S-Rock Methods for Stiff Ito Sdes

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

A Derivative-Free Explicit Method with Order 1.0 for Solving Stochastic Delay Differential Equations

A fully adaptive explicit stabilized integrator for advection-diffusion-reaction problems

Steady State and Sign Preserving Semi-Implicit Runge-Kutta Methods for ODEs with Stiff Damping Term.

Strong stability preserving multiderivative time marching methods for stiff reaction–diffusion systems

Stability of the split-step backward Euler scheme for stochastic delay integro-differential equations with Markovian switching

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

Effectiveness of Implicit Methods for Stiff Stochastic Differential Equations

Exponential Adams Bashforth integrators for stiff ODEs, application to cardiac electrophysiology

High order implicit time integration schemes on multiresolution adaptive grids for stiff PDEs

A COMPUTATIONAL METHOD FOR HIGHLY STIFF NONLINEAR SPATIOTEMPORAL SYSTEMS

An Improved Milstein Method for Stiff Stochastic Differential Equations

Stable semi-implicit SDC methods for conservation laws

Tamed Runge-Kutta Methods for SDEs with Super-Linearly Growing Drift and Diffusion Coefficients

A note on stability of the split-step backward Euler method for linear stochastic delay integro-differential equations