A note on the numerical approximation of Greeks for American-style options

In this note we consider the approximation of the Greeks Delta and Gamma of American-style options through the numerical solution of time-dependent partial differential complementarity problems (PDCPs). This approach is very attractive as it can yield accurate approximations to these Greeks at essentially no additional computational cost during the numerical solution of the PDCP for the pertinent option value function. For the temporal discretization, the Crank-Nicolson method is arguably the most popular method in computational finance. It is well-known, however, that this method can have an undesirable convergence behaviour in the approximation of the Greeks Delta and Gamma for American-style options, even when backward Euler damping (Rannacher smoothing) is employed. In this note we study for the temporal discretization an interesting family of diagonally implicit Runge-Kutta (DIRK) methods together with the two-stage Lobatto IIIC method. Through ample numerical experiments for one- and two-asset American-style options, it is shown that these methods can yield a regular second-order convergence behaviour for the option value as well as for the Greeks Delta and Gamma. A mutual comparison reveals that the DIRK method with suitably chosen parameter $\theta$ is preferable

Stability of central finite difference schemes for the Heston PDE

This paper deals with stability in the numerical solution of the prominent Heston partial differential equation from mathematical finance. We study the well-known central second-order finite difference discretization, which leads to large semi-discrete systems with non-normal matrices A. By employing the logarithmic spectral norm we prove practical, rigorous stability bounds. Our theoretical stability results are illustrated by ample numerical experiments

ADI schemes for valuing European options under the Bates model

This paper is concerned with the adaptation of alternating direction implicit (ADI) time discretization schemes for the numerical solution of partial integro-differential equations (PIDEs) with application to the Bates model in finance. Three different adaptations are formulated and their (von Neumann) stability is analyzed. Ample numerical experiments are provided for the Bates PIDE, illustrating the actual stability and convergence behaviour of the three adaptations.peerReviewe

ADI finite difference schemes for the Heston-Hull-White PDE

In this paper we investigate the effectiveness of Alternating Direction Implicit (ADI) time discretization schemes in the numerical solution of the three-dimensional Heston-Hull-White partial differential equation, which is semidiscretized by applying finite difference schemes on nonuniform spatial grids. We consider the Heston-Hull-White model with arbitrary correlation factors, with time-dependent mean-reversion levels, with short and long maturities, for cases where the Feller condition is satisfied and for cases where it is not. In addition, both European-style call options and up-and-out call options are considered. It is shown through extensive tests that ADI schemes, with a proper choice of their parameters, perform very well in all situations - in terms of stability, accuracy and efficiency.