High-performance applications of the nonuniform fast Fourier transform to option pricing


















































High-performance applications of the nonuniform fast Fourier transform to option pricing – Journal of Computational Finance



Skip to main content



Risk.net


  • We develop a fast, flexible framework for Fourier-based option pricing that overcomes key limitations of standard FFT methods.
  • The method combines very high accuracy with high throughput, achieving pricing errors around 10-13 to 10-16 while computing tens of thousands to millions of option prices per second in demanding test cases.
  • The framework is robust across challenging models and payoffs, and its general design makes it suitable for a wide range of pricing applications.

Many sophisticated option pricing models involve random variables whose transition density functions are only tractable in Fourier space. For such models, fast Fourier transform algorithms are key for accelerating the repeated Fourier-space integrations that are common in applications. While simple algorithms have been devised in special cases, these methods tend to be both delicate and strongly model- and payoff-specific. Aiming for much broader applicability, this paper demonstrates how the type 3 nonuniform fast Fourier transform can be combined with modern quadrature methods to construct generic, robust and highly performant computing schemes for Fourier-based option pricing. Two approaches are considered: one (broadly applicable) based on density or cumulative distribution recovery, and one (more specialized) based on the Plancherel theorem. Both are handled by a single algorithm designed to fully use the node placement freedom of the nonuniform fast Fourier transform (NUFFT) and to efficiently accommodate adaptive integration strategies. Numerical tests in the paper include the pricing of European options in a variance gamma Lévy-jump model with a singular density, and in a constant elasticity of variance model extended with stochastic volatility. For both cases, we can achieve high-accuracy results (with errors of around 10-13 to 10-16, say) while pricing tens of thousands – and sometimes even millions – of options per second on a single core.

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

Similar Posts

Leave a Reply

Your email address will not be published. Required fields are marked *