Annals of Applied Probability

Multilevel Monte Carlo algorithms for Lévy-driven SDEs with Gaussian correction

Steffen Dereich

Full-text: Open access


We introduce and analyze multilevel Monte Carlo algorithms for the computation of $\mathbb {E}f(Y)$, where Y=(Yt)t∈[0, 1] is the solution of a multidimensional Lévy-driven stochastic differential equation and f is a real-valued function on the path space. The algorithm relies on approximations obtained by simulating large jumps of the Lévy process individually and applying a Gaussian approximation for the small jump part. Upper bounds are provided for the worst case error over the class of all measurable real functions f that are Lipschitz continuous with respect to the supremum norm. These upper bounds are easily tractable once one knows the behavior of the Lévy measure around zero.

In particular, one can derive upper bounds from the Blumenthal–Getoor index of the Lévy process. In the case where the Blumenthal–Getoor index is larger than one, this approach is superior to algorithms that do not apply a Gaussian approximation. If the Lévy process does not incorporate a Wiener process or if the Blumenthal–Getoor index β is larger than 4∕3, then the upper bound is of order τ−(4−β)∕(6β) when the runtime τ tends to infinity. Whereas in the case, where β is in [1, 4∕3] and the Lévy process has a Gaussian component, we obtain bounds of order τβ∕(6β−4). In particular, the error is at most of order τ−1∕6.

Article information

Ann. Appl. Probab., Volume 21, Number 1 (2011), 283-311.

First available in Project Euclid: 17 December 2010

Permanent link to this document

Digital Object Identifier

Mathematical Reviews number (MathSciNet)

Zentralblatt MATH identifier

Primary: 60H35: Computational methods for stochastic equations [See also 65C30]
Secondary: 60H05: Stochastic integrals 60H10: Stochastic ordinary differential equations [See also 34F05] 60J75: Jump processes

Multilevel Monte Carlo Komlós–Major–Tusnády coupling weak approximation numerical integration Lévy-driven stochastic differential equation


Dereich, Steffen. Multilevel Monte Carlo algorithms for Lévy-driven SDEs with Gaussian correction. Ann. Appl. Probab. 21 (2011), no. 1, 283--311. doi:10.1214/10-AAP695.

Export citation


  • [1] Applebaum, D. (2004). Lévy Processes and Stochastic Calculus. Cambridge Studies in Advanced Mathematics 93. Cambridge Univ. Press, Cambridge.
  • [2] 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.
  • [3] Bertoin, J. (1998). Lévy Processes. Cambridge Univ. Press.
  • [4] Cohen, S. and Rosiński, J. (2007). Gaussian approximation of multivariate Lévy processes with applications to simulation of tempered stable processes. Bernoulli 13 195–210.
  • [5] Creutzig, J., Dereich, S., Müller-Gronbach, T. and Ritter, K. (2009). Infinite-dimensional quadrature and approximation of distributions. Found. Comput. Math. 9 391–429.
  • [6] Dereich, S. and Heidenreich, F. (2009). A multilevel Monte Carlo algorithm for Lévy driven stochastic differential equations. Preprint.
  • [7] Giles, M. B. (2008). Improved multilevel Monte Carlo convergence using the Milstein scheme. In Monte Carlo and Quasi-Monte Carlo Methods 2006 343–358. Springer, Berlin.
  • [8] Giles, M. B. (2008). Multilevel Monte Carlo path simulation. Oper. Res. 56 607–617.
  • [9] Heinrich, S. (1998). Monte Carlo complexity of global solution of integral equations. J. Complexity 14 151–175.
  • [10] Jacod, J. (2004). The Euler scheme for Lévy driven stochastic differential equations: Limit theorems. Ann. Probab. 32 1830–1872.
  • [11] 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. H. Poincaré Probab. Statist. 41 523–558.
  • [12] Kloeden, P. E. and Platen, E. (1992). Numerical Solution of Stochastic Differential Equations. Applications of Mathematics (New York) 23. Springer, Berlin.
  • [13] Komlós, J., Major, P. and Tusnády, G. (1975). An approximation of partial sums of independent RV’s and the sample DF. I. Z. Wahrsch. Verw. Gebiete 32 111–131.
  • [14] Komlós, J., Major, P. and Tusnády, G. (1976). An approximation of partial sums of independent RV’s, and the sample DF. II. Z. Wahrsch. Verw. Gebiete 34 33–58.
  • [15] Novak, E. (1995). The real number model in numerical analysis. J. Complexity 11 57–73.
  • [16] Protter, P. (2005). Stochastic Integration and Differential Equations, 2nd ed. Stochastic Modelling and Applied Probability 21. Springer, Berlin.
  • [17] Protter, P. and Talay, D. (1997). The Euler scheme for Lévy driven stochastic differential equations. Ann. Probab. 25 393–423.
  • [18] 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.
  • [19] Rubenthaler, S. and Wiktorsson, M. (2003). Improved convergence rate for the simulation of stochastic differential equations driven by subordinated Lévy processes. Stochastic Process. Appl. 108 1–26.
  • [20] Sato, K. (1999). Lévy Processes and Infinitely Divisible Distributions. Cambridge Studies in Advanced Mathematics 68. Cambridge Univ. Press, Cambridge.
  • [21] Talay, D. and Tubaro, L. (1990). Expansion of the global error for numerical schemes solving stochastic differential equations. Stochastic Anal. Appl. 8 483–509.
  • [22] Zaitsev, A. Y. (1998). Multidimensional version of the results of Komlós, Major and Tusnády for vectors with finite exponential moments. ESAIM Probab. Statist. 2 41–108.