Efficient solution of backward jump-diffusion partial integro-differential equations with splitting and matrix exponentials

ABSTRACT We propose a new, unified approach to solving jump-diffusion partial integrodifferential equations (PIDEs), which often appear in mathematical finance. Our method consists of the following steps. First, a second-order operator splitting on financial processes (diffusion and jumps) is applied to these PIDEs. To solve the diffusion equation, we use standard finite-difference methods, which for multidimensional problems could also include splitting on various dimensions. For the jump part, we transform the jump integral into a pseudo-differential operator. Then, for various jump models, we show how to construct an appropriate first- and secondorder approximation on a grid that supersets the grid we used for the diffusion part. These approximations make the scheme unconditionally stable in time and preserve the positivity of the solution, which is computed either via a matrix exponential or via a Padé approximation of the matrix exponent.Various numerical experiments are provided to justify these results.