• Bernoulli
  • Volume 23, Number 2 (2017), 927-950.

Multilevel path simulation for weak approximation schemes with application to Lévy-driven SDEs

Denis Belomestny and Tigran Nagapetyan

Full-text: Access denied (no subscription detected)

We're sorry, but we are unable to provide you with the full text of this article because we are not able to identify you as a subscriber. If you have a personal subscription to this journal, then please login. If you are already logged in, then you may need to update your profile to register your subscription. Read more about accessing full-text


In this paper, we discuss the possibility of using multilevel Monte Carlo (MLMC) approach for weak approximation schemes. It turns out that by means of a simple coupling between consecutive time discretisation levels, one can achieve the same complexity gain as under the presence of a strong convergence. We exemplify this general idea in the case of weak Euler schemes for Lévy-driven stochastic differential equations. The numerical performance of the new “weak” MLMC method is illustrated by several numerical examples.

Article information

Bernoulli, Volume 23, Number 2 (2017), 927-950.

Received: September 2014
Revised: August 2015
First available in Project Euclid: 4 February 2017

Permanent link to this document

Digital Object Identifier

Mathematical Reviews number (MathSciNet)

Zentralblatt MATH identifier

Lévy-driven stochastic differential equations multilevel Monte Carlo weak approximation schemes


Belomestny, Denis; Nagapetyan, Tigran. Multilevel path simulation for weak approximation schemes with application to Lévy-driven SDEs. Bernoulli 23 (2017), no. 2, 927--950. doi:10.3150/15-BEJ764.

Export citation


  • [1] Asmussen, S. and Rosiński, J. (2001). Approximations of small jumps of Lévy processes with a view towards simulation. J. Appl. Probab. 38 482–493.
  • [2] Bally, V. and Talay, D. (1995). The Euler scheme for stochastic differential equations: Error analysis with Malliavin calculus. Math. Comput. Simulation 38 35–41.
  • [3] Dereich, S. (2011). Multilevel Monte Carlo algorithms for Lévy-driven SDEs with Gaussian correction. Ann. Appl. Probab. 21 283–311.
  • [4] Devroye, L. (1986). Nonuniform Random Variate Generation. New York: Springer.
  • [5] Dia, E.H.A. (2013). Error bounds for small jumps of Lévy processes. Adv. in Appl. Probab. 45 86–105.
  • [6] Fournier, N. (2011). Simulation and approximation of Lévy-driven stochastic differential equations. ESAIM Probab. Stat. 15 233–248.
  • [7] Giles, M. and Xia, Y. (2014). Multilevel Monte Carlo for exponential Lévy models. Preprint. Available at arXiv:1403.5309.
  • [8] Giles, M.B. (2008). Multilevel Monte Carlo path simulation. Oper. Res. 56 607–617.
  • [9] Glasserman, P. (2004). Monte Carlo Methods in Financial Engineering. Applications of Mathematics (New York) 53. New York: Springer.
  • [10] Heinrich, S. (1998). Monte Carlo complexity of global solution of integral equations. J. Complexity 14 151–175.
  • [11] Ikeda, N. and Watanabe, S. (1981). Stochastic Differential Equations and Diffusion Processes. North-Holland Mathematical Library 24. North-Holland: Amsterdam.
  • [12] Jacod, J., Kurtz, T.G., Méléard, S. and Protter, P. (2005). The approximate Euler method for Lévy driven stochastic differential equations. Ann. Inst. Henri Poincaré Probab. Stat. 41 523–558.
  • [13] Jourdain, B. and Kohatsu-Higa, A. (2011). A review of recent results on approximation of solutions of stochastic differential equations. In Stochastic Analysis with Financial Applications. Progress in Probability 65 121–144. Basel: Birkhäuser/Springer Basel AG.
  • [14] Kusuoka, S. (2004). Approximation of expectation of diffusion processes based on Lie algebra and Malliavin calculus. In Advances in Mathematical Economics. Vol. 6. Adv. Math. Econ. 6 69–83. Tokyo: Springer.
  • [15] Lyons, T. and Victoir, N. (2004). Cubature on Wiener space. Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 460 169–198.
  • [16] Ninomiya, M. and Ninomiya, S. (2009). A new higher-order weak approximation scheme for stochastic differential equations and the Runge–Kutta method. Finance Stoch. 13 415–443.
  • [17] Ninomiya, S. and Victoir, N. (2008). Weak approximation of stochastic differential equations and application to derivative pricing. Appl. Math. Finance 15 107–121.
  • [18] Omland, S., Hefter, M., Ritter, K., Brugger, C., De Schryver, C., Wehn, N. and Kostiuk, A. (2015). Exploiting mixed-precision arithmetics in a Multilevel Monte Carlo Approach on FPGAs. In FPGA Based Accelerators for Financial Applications 191–220. Berlin: Springer.
  • [19] Platen, E. and Bruti-Liberati, N. (2010). Numerical Solution of Stochastic Differential Equations with Jumps in Finance. Stochastic Modelling and Applied Probability 64. Berlin: Springer.
  • [20] Protter, P. and Talay, D. (1997). The Euler scheme for Lévy driven stochastic differential equations. Ann. Probab. 25 393–423.
  • [21] Rubenthaler, S. (2003). Numerical simulation of the solution of a stochastic differential equation driven by a Lévy process. Stochastic Process. Appl. 103 311–349.
  • [22] Talay, D. and Tubaro, L. (1990). Expansion of the global error for numerical schemes solving stochastic differential equations. Stoch. Anal. Appl. 8 483–509.
  • [23] Tanaka, H. and Kohatsu-Higa, A. (2009). An operator approach for Markov chain weak approximations with an application to infinite activity Lévy driven SDEs. Ann. Appl. Probab. 19 1026–1062.