Fourier-Cosine Expansion

A second order finite volume IMEX Runge-Kutta scheme for two dimensional PDEs in finance ArXiv ID: 2410.02925 “View on arXiv” Authors: Unknown Abstract In this article we present a novel and general methodology for building second order finite volume implicit-explicit (IMEX) numerical schemes for solving two dimensional financial parabolic PDEs with mixed derivatives. In particular, applications to basket and Heston models are presented. The obtained numerical schemes have excellent properties and are able to overcome the well-documented difficulties related with numerical approximations in the financial literature. The methods achieve true second order convergence with non-regular initial conditions. Besides, the IMEX time integrator allows to overcome the tiny time-step induced by the diffusive term in the explicit schemes, also providing very accurate and non-oscillatory approximations of the Greeks. Finally, in order to assess all the aforementioned good properties of the developed numerical schemes, we compute extremely accurate semi-analytic solutions using multi-dimensional Fourier cosine expansions. A novel technique to truncate the Fourier series for basket options is presented and it is efficiently implemented using multi-GPUs. ...

From characteristic functions to multivariate distribution functions and European option prices by the damped COS method ArXiv ID: 2307.12843 “View on arXiv” Authors: Unknown Abstract We provide a unified framework to obtain numerically certain quantities, such as the distribution function, absolute moments and prices of financial options, from the characteristic function of some (unknown) probability density function using the Fourier-cosine expansion (COS) method. The classical COS method is numerically very efficient in one-dimension, but it cannot deal very well with certain integrands in general dimensions. Therefore, we introduce the damped COS method, which can handle a large class of integrands very efficiently. We prove the convergence of the (damped) COS method and study its order of convergence. The method converges exponentially if the characteristic function decays exponentially. To apply the (damped) COS method, one has to specify two parameters: a truncation range for the multivariate density and the number of terms to approximate the truncated density by a cosine series. We provide an explicit formula for the truncation range and an implicit formula for the number of terms. Numerical experiments up to five dimensions confirm the theoretical results. ...