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.

A Class of Adams Hybrid Methods for Solving Problems of Stiff ODEs

A CLASS OF STIFFLY STABLE K-STEP SECOND DERIVATIVE METHODS OF ORDER 2K+1

A kind of semi-implicit methods for nonlinear equations

An efficient hybrid implicit time integration method with high accuracy for linear and nonlinear dynamics

Numerical integration of stiff problems using a new time-efficient hybrid block solver based on collocation and interpolation techniques

A Hybrid High-Order Method for a Coupled Stokes-Darcy Problem on General Meshes.

Adaptive second derivative multistep methods for solving stiff chemical problems

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

A Class of Two-Step Steffensen Type Methods with Fourth-Order Convergence.

High Order Stable Multi-Domain Hybrid RKDG and WENO-FD Methods

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

Almost sure exponential stability of θ-method for hybrid stochastic differential equations

A Hybrid High-Order Method for a Class of Quasi-Newtonian Stokes Equations on General Meshes.

Discontinuous Galerkin Method with Staggered Hybridization for a Class of Nonlinear Stokes Equations

A Class of Hybrid Methods of SDBDF and Its Application in Simulation

High Order Iterative Methods Without Derivatives for Solving Nonlinear Equations

Stability of multistep Runge-Kutta methods for systems of functional-differential and functional equations

Algebraically Stable SDIRK Methods with Controllable Numerical Dissipation for First/second-Order Time-Dependent Problems

Explicit stabilized multirate method for stiff differential equations

A Path-Conservative ADER Discontinuous Galerkin Method for Non-Conservative Hyperbolic Systems: Applications to Shallow Water Equations

A Family of Conditionally Explicit Methods for Second-Order ODEs and DAEs: Application in Multibody Dynamics