Multiplicative noise, fast convolution and pricing

In this work we detail the application of a fast convolution algorithm to compute high-dimensional integrals in the context of multiplicative noise stochastic processes. The algorithm provides a numerical solution to the problem of characterizing conditional probability density functions at arbitrary times, and we apply it successfully to quadratic and piecewise linear diffusion processes. The ability to reproduce statistical features of financial return time series, such as thickness of the tails and scaling properties, makes these processes appealing for option pricing. Since exact analytical results are lacking, we exploit the fast convolution as a numerical method alternative to Monte Carlo simulation both in the objective and risk-neutral settings. In numerical sections we document how fast convolution outperforms Monte Carlo both in speed and efficiency terms.