Skip to content
Licensed Unlicensed Requires Authentication Published by De Gruyter November 13, 2013

Implementation and analysis of an adaptive multilevel Monte Carlo algorithm

  • Håkon Hoel EMAIL logo , Erik von Schwerin , Anders Szepessy and Raúl Tempone


We present an adaptive multilevel Monte Carlo (MLMC) method for weak approximations of solutions to Itô stochastic differential equations (SDE). The work [Oper. Res. 56 (2008), 607–617] proposed and analyzed an MLMC method based on a hierarchy of uniform time discretizations and control variates to reduce the computational effort required by a single level Euler–Maruyama Monte Carlo method from 𝒪( TOL -3) to 𝒪( TOL -2log( TOL -1)2) for a mean square error of 𝒪( TOL 2). Later, the work [Lect. Notes Comput. Sci. Eng. 82, Springer-Verlag, Berlin (2012), 217–234] presented an MLMC method using a hierarchy of adaptively refined, non-uniform time discretizations, and, as such, it may be considered a generalization of the uniform time discretization MLMC method. This work improves the adaptive MLMC algorithms presented in [Lect. Notes Comput. Sci. Eng. 82, Springer-Verlag, Berlin (2012), 217–234] and it also provides mathematical analysis of the improved algorithms. In particular, we show that under some assumptions our adaptive MLMC algorithms are asymptotically accurate and essentially have the correct complexity but with improved control of the complexity constant factor in the asymptotic analysis. Numerical tests include one case with singular drift and one with stopped diffusion, where the complexity of a uniform single level method is 𝒪( TOL -4). For both these cases the results confirm the theory, exhibiting savings in the computational cost for achieving the accuracy 𝒪( TOL ) from 𝒪( TOL -3) for the adaptive single level algorithm to essentially 𝒪( TOL -2log( TOL -1)2) for the adaptive MLMC algorithm.

Funding source: Royal Institute of Technology in Stockholm

Award Identifier / Grant number: Dahlquist fellowship

Funding source: Department of Scientific Computing in Florida State University

Funding source: University of Austin Subcontract

Award Identifier / Grant number: 024550

Funding source: VR project

Award Identifier / Grant number: “Effektiva numeriska metoder för stokastiska differentialekvationer med tillämpningar”

Funding source: Center for Industrial and Applied Mathematics (CIAM) at the Royal Institute of Technology

Funding source: King Abdullah University of Science and Technology (KAUST)

The authors would like to thank Mike Giles and the two anonymous reviewers for valuable comments.

Received: 2012-11-12
Accepted: 2013-09-12
Published Online: 2013-11-13
Published in Print: 2014-03-01

© 2014 by Walter de Gruyter Berlin/Boston

Downloaded on 23.2.2024 from
Scroll to top button