Penalty methods for American options with stochastic volatility

The American early exercise constraint can be viewed as transforming the original linear two dimensional stochastic volatility option pricing PDE into a PDE with a nonlinear source term. Several methods are described for enforcing the early exercise constraint by using a penalty source term in the discrete equations. The resulting nonlinear algebraic equations are solved using an approximate Newton iteration. The solution of the Jacobian is obtained using an incomplete LU (ILU) preconditioned conjugate gradient-like (PCG) method. Some example computations are presented for option pricing problems based on a stochastic volatility model, including an exotic American chooser option written on a put and call with discrete double knockout barriers and discrete dividends.

[1]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[2]  K. Vetzal.  Stochastic volatility, movements in short term interest rates, and bond option values , 1997 .

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

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

[5]  P. Forsyth,et al.  PDE methods for pricing barrier options , 2000 .

[6]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[7]  Peter A. Forsyth,et al.  A finite element approach to the pricing of discrete lookbacks with stochastic volatility , 1999 .

[8]  Dawn Hunter.  Robust numerical methods for PDE models of Asian options , 2000 .

[9]  R. LeVeque.  Numerical methods for conservation laws , 1990 .

[10]  W. K. Anderson,et al.  An implicit upwind algorithm for computing turbulent flows on unstructured grids , 1994 .

[11]  Alan G. White,et al.  The Pricing of Options on Assets with Stochastic Volatilities , 1987 .

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

[13]  G. Golub,et al.  Gmres: a Generalized Minimum Residual Algorithm for Solving , 2022 .

[14]  C. Los.  OPTION PRICING I , 2000 .

[15]  S. Turnbull,et al.  Pricing foreign currency options with stochastic volatility , 1990 .

[16]  B. Rubin,et al.  Practical control of timestep selection in thermal simulation , 1986 .

[17]  E. D'Azevedo,et al.  Towards a cost-effective ILU preconditioner with high level fill , 1992 .

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

[19]  Wei-Pai Tang,et al.  Ordering Methods for Preconditioned Conjugate Gradient Methods Applied to Unstructured Grid Problems , 1992, SIAM J. Matrix Anal. Appl..

[20]  M. Blunt,et al.  Implicit flux limiting schemes for petroleum reservoir simulation , 1992 .

[21]  Peter A. Forsyth,et al.  A Control Volume Finite Element Approach to NAPL Groundwater Contamination , 1991, SIAM J. Sci. Comput..

[22]  Mary Catherine A. Kropinski,et al.  Monotonicity Considerations for Saturated-Unsaturated Subsurface Flow , 1997, SIAM J. Sci. Comput..

[23]  Michael A. H. Dempster,et al.  Fast Numerical Valuation of American, Exotic and Complex Options , 1997 .

[24]  Peter A. Forsyth,et al.  Robust linear and nonlinear strategies for solution of the transonic Euler equations , 1995 .

[25]  J. Meijerink,et al.  An iterative solution method for linear systems of which the coefficient matrix is a symmetric -matrix , 1977 .

[26]  C. Angelopoulos.  High resolution schemes for hyperbolic conservation laws , 1992 .

[27]  P. Forsyth,et al.  Nonlinear iteration methods for high speed laminar compressible Navier-Stokes equations , 1997 .

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

[29]  W. K. Anderson,et al.  Comparison of Finite Volume Flux Vector Splittings for the Euler Equations , 1985 .