Abstract:Abstract We construct a higher-order adaptive method for strong approximations of exit times of Itô stochastic differential equations (SDEs). The method employs a strong Itô–Taylor scheme for simulating SDE paths, and adaptively decreases the step size in the numerical integration as the solution approaches the boundary of the domain. These techniques complement each other nicely: adaptive timestepping improves the accuracy of the exit time by reducing the magnitude of the overshoot of the numerical solution when it exits the domain, and higher-order schemes improve the approximation of the state of the diffusion process. We present two versions of the higher-order adaptive method. The first one uses the Milstein scheme as the numerical integrator and two step sizes for adaptive timestepping: $h$ when far away from the boundary and $h^2$ when close to the boundary. The second method is an extension of the first one using the strong Itô–Taylor scheme of order 1.5 as the numerical integrator and three step sizes for adaptive timestepping. Under some regularity assumptions, we show that for any $\xi>0$, the strong error is ${\mathcal{O}}(h^{1-\xi })$ and ${\mathcal{O}}(h^{3/2-\xi })$ for the first and second method, respectively. Provided quite restrictive commutativity conditions hold for the diffusion coefficient, we further show that the expected computational cost for both methods is ${\mathcal{O}}(h^{-1} \log (h^{-1}))$. This results in a near doubling/trebling of the strong error rate compared to the standard Euler–Maruyama-based approach, while the computational cost rate is kept close to order one. Numerical examples that support the theoretical results are provided, and we discuss the potential for extensions that would further improve the strong convergence rate of the method.

Adaptive Time-Step Semi-Implicit One-Step Taylor Scheme for Stiff Ordinary Differential Equations

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

Semi-Implicit-type order-adaptive CAT2 schemes for systems of balance laws with relaxed source term

On approximate implicit Taylor methods for ordinary differential equations

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

Study of General Taylor-like Explicit Methods in Solving Stiff Ordinary Differential Equations.

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

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

Implicit and semi-implicit second-order time stepping methods for the Richards equation

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

High-Order Schemes of Exponential Time Differencing for Stiff Systems with Nondiagonal Linear Part

A Strong Stability Preserving Analysis for Multistage Two-Derivative Time-Stepping Schemes Based on Taylor Series Conditions

A multistep algorithm for ODEs

Adaptive second derivative multistep methods for solving stiff chemical problems

High-order implicit time integration scheme based on Padé expansions

A Strong Stability Preserving Analysis for Explicit Multistage Two-Derivative Time-Stepping Schemes Based on Taylor Series Conditions

On a new class of BDF and IMEX schemes for parabolic type equations

High-order Accurate Implicit-Explicit Time-Stepping Schemes for Wave Equations on Overset Grids

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

Higher-order adaptive methods for exit times of Itô diffusions

Effectiveness of Implicit Methods for Stiff Stochastic Differential Equations