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

08/24/2022
by   Håkon Hoel, et al.
0

We construct a higher-order adaptive method for strong approximations of exit times of Itô stochastic differential equations (SDE). 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 turn out to complement each other nicely: adaptive time-stepping 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 numerical integrator and two step-sizes for adaptive time-stepping: 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 numerical integrator and three step-sizes for adaptive time-stepping. For any ξ>0, we prove that the strong error is bounded by 𝒪(h^1-ξ) and 𝒪(h^3/2-ξ) for the first and second method, respectively. Under some conditions, we show that the expected computational cost of both methods are bounded by 𝒪(h^-1 |log(h)|), indicating that both methods are tractable. The theoretical results are supported by numerical examples, and we discuss the potential for extensions that improve the strong convergence rate even further.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
11/05/2020

Magnetic Field Simulations Using Explicit Time Integration With Higher Order Schemes

A transient magneto-quasistatic vector potential formulation involving n...
research
02/08/2021

Higher Strong Order Methods for Itô SDEs on Matrix Lie Groups

In this paper we present a general procedure for designing higher strong...
research
08/08/2023

Boundary-preserving Lamperti-splitting scheme for some Stochastic Differential Equations

We propose and analyse an explicit boundary-preserving scheme for the st...
research
03/30/2022

Semi-explicit integration of second order for weakly coupled poroelasticity

We introduce a semi-explicit time-stepping scheme of second order for li...
research
10/11/2022

Higher Order Far-Field Boundary Conditions for Crystalline Defects

Lattice defects in crystalline materials create long-range elastic field...
research
12/31/2019

Higher-order algorithms for nonlinearly parameterized adaptive control

A set of new adaptive control algorithms is presented. The algorithms ar...
research
10/02/2014

Scalable Nonlinear Learning with Adaptive Polynomial Expansions

Can we effectively learn a nonlinear representation in time comparable t...

Please sign up or login with your details

Forgot password? Click here to reset