Adaptive multilevel Monte Carlo simulation

by Håkon Hoel, Erik Von Schwerin, Anders Szepessy, Raul Tempone
Year: 2012

Bibliography

Hoel, Håkon; von Schwerin, Erik; Szepessy, Anders; Tempone, Raúl "Adaptive multilevel Monte Carlo simulation." Numerical analysis of multiscale computations, 217–234, Lect. Notes Comput. Sci. Eng., 82, Springer, Heidelberg, 2012.

Abstract

This work generalizes a multilevel forward Euler Monte Carlo method introduced in Michael B. Giles. (Michael Giles. Oper. Res. 56(3):607–617, 2008.) for the approximation of expected values depending on the solution to an Itô stochastic differential equation. The work (Michael Giles. Oper. Res. 56(3):607– 617, 2008.) proposed and analyzed a forward Euler multilevelMonte Carlo method based on a hierarchy of uniform time discretizations and control variates to reduce the computational effort required by a standard, single level, Forward Euler Monte Carlo method. This work introduces an adaptive hierarchy of non uniform time discretizations, generated by an adaptive algorithmintroduced in (AnnaDzougoutov et al. Raùl Tempone. Adaptive Monte Carlo algorithms for stopped diffusion. In Multiscale methods in science and engineering, volume 44 of Lect. Notes Comput. Sci. Eng., pages 59–88. Springer, Berlin, 2005; Kyoung-Sook Moon et al. Stoch. Anal. Appl. 23(3):511–558, 2005; Kyoung-Sook Moon et al. An adaptive algorithm for ordinary, stochastic and partial differential equations. In Recent advances in adaptive computation, volume 383 of Contemp. Math., pages 325–343. Amer. Math. Soc., Providence, RI, 2005.). This form of the adaptive algorithm generates stochastic, path dependent, time steps and is based on a posteriori error expansions first developed in (Anders Szepessy et al. Comm. Pure Appl. Math. 54(10):1169– 1214, 2001). Our numerical results for a stopped diffusion problem, exhibit savings in the computational cost to achieve an accuracy of ϑ(TOL),from(TOL−3) , from using a single level version of the adaptive algorithm to ϑ(((TOL−1)log(TOL))2).

1439-7358

 

Keywords

Computational finance Monte Carlo multilevel adaptivity Weak approximation error control Euler–Maruyama method A posteriori error estimates backward dual functions Adjoints