A Spectral Element Method to Price European Options. I. Single Asset with and without Jump Diffusion

  • Authors:
  • Wuming Zhu;David A. Kopriva

  • Affiliations:
  • Department of Mathematics, The Florida State University, Tallahassee, USA;Department of Mathematics, The Florida State University, Tallahassee, USA

  • Venue:
  • Journal of Scientific Computing
  • Year:
  • 2009

Quantified Score

Hi-index 0.00

Visualization

Abstract

We develop a spectral element method to price single factor European options with and without jump diffusion. The method uses piecewise high order Legendre polynomial expansions to approximate the option price represented pointwise on a Gauss-Lobatto mesh within each element, which allows an exact representation of the non-smooth payoff function. The convolution integral is approximated by high order Gauss-Lobatto quadratures. A second order implicit/explicit (IMEX) approximation is used to integrate in time, with the convolution integral integrated explicitly. The method is spectrally accurate (exponentially convergent) in space for the solution and Greeks, and second-order accurate in time. The spectral element solution to the Black-Scholes equation is ten to one hundred times faster than commonly used second order finite difference methods.