Journal of Computational Finance

Risk.net

Optimal damping with a hierarchical adaptive quadrature for efficient Fourier pricing of multi-asset options in Lévy models

Christian Bayer, Chiheb Ben Hammouda, Antonis Papapantoleon, Michael Samet and Raul Tempone

  • We propose a rule for the choice of the damping parameters which results in a smoother integrand and accelerates the convergence of numerical quadrature methods.
  • A design for dimension-adaptive sparse grids quadrature (ASGQ) method to alleviate the curse of dimensionality when pricing multi-asset options using a Fourier-based approach is proposed.
  • The advantage of using the proposed damping rule and employing ASGQ for basket and rainbow options under Geometric Brownian motion, Variance Gamma and Normal Inverse Gaussian is shown.
  • We demonstrate the achieved computational gains when using our proposed approach compared to applying the COS method or using Monte Carlo in the physical space.

Efficiently pricing multi-asset options is a challenging problem in quantitative finance. When the characteristic function is available, Fourier-based methods are more competitive than alternative techniques because the integrand in the frequency space often has a greater regularity than that in the physical space. However, when designing a numerical quadrature method for most Fourier pricing approaches, two key aspects affecting the numerical complexity should be carefully considered: the choice of damping parameters to ensure integrability and control the regularity class of the integrand, and the effective treatment of high dimensionality. To address these challenges we propose an efficient numerical method for pricing European multi-asset options that is based on two complementary ideas. First, we smooth the Fourier integrand via an optimized choice of damping parameters based on a proposed optimization rule. Then, we employ sparsification and dimension-adaptivity techniques to accelerate the convergence of the quadrature in high dimensions. The extensive numerical study on basket and rainbow options under the multivariate geometric Brownian motion and some Lévy models demonstrates the advantages of adaptivity and the damping rule on the numerical complexity of quadrature methods. Moreover, for the tested two-asset examples, the proposed approach outperforms the Fourier cosine expansion (COS) method in terms of computation time. Finally, we show a greater speed-up for up to six dimensions compared with the Monte Carlo method.

Sorry, our subscription options are not loading right now

Please try again later. Get in touch with our customer services team if this issue persists.

New to Risk.net? View our subscription options

You need to sign in to use this feature. If you don’t have a Risk.net account, please register for a trial.

Sign in
You are currently on corporate access.

To use this feature you will need an individual account. If you have one already please sign in.

Sign in.

Alternatively you can request an individual account here