Abstract:We study the binomial, trinomial, and Black-Scholes-Merton models of option pricing. We present fast parallel discrete-time finite-difference algorithms for American call option pricing under the binomial and trinomial models and American put option pricing under the Black-Scholes-Merton model. For $T$-step finite differences, each algorithm runs in $O(\left(T\log^2{T}\right)/p + T)$ time under a greedy scheduler on $p$ processing cores, which is a significant improvement over the $\Theta({T^2}/{p}) + \Omega(T\log{T})$ time taken by the corresponding state-of-the-art parallel algorithm. Even when run on a single core, the $O(T\log^2{T})$ time taken by our algorithms is asymptotically much smaller than the $\Theta(T^2)$ running time of the fastest known serial algorithms. Implementations of our algorithms significantly outperform the fastest implementations of existing algorithms in practice, e.g., when run for $T \approx 1000$ steps on a 48-core machine, our algorithm for the binomial model runs at least $15\times$ faster than the fastest existing parallel program for the same model with the speed-up factor gradually reaching beyond $500\times$ for $T \approx 0.5 \times 10^6$. It saves more than 80\% energy when $T \approx 4000$, and more than 99\% energy for $T > 60,000$. Our option pricing algorithms can be viewed as solving a class of nonlinear 1D stencil (i.e., finite-difference) computation problems efficiently using the Fast Fourier Transform (FFT). To our knowledge, ours are the first algorithms to handle such stencils in $o(T^2)$ time. These contributions are of independent interest as stencil computations have a wide range of applications beyond quantitative finance.

A fast Fourier transform method for Mellin-type option pricing

European Option Pricing With a Fast Fourier Transform Algorithm for Big Data Analysis

Multi-Asset Option Pricing with Exponential Lévy Processes and the Mellin Transform

Pricing American Options by a Fourier Transform Multinomial Tree in a Conic Market

Distributional Mellin calculus in $\mathbb{C}^n$, with applications to option pricing

A Fitted L-Multi-Point Flux Approximation Method for Pricing Options

Generalized Finite Integration Method with Laplace transform for European option pricing under Black–Scholes and Heston models

Efficient Pricing of European-Style Asian Options under Exponential Lévy Processes Based on Fourier Cosine Expansions

Fast and Simple Method for Pricing Exotic Options using Gauss-Hermite Quadrature on a Cubic Spline Interpolation

Option Valuation using Fourier Space Time Stepping

A new radial basis functions method for pricing American options under Merton's jump-diffusion model

An Adaptive and Explicit Fourth Order Runge-Kutta-Fehlberg Method Coupled with Compact Finite Differencing for Pricing American Put Options

Failure of Fourier pricing techniques to approximate the Greeks

Convergence of the Two Point Flux Approximation method and the fitted Two Point Flux Approximation method for options pricing with local volatility function

Fast American Option Pricing using Nonlinear Stencils

Efficient European and American option pricing under a jump-diffusion process

American-type basket option pricing: a simple two-dimensional partial differential equation

Pricing of early-exercise Asian options under Lévy processes based on Fourier cosine expansions

Fast Fourier Transform Option Pricing with Stochastic Interest Rate, Stochastic Volatility and Double Jumps

European Option Pricing Under Generalized Tempered Stable Process: Empirical Analysis

A Fast Fourier Transform Technique for Pricing European Options with Stochastic Volatility and Jump Risk