The aim is to apply the techniques used in Arbitrary Number of Stencil Points to the partial differential equation (PDE) of the Heston model with

to get – for a given accuracy – to the smallest possible lattice grids. The Heston model parameters at hand are

,

the strikes of the benchmark options are given by

and the maturity of the options should be one year. First let’s use a five point stencil in direction. As expected the order of convergence is increasing from second to fourth order.

The same experiment but now with five point stencil also in direction is shown below. Finite lattice effects show-up already at relatively small lattices due to the high convergence speed.

In time direction Richardson extrapolation can be used to increase the order of convergence from two to three.

So finally for this parameter set if the aim is to get the average pricing error below then a lattice size of

seems to be sufficient when using five point stencil operators in and direction and Richardson extrapolation in time direction. Unfortunately for a given accuracy this technique is slower than the usual suspects, e.g. Operator splitting because the matrix inversion for the Crank-Nicolson scheme has to be carried out by an iterative solver, namely by BiCGstab. The code for the numerical experiments can be found in the test case

`NthOrderDerivativeOpTest::testHigerOrderHestonOptionPricing()`

.