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

The Evaluation of American Option Prices Under Stochastic Volatility and Jump-Diffusion Dynamics Using the Method of Lines

This paper considers the problem of numerically evaluating American option prices when the dynamics of the underlying are driven by both stochastic volatility following the square root process of Heston (1993), and by a Poisson jump process of the type originally introduced by Merton (1976). We develop a method of lines algorithm to evaluate the price as well as the delta and gamma of the option, thereby extending the method developed by Meyer (1998) for the case of jump-diffusion dynamics. The accuracy of the method is tested against two numerical methods that directly solve the integro-partial differential pricing equation. The first is an extension to the jump-diffusion situation of the componentwise splitting method of Ikonen & Toivanen (2007). The second method is a Crank-Nicolson scheme that is solved using projected successive over relaxation which is taken as the benchmark. The relative efficiency of these methods for computing the American call option price, delta, gamma and free boundary is analysed. If one seeks an algorithm that gives not only the price but also the delta and gamma to the same level of accuracy for a given computational effort then the method of lines seems to perform best amongst the methods considered.American options; stochastic volatility; jump-diffusion processes; Volterra integral equations; free boundary problem; method of lines

A hybrid tree/finite-difference approach for Heston-Hull-White type models

We study a hybrid tree-finite difference method which permits to obtain efficient and accurate European and American option prices in the Heston Hull-White and Heston Hull-White2d models. Moreover, as a by-product, we provide a new simulation scheme to be used for Monte Carlo evaluations. Numerical results show the reliability and the efficiency of the proposed method

Reduced Order Models for Pricing European and American Options under Stochastic Volatility and Jump-Diffusion Models

European options can be priced by solving parabolic partial(-integro) differential equations under stochastic volatility and jump-diffusion models like Heston, Merton, and Bates models. American option prices can be obtained by solving linear complementary problems (LCPs) with the same operators. A finite difference discretization leads to a so-called full order model (FOM). Reduced order models (ROMs) are derived employing proper orthogonal decomposition (POD). The early exercise constraint of American options is enforced by a penalty on subset of grid points. The presented numerical experiments demonstrate that pricing with ROMs can be orders of magnitude faster within a given model parameter variation range

Nonlinear Parabolic Equations arising in Mathematical Finance

This survey paper is focused on qualitative and numerical analyses of fully nonlinear partial differential equations of parabolic type arising in financial mathematics. The main purpose is to review various non-linear extensions of the classical Black-Scholes theory for pricing financial instruments, as well as models of stochastic dynamic portfolio optimization leading to the Hamilton-Jacobi-Bellman (HJB) equation. After suitable transformations, both problems can be represented by solutions to nonlinear parabolic equations. Qualitative analysis will be focused on issues concerning the existence and uniqueness of solutions. In the numerical part we discuss a stable finite-volume and finite difference schemes for solving fully nonlinear parabolic equations.Comment: arXiv admin note: substantial text overlap with arXiv:1603.0387