Efficient L-stable method for parabolic problems with application to pricing American options under stochastic volatility

Efficient L-stable numerical method for semilinear parabolic problems with nonsmooth initial data is proposed and implemented to solve Heston's stochastic volatility model based PDE for pricing American options under stochastic volatility. The proposed new method is also used to solve two asset American options pricing problem. Cox and Matthews [S.M. Cox, P.C. Matthews, Exponential time differencing for stiff systems, Journal of Computational Physics 176 (2002) 430-455] developed a class of exponential time differencing Runge-Kutta schemes (ETDRK) for nonlinear parabolic problems. Kassam and Trefethen [A.K. Kassam, L.N. Trefethen, Fourth-order time stepping for stiff PDEs, SIAM Journal on Scientific Computing 26 (4) (2005) 1214-1233] showed that while computing certain functions involved in the Cox-Matthews schemes, severe cancelation errors can occur which affect the accuracy and stability of the schemes. Kassam and Trefethen proposed complex contour integration technique to implement these schemes in a way that avoids these cancelation errors. But this approach creates new difficulties in choosing and evaluating the contour integrals for larger problems. We modify the ETDRK schemes using positivity preserving Pade approximations of the matrix exponential functions and construct computationally efficient parallel version using splitting technique. As a result of this approach it is required only to solve several backward Euler linear problems in serial or parallel.

[1]  S. Cox,et al.  Exponential Time Differencing for Stiff Systems , 2002 .

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

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

[4]  J. Lambert.  Numerical Methods for Ordinary Differential Systems: The Initial Value Problem , 1991 .

[5]  G. Smith,et al.  Numerical Solution of Partial Differential Equations: Finite Difference Methods , 1978 .

[6]  B. A. Wade,et al.  High order smoothing schemes for inhomogeneous parabolic problems with applications in option pricing , 2007 .

[7]  Abdul-Qayyum M. Khaliq,et al.  Time-stepping algorithms for semidiscretized linear parabolic PDEs based on rational approximants with distinct real poles , 1996, Adv. Comput. Math..

[8]  Dawn Hunter,et al.  Penalty and front-fixing methods for the numerical solution of American option problems , 2002 .

[9]  E. H. Twizell,et al.  On parallel algorithms for semidiscretized parabolic partial differential equations based on subdiagonal Padé approximations , 1993 .

[10]  乔花玲,et al.  关于Semigroups of Linear Operators and Applications to Partial Differential Equations的两个注解 , 2003 .

[11]  D. A. Voss,et al.  A linearly implicit predictor–corrector scheme for pricing American options using a penalty method approach , 2006 .

[12]  Gene H. Golub,et al.  Scientific computing: an introduction with parallel computing , 1993 .

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

[14]  Lloyd N. Trefethen,et al.  Fourth-Order Time-Stepping for Stiff PDEs , 2005, SIAM J. Sci. Comput..

[15]  M. Yousuf,et al.  Pricing exotic options with L-stable Padé schemes , 2007 .

[16]  M. Yousuf,et al.  Smoothing schemes for reaction-diffusion systems with nonsmooth data , 2009 .

[17]  Paul Wilmott,et al.  Paul Wilmott on Quantitative Finance , 2010 .

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

[19]  D. Higham.  An Introduction to Financial Option Valuation: Option , 2004 .

[20]  Philip W. Livermore.  An implementation of the exponential time differencing scheme to the magnetohydrodynamic equations in a spherical shell , 2007, J. Comput. Phys..

[21]  Yousef Saad,et al.  On the parallel solution of parabolic equations , 1989, ICS '89.