Search

This paper concerns the numerical solution of the two-dimensional time-dependent partial integro-differential equation (PIDE) that holds for the values of European-style options under the two-asset Kou jump-diffusion model. A main feature of this equation is the presence of a nonlocal double integral term. For its numerical evaluation, we extend a highly efficient algorithm derived by Toivanen (2008) in the case of the one-dimensional Kou integral. The acquired algorithm for the two-dimensional Kou integral has optimal computational cost: the number of basic arithmetic operations is directly proportional to the number of spatial grid points in the semidiscretization. For the effective discretization in time, we study seven contemporary operator splitting schemes of the implicit-explicit (IMEX) and the alternating direction implicit (ADI) kind. All these schemes allow for a convenient, explicit treatment of the integral term. We analyze their (von Neumann) stability. By ample numerical experiments for put-on-the-average option values, the actual convergence behavior as well as the mutual performance of the seven operator splitting schemes are investigated. Moreover, the Greeks Delta and Gamma are considered.

We study the principal component analysis based approach introduced by Reisinger & Wittum (2007) and the comonotonic approach considered by Hanbali & Linders (2019) for the approximation of American basket option values via multidimensional partial differential complementarity problems (PDCPs). Both approximation approaches require the solution of just a limited number of low-dimensional PDCPs. It is demonstrated by ample numerical experiments that they define approximations that lie close to each other. Next, an efficient discretisation of the pertinent PDCPs is presented that leads to a favourable convergence behaviour.

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.

In this note we consider splitting methods based on linear multistep methods and stabilizing corrections. To enhance the stability of the methods, we employ an idea of Bruno & Cubillos (2016) who combine a high-order extrapolation formula for the explicit term with a formula of one order lower for the implicit terms. Several examples of the obtained multistep stabilizing correction methods are presented, and results on linear stability and convergence are derived. The methods are tested in the application to the well-known Heston model arising in financial mathematics and are found to be competitive with well-established one-step splitting methods from the literature.

This paper deals with the numerical approximation of American-style option values governed by partial differential complementarity problems. For a variety of one- and two-asset American options we investigate by ample numerical experiments the temporal convergence behaviour of three modern splitting methods: the explicit payoff approach, the Ikonen-Toivanen approach and the Peaceman-Rachford method. In addition, the temporal accuracy of these splitting methods is compared to that of the penalty approach.

This paper deals with the exact calibration of semidiscretized stochastic local volatility (SLV) models to their underlying semidiscretized local volatility (LV) models. Under an SLV model, it is common to approximate the fair value of European-style options by semidiscretizing the backward Kolmogorov equation using finite differences. In the present paper we introduce an adjoint semidiscretization of the corresponding forward Kolmogorov equation. This adjoint semidiscretization is used to obtain an expression for the leverage function in the pertinent SLV model such that the approximated fair values defined by the LV and SLV models are identical for non-path-dependent European-style options. In order to employ this expression, a large non-linear system of ODEs needs to be solved. The actual numerical calibration is performed by combining ADI time stepping with an inner iteration to handle the non-linearity. Ample numerical experiments are presented that illustrate the effectiveness of the calibration procedure.

Financial derivatives pricing aims to find the fair value of a financial contract on an underlying asset. Here we consider option pricing in the partial differential equations framework. The contemporary models lead to one-dimensional or multidimensional parabolic problems of the convection-diffusion type and generalizations thereof. An overview of various operator splitting methods is presented for the efficient numerical solution of these problems. Splitting schemes of the Alternating Direction Implicit (ADI) type are discussed for multidimensional problems, e.g. given by stochastic volatility (SV) models. For jump models Implicit-Explicit (IMEX) methods are considered which efficiently treat the nonlocal jump operator. For American options an easy-to-implement operator splitting method is described for the resulting linear complementarity problems. Numerical experiments are presented to illustrate the actual stability and convergence of the splitting schemes. Here European and American put options are considered under four asset price models: the classical Black-Scholes model, the Merton jump-diffusion model, the Heston SV model, and the Bates SV model with jumps.

We consider the Modified Craig-Sneyd (MCS) scheme which forms a prominent time stepping method of the Alternating Direction Implicit type for multidimensional time-dependent convection-diffusion equations with mixed spatial derivative terms. Such equations arise often, notably, in the field of financial mathematics. In this paper a first convergence theorem for the MCS scheme is proved where the obtained bound on the global temporal discretization errors has the essential property that it is independent of the (arbitrarily small) spatial mesh width from the semidiscretization. The obtained theorem is directly pertinent to two-dimensional convection-diffusion equations with mixed derivative term. Numerical experiments are provided that illustrate our result.

In this paper a simple, effective adaptation of Alternating Direction Implicit (ADI) time discretization schemes is proposed for the numerical pricing of American-style options under the Heston model via a partial differential complementarity problem. The stability and convergence of the new methods are extensively investigated in actual, challenging applications. In addition a relevant theoretical result is proved.

In this paper we consider the stability and convergence of numerical discretizations of the Black-Scholes partial differential equation (PDE) when complemented with the popular linear boundary condition. This condition states that the second derivative of the option value vanishes when the underlying asset price gets large and is often applied in the actual numerical solution of PDEs in finance. To our knowledge, the only theoretical stability result in the literature up to now pertinent to the linear boundary condition has been obtained by Windcliff, Forsyth and Vetzal (2004) who showed that for a common discretization a necessary eigenvalue condition for stability holds. In this paper, we shall present sufficient conditions for stability and convergence when the linear boundary condition is employed. We deal with finite difference discretizations in the spatial (asset) variable and a subsequent implicit discretization in time. As a main result we prove that even though the maximum norm of exp(tM) ($t\ge 0$) can grow with the dimension of the semidiscrete matrix M, this generally does not impair the convergence behavior of the numerical discretizations. Our theoretical results are illustrated by ample numerical experiments.

In this paper the unconditional stability of four well-known ADI schemes is analyzed in the application to time-dependent multidimensional diffusion equations with mixed derivative terms. Necessary and sufficient conditions on the parameter theta of each scheme are obtained that take into account the actual size of the mixed derivative coefficients. Our results generalize results obtained previously by Craig & Sneyd (1988) and In 't Hout & Welfert (2009). Numerical experiments are presented illustrating our main theorems.