In this paper, we consider the fast numerical valuation of European and American options under Merton’s jump-diffusion model, which is given by a partial integro-differential equations. Due to the singularities and discontinuities of the model, the time-space grids are nonuniform with refinement near the strike price and expiry. On such nonuniform grids, the spatial differential operators are discretized by finite difference methods, and time stepping is performed using the discontinuous Galerkin finite element method. Owing to the nonuniform ...
