This is the accepted version of the paper.This version of the publication may differ from the final published version. We present a joint Monte Carlo-Fourier transform sampling scheme for pricing derivative products under a CGMY model exhibiting jumps of infinite activity and finite or infinite variation. The approach relies on numerical transform inversion with computable error estimates, which allow generating the unknown cumulative distribution function (CDF) of the CGMY process increments at the desired accuracy level. We use this to generate samples and simulate the entire trajectory of the process without need of truncating the process small jumps. We illustrate the computational efficiency of the proposed method by comparing it to the existing methods in the literature on pricing a wide range of option contracts, including path-dependent univariate and multivariate products.
Permanent repository link:The authors would like to thank Gianluca Fusai for interesting comments to a previous version of this paper, Russell Gerrard for his valuable contribution that helped improve the paper, Gerald Rickayzen for useful discussions, and Michele Bianchi, Reiichiro Kawai and Hiroki Masuda for useful suggestions on the implementation of the SR and AR sampling schemes. Usual caveat applies.JEL Classification: G12, G13, C63 †Corresponding author.