Option valuation under the VG process by a DG method

The paper presents a discontinuous Galerkin method for solving partial integrodifferential equations arising from the European as well as American option pricing when the underlying asset follows an exponential variance gamma process. For practical purposes of numerical solving we introduce the modified option pricing problem resulting from a localization to a bounded domain and an approximation of small jumps, and we discuss the related error estimates. Then we employ a robust numerical procedure based on piecewise polynomial generally discontinuous approximations in the spatial domain. This technique enables a simple treatment of the American early exercise constraint by a direct encompassing it as an additional nonlinear source term to the governing equation. Special attention is paid to the proper discretization of non-local jump integral components, which is based on splitting integrals with respect to the domain according to the size of the jumps. Moreover, to preserve sparsity of resulting linear algebraic systems the pricing equation is integrated in the temporal variable by a semi-implicit Euler scheme. Finally, the numerical results demonstrate the capability of the numerical scheme presented within the reference benchmarks.