This discussion paper resulted in a publication in 'Quantitative Finance', 2010, 10, 177-194.<P> When using an Euler discretisation to simulate a mean-reverting square root process, one runs into the problem that while the process itself is guaranteed to be nonnegative, the discretisation is not....</p>

At the time of writing this article, Fourier inversion is the computational method of choice for a fast and accurate calculation of plain vanilla option prices in models with an analytically available characteristic function. Shifting the contour of integration along the complex plane allows for...

The characteristic functions of many affine jump-diffusion models, such as Heston’s stochastic volatility model and all of its extensions, involve multivalued functions such as the complex logarithm. If we restrict the logarithm to its principal branch, as is done in most software packages,...

Using an Euler discretization to simulate a mean-reverting CEV process gives rise to the problem that while the process itself is guaranteed to be nonnegative, the discretization is not. Although an exact and efficient simulation algorithm exists for this process, at present this is not the case...

Guo and Hung (2007) recently studied the complex logarithm present in the characteristic function of Heston's stochastic volatility model. They proposed an algorithm for the evaluation of the characteristic function that is claimed to preserve its continuity. We show their algorithm is correct,...

The characteristic functions of many affine jump-diffusion models, such as Heston’s stochastic volatility model and all of its extensions, involve multivalued functions such as the complex logarithm. If we restrict the logarithm to its principal branch, as is done in most software packages,...

When using an Euler discretisation to simulate a mean-reverting square root process, one runs into the problem that while the process itself is guaranteed to be nonnegative, the discretisation is not. Although an exact and efficient simulation algorithm exists for this process, at present this...

At the time of writing this article, Fourier inversion is the computational method of choice for a fast and accurate calculation of plain vanilla option prices in models with an analytically available characteristic function. Shifting the contour of integration along the complex plane allows for...

When using an Euler discretisation to simulate a mean-reverting square root process, one runs into the problem that while the process itself is guaranteed to be nonnegative, the discretisation is not. Although an exact and efficient simulation algorithm exists for this process, at present this...

The characteristic functions of many affine jump-diffusion models, such as Heston’s stochastic volatility model and all of its extensions, involve multivalued functions such as the complex logarithm. If we restrict the logarithm to its principal branch, as is done in most software packages,...