Primal-dual active set method for pricing American better-of option on two assets

Abstract An American better-of option is a rainbow option with two underlying assets. It can be described by a parabolic variational inequality (VI) on a two-dimensional unbounded domain, which can also be characterized by a two-dimensional free boundary problem. Based on the numeraire transformation and the known information on the free boundary, we derive a one-dimensional linear complementarity problem (LCP) related to options on a bounded domain. Moreover, the full discretization scheme of LCP is constructed by finite difference and finite element methods in temporal and spatial directions, respectively. The primal-dual active-set (PDAS) method is adopted for solving the resulting large-scale discretized system. In each PDAS iteration, a unique index set of primal-dual variables will be computed by solving a linear system. A systematical convergent analysis is presented on our proposed method for pricing American better-of options. One of the desirable features of our method is that we can get the option value and two free boundaries simultaneously. Numerical simulations are performed to verify the efficiency of our proposed method.

[1]  F. Black,et al.  The Pricing of Options and Corporate Liabilities , 1973, Journal of Political Economy.

[2]  Stephen J. Wright,et al.  Numerical Optimization , 2018, Fundamental Statistical Inference.

[3]  Y. Kwok.  Mathematical models of financial derivatives , 2008 .

[4]  Nassim Nicholas Taleb,et al.  Dynamic Hedging: Managing Vanilla and Exotic Options , 1997 .

[5]  S. Ross,et al.  Option pricing: A simplified approach☆ , 1979 .

[6]  Kazufumi Ito,et al.  Augmented Lagrangian methods for nonsmooth, convex optimization in Hilbert spaces , 2000 .

[7]  Claes Johnson,et al.  A Convergence Estimate for an Approximation of a Parabolic Variational Inequality , 1976 .

[8]  Kazufumi Ito,et al.  Optimal Control of Elliptic Variational Inequalities , 2000 .

[9]  Sunil Kumar,et al.  MULTIDIMENSIONAL PORTFOLIO OPTIMIZATION WITH PROPORTIONAL TRANSACTION COSTS , 2006 .

[10]  J. Lions,et al.  Non-homogeneous boundary value problems and applications , 1972 .

[11]  René M. Stulz,et al.  Options on the minimum or the maximum of two risky assets : Analysis and applications , 1982 .

[13]  L. Rogers.  Monte Carlo valuation of American options , 2002 .

[14]  Kazufumi Ito,et al.  The Primal-Dual Active Set Strategy as a Semismooth Newton Method , 2002, SIAM J. Optim..

[15]  J. Hull.  Options, Futures, and Other Derivatives , 1989 .

[16]  Christian Kanzow,et al.  Inexact semismooth Newton methods for large-scale complementarity problems , 2004, Optim. Methods Softw..

[17]  C. Vuik,et al.  AnL2-error estimate for an approximation of the solution of a parabolic variational inequality , 1990 .

[18]  K. Kunisch,et al.  Primal-Dual Strategy for Constrained Optimal Control Problems , 1999 .

[19]  Lishang Jiang.  Mathematical Modeling and Methods of Option Pricing , 2005 .

[20]  Haiming Song,et al.  Front-fixing FEMs for the pricing of American options based on a PML technique , 2015 .

[21]  Walter Allegretto,et al.  Finite Element Error Estimates for a Nonlocal Problem in American Option Valuation , 2001, SIAM J. Numer. Anal..

[22]  M. Broadie,et al.  The Valuation of American Options on Multiple Assets , 1997 .

[23]  Xiaonan Wu,et al.  A Fast Numerical Method for the Black-Scholes Equation of American Options , 2003, SIAM J. Numer. Anal..

[24]  Haiming Song,et al.  A fast numerical method for the valuation of American lookback put options , 2015, Commun. Nonlinear Sci. Numer. Simul..

[25]  Gilbert Strang,et al.  Approximation in the finite element method , 1972 .