COMPONENTWISE SPLITTING METHODS FOR PRICING AMERICAN OPTIONS UNDER STOCHASTIC VOLATILITY

Efficient numerical methods for pricing American options using Heston's stochastic volatility model are proposed. Based on this model the price of a European option can be obtained by solving a two-dimensional parabolic partial differential equation. For an American option the early exercise possibility leads to a lower bound for the price of the option. This price can be computed by solving a linear complementarity problem. The idea of operator splitting methods is to divide each time step into fractional time steps with simpler operators. This paper proposes componentwise splitting methods for solving the linear complementarity problem. The basic componentwise splitting decomposes the discretized problem into three linear complementarity problems with tridiagonal matrices. These problems can be efficiently solved using the Brennan and Schwartz algorithm, which was originally introduced for American options under the Black and Scholes model. The accuracy of the componentwise splitting method is increased by applying the Strang symmetrization. The good accuracy and the computational efficiency of the proposed symmetrized splitting method are demonstrated by numerical experiments.

[1]  Peter A. Forsyth,et al.  Convergence remedies for non-smooth payoffs in option pricing , 2003 .

[2]  A. R. Mitchell,et al.  The Finite Difference Method in Partial Differential Equations , 1980 .

[3]  H. H. Rachford,et al.  On the numerical solution of heat conduction problems in two and three space variables , 1956 .

[4]  C. M. Elliott,et al.  Weak and variational methods for moving boundary problems , 1982 .

[5]  Antonino Zanette,et al.  Parabolic ADI Methods for Pricing America Options on Two Stocks , 2002, Math. Oper. Res..

[6]  D. Duffy.  Financial Instrument Pricing Using C , 2004 .

[7]  C. Cryer.  The Solution of a Quadratic Programming Problem Using Systematic Overrelaxation , 1971 .

[8]  G. Papanicolaou,et al.  Derivatives in Financial Markets with Stochastic Volatility , 2000 .

[9]  C W Cryer,et al.  Minkowski matrices. , 1983, TOMS.

[10]  N. N. Yanenko,et al.  The Method of Fractional Steps , 1971 .

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

[12]  G. Marchuk.  Splitting and alternating direction methods , 1990 .

[13]  Eduardo S. Schwartz,et al.  The Valuation of American Put Options , 1977 .

[14]  Cornelis W. Oosterlee,et al.  On multigrid for linear complementarity problems with application to American-style options. , 2003 .

[15]  Jari Toivanen,et al.  Operator Splitting Methods for Pricing American Options with Stochastic Volatility , 2004 .

[16]  Shirley Dex,et al.  JR 旅客販売総合システム（マルス）における運用及び管理について , 1991 .

[17]  Brian J. McCartin,et al.  Accurate and efficient pricing of vanilla stock options via the Crandall-Douglas scheme , 2003, Appl. Math. Comput..

[18]  K. Manjunatha,et al.  Derivatives , 2006 .

[19]  F. Krogh,et al.  Solving Ordinary Differential Equations , 2019, Programming for Computations - Python.

[20]  Dawn Hunter.  Negative coefficients in two-factor option pricing models , 2003 .

[21]  Rüdiger U. Seydel.  Tools for Computational Finance (Universitext) , 2006 .

[22]  Sean McKee,et al.  An Alternating Direction Implicit Scheme for Parabolic Equations with Mixed Derivative and Convective Terms , 1996 .

[23]  Pricing European multi-asset options using a space–time adaptive FD-method , 2007 .

[24]  H. H. Rachford,et al.  The Numerical Solution of Parabolic and Elliptic Differential Equations , 1955 .

[25]  Kevin Parrott,et al.  Multigrid for American option pricing with stochastic volatility , 1999 .

[26]  A. Brandt,et al.  Multigrid Algorithms for the Solution of Linear Complementarity Problems Arising from Free Boundary Problems , 1983 .

[27]  Frank Cuypers.  Tools for Computational Finance , 2003 .

[28]  G. Strang.  On the Construction and Comparison of Difference Schemes , 1968 .

[29]  S. Heston.  A Closed-Form Solution for Options with Stochastic Volatility with Applications to Bond and Currency Options , 1993 .

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

[31]  Peter A. Forsyth,et al.  Penalty methods for American options with stochastic volatility , 1998 .

[32]  Peter A. Forsyth,et al.  Quadratic Convergence for Valuing American Options Using a Penalty Method , 2001, SIAM J. Sci. Comput..

[33]  Andrew R. Mitchell,et al.  Alternating Direction Methods for Parabolic Equations in Two Space Dimensions with a Mixed Derivative , 1970, Comput. J..

[34]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[35]  P. Wilmott,et al.  The Mathematics of Financial Derivatives: Preface , 1995 .

[36]  Curt Randall,et al.  Pricing Financial Instruments: The Finite Difference Method , 2000 .

[37]  R. Rannacher.  Finite element solution of diffusion problems with irregular data , 1984 .

[38]  Gabriel Wittum,et al.  On multigrid for anisotropic equations and variational inequalities “Pricing multi-dimensional European and American options” , 2004 .

[39]  D. Lamberton,et al.  Variational inequalities and the pricing of American options , 1990 .

[40]  Les Clewlow,et al.  Implementing derivatives models , 1998 .

[41]  N. N. I︠A︡nenko.  The Method of Fractional Steps: The Solution of Problems of Mathematical Physics in Several Variables , 1971 .

[42]  R. Glowinski.  Finite element methods for incompressible viscous flow , 2003 .

[43]  G. Windisch.  M-matrices in numerical analysis , 1989 .