We consider interior penalty discontinuous Galerkin finite element (dGFEM) method for variable coefficient diffusion-convection-reaction equation to discretize the Heston PDE for the numerical pricing of European options. The mixed derivatives in the cross diffusion term are handled in a natural way compared to the finite difference methods. The advantages of dGFEM space discretization and Crank-Nicolson method with Rannacher smoothing as time integrator for Heston model with non-smooth initial and boundary conditions are illustrated in several numerical examples for European call, butterfly spread and digital options. The convection dominated Heston PDE for vanishing volatility is efficiently solved utilizing the adaptive dGFEM algorithm. Numerical experiments illustrate that dGFEM is highly accurate and very efficient for pricing financial options.
read more...↧