Mean exit times and the multilevel Monte Carlo method

Higham, Desmond and Mao, Xuerong and Roj, Mikolaj and Song, Qingshuo and Yin, George (2013) Mean exit times and the multilevel Monte Carlo method. SIAM/ASA Journal on Uncertainty Quantification (JUQ), 1 (1). pp. 2-18. (

[thumbnail of 88380.pdf] PDF. Filename: 88380.pdf
Final Published Version

Download (304kB)


Numerical methods for stochastic differential equations are relatively inefficient when used to approximate mean exit times. In particular, although the basic Euler–Maruyama method has weak order equal to one for approximating the expected value of the solution, the order reduces to one half when it is used in a straightforward manner to approximate the mean value of a (stopped) exit time. Consequently, the widely used standard approach of combining an Euler–Maruyama discretization with a Monte Carlo simulation leads to a computationally expensive procedure. In this work, we show that the multilevel approach developed by Giles [Oper. Res., 56 (2008), pp. 607–617] can be adapted to the mean exit time context. In order to justify the algorithm, we analyze the strong error of the discretization method in terms of its ability to approximate the exit time. We then show that the resulting multilevel algorithm improves the expected computational complexity by an order of magnitude, in terms of the required accuracy. Numerical results are provided to illustrate the analysis.


Higham, Desmond ORCID logoORCID:, Mao, Xuerong ORCID logoORCID:, Roj, Mikolaj, Song, Qingshuo and Yin, George;