Antithetic multilevel Monte Carlo estimation for multidimensional SDEs without Lévy area simulation. Arxiv preprint arXiv:1202.6283

Abstract In this paper we develop antithetic multilevel Monte Carlo (MLMC) esti-mators for multidimensional SDEs driven by Brownian motion. Giles has previously shown that if we combine a numerical approximation with strong order of conver-gence O( ∆ t) with MLMC we can reduce the computational complexity to estimate expected values of Lipschitz functionals of SDE solutions with a root-mean-square error of ε from O(ε−3) to O(ε−2). However, in general, to obtain a rate of strong convergence higher than O( ∆ t1/2) requires simulation, or approximation, of Lévy areas. Recently, Giles and Szpruch [5] constructed an antithetic multilevel estimator that avoids the simulation of Lévy areas and still achieves an MLMC correction vari-ance which is O( ∆ t2) for smooth payoffs and almost O( ∆ t3/2) for piecewise smooth payoffs, even though there is only O( ∆ t1/2) strong convergence. This results in an O(ε−2) complexity for estimating the value of financial European and Asian put and call options. In this paper, we extend these results to more complex payoffs based on the path minimum. To achieve this, an approximation of the Lévy areas is needed, resulting in O( ∆ t3/4) strong convergence. By modifying the antithetic MLMC estimator we are able to obtain O(ε−2 log(ε)2) complexity for estimating financial barrier and lookback options.