Higher-order adaptive methods for exit times of It\^o diffusions
Abstract
We construct a higher-order adaptive method for strong approximations of exit times of It\^o stochastic differential equations (SDE). The method employs a strong It\^o--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: when far away from the boundary and when close to the boundary. The second method is an extension of the first one using the strong It\^o--Taylor scheme of order 1.5 as numerical integrator and three step-sizes for adaptive time-stepping. For any , we prove that the strong error is bounded by and for the first and second method, respectively, and the expected computational cost for both methods is . Theoretical results are supported by numerical examples, and we discuss the potential for extensions that improve the strong convergence rate even further.
Cite
@article{arxiv.2208.11288,
title = {Higher-order adaptive methods for exit times of It\^o diffusions},
author = {Håkon Hoel and Sankarasubramanian Ragunathan},
journal= {arXiv preprint arXiv:2208.11288},
year = {2022}
}
Comments
The computational cost results in Theorems 2.8 and 2.11 are improved through using a connection between the SDE and an absorbing-boundary Fokker--Planck equation. Correction of a few typos