Multilevel Monte Carlo for exponential Lévy models

We apply the multilevel Monte Carlo method for option pricing problems using exponential Lévy models with a uniform timestep discretisation. For lookback and barrier options, we derive estimates of the convergence rate of the error introduced by the discrete monitoring of the running supremum of a broad class of Lévy processes. We then use these to obtain upper bounds on the multilevel Monte Carlo variance convergence rate for the Variance Gamma, NIG and a-stable processes. We also provide analysis of a trapezoidal approximation for Asian options. Our method is illustrated by numerical experiments.