Speeding up the Euler scheme for killed diffusions
Abstract
Let be a linear diffusion taking values in and consider the standard Euler scheme to compute an approximation to for a given function and a deterministic , where . It is well-known since \cite{GobetKilled} that the presence of killing introduces a loss of accuracy and reduces the weak convergence rate to with being the number of discretisatons. We introduce a drift-implicit Euler method to bring the convergence rate back to , i.e. the optimal rate in the absence of killing, using the theory of recurrent transformations developed in \cite{rectr}. Although the current setup assumes a one-dimensional setting, multidimensional extension is within reach as soon as a systematic treatment of recurrent transformations is available in higher dimensions.
Cite
@article{arxiv.2107.03534,
title = {Speeding up the Euler scheme for killed diffusions},
author = {Umut Çetin and Julien Hok},
journal= {arXiv preprint arXiv:2107.03534},
year = {2021}
}
Comments
Some typos and errors in the earlier version are corrected