Fast Exponential Time Integration for Pricing Options in Stochastic Volatility Jump Diffusion Models

The stochastic volatility jump diffusion model with jumps in both return and volatility leads to a two-dimensional partial integro-differential equation (PIDE). We exploit a fast exponential time integration scheme to solve this PIDE. After spatial discretization and temporal integration, the solution of the PIDE can be formulated as the action of an exponential of a block Toeplitz matrix on a vector. The shift-invert Arnoldi method is employed to approximate this product. To reduce the computational cost, matrix splitting is combined with the multigrid method to deal with the shift-invert matrix-vector product in each inner iteration. Numerical results show that our proposed scheme is more robust and efficient than the existing high accurate implicit-explicit Euler-based extrapolation scheme.

Exponential Time Integration and Second-Order Difference Scheme for a Generalized Black-Scholes Equation

We apply an exponential time integration scheme combined with a central difference scheme on a piecewise uniform mesh with respect to the spatial variable to evaluate a generalized Black-Scholes equation. We show that the scheme is second-order convergent for both time and spatial variables. It is proved that the scheme is unconditionally stable. Numerical results support the theoretical results.