Pricing European and American Options under Heston Model using Discontinuous Galerkin Finite Elements

This paper deals with pricing of European and American options, when the underlying asset price follows Heston model, via the interior penalty discontinuous Galerkin finite element method (dGFEM). The advantages of dGFEM space discretization with Rannacher smoothing as time integrator with nonsmooth initial and boundary conditions are illustrated for European vanilla options, digital call and American put options. The convection dominated Heston model for vanishing volatility is efficiently solved utilizing the adaptive dGFEM. For fast solution of the linear complementary problem of the American options, a projected successive over relaxation (PSOR) method is developed with the norm preconditioned dGFEM. We show the efficiency and accuracy of dGFEM for option pricing by conducting comparison analysis with other methods and numerical experiments

Application of Operator Splitting Methods in Finance

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

B-Spline Based Methods: From Monotone Multigrid Schemes for American Options to Uncertain Volatility Models

In the first part of this thesis, we consider B-spline based methods for pricing American options in the Black-Scholes and Heston model. The difference between these two models is the assumption on the volatility of the underlying asset. While in the Black-Scholes model the volatility is assumed to be constant, the Heston model includes a stochastic volatility variable. The underlying problems are formulated as parabolic variational inequalities. Recall that, in finance, to determine optimal risk strategies, one is not only interested in the solution of the variational inequality, i.e., the option price, but also in its partial derivatives up to order two, the so-called Greeks. A special feature for these option price problems is that initial conditions are typically given as piecewise linear continuous functions. Consequently, we have derived a spatial discretization based on cubic B-splines with coinciding knots at the points where the initial condition is not differentiable. Together with an implicit time stepping scheme, this enables us to achieve an accurate pointwise approximation of the partial derivatives up to order two. For the efficient numerical solution of the discrete variational inequality, we propose a monotone multigrid method for (tensor product) B-splines with possible internal coinciding knots. Corresponding numerical results show that the monotone multigrid method is robust with respect to the refinement level and mesh size. In the second part of this thesis, we consider the pricing of a European option in the uncertain volatility model. In this model the volatility of the underlying asset is a priori unknown and is assumed to lie within a range of extreme values. Mathematically, this problem can be formulated as a one dimensional parabolic Hamilton-Jacobi-Bellman equation and is also called Black-Scholes-Barenblatt equation. In the resulting non-linear equation, the diffusion coefficient is given by a volatility function which depends pointwise on the second derivative. This kind of non-linear partial differential equation does not admit a weak H^1-formulation. This is due to the fact that the non-linearity depends pointwise on the second derivative of the solution and, thus, no integration by parts is possible to pass the partial derivative onto a test function. But in the discrete setting this pointwise second derivative can be approximated in H^1 by L^1-normalized B-splines. It turns out that the approximation of the volatility function leads to discontinuities in the partial derivatives. In order to improve the approximation of the solution and its partial derivatives for cubic B-splines, we develop a Newton like algorithm within a knot insertion step. Corresponding numerical results show that the convergence of the solution and its partial derivatives are nearly optimal in the L^2-norm, when the location of volatility change is approximated with desired accuracy

Analytic approach to solve a degenerate parabolic PDE for the Heston Model

We present an analytic approach to solve a degenerate parabolic problem associated to the Heston model, which is widely used in mathematical finance to derive the price of an European option on an risky asset with stochastic volatility. We give a variational formulation, involving weighted Sobolev spaces, of the second order degenerate elliptic operator of the parabolic PDE. We use this approach to prove, under appropriate assumptions on some involved unknown parameters, the existence and uniqueness of weak solutions to the parabolic problem on unbounded subdomains of the half-plane