Mixing monte-carlo and partial differential equations for pricing options

Abstract

There is a need for very fast option pricers when the financial objects are modeled by complex systems of stochastic differential equations. Here the authors investigate option pricers based on mixed Monte-Carlo partial differential solvers for stochastic volatility models such as Heston’s. It is found that orders of magnitude in speed are gained on full Monte-Carlo algorithms by solving all equations but one by a Monte-Carlo method, and pricing the underlying asset by a partial differential equation with random coefficients, derived by Itô calculus. This strategy is investigated for vanilla options, barrier options and American options with stochastic volatilities and jumps optionally.

Keywords

Monte-Carlo / Partial differential equations / Heston model / Financial mathematics / Option pricing

References

[1] Achdou Y, Pironneau O. Numerical Methods for Option Pricing, 2005, Philadelphia, USA: SIAM series

[2] Amin K, Khanna A. Convergence of American option values from discrete- to continuous-time financial models. Mathematical Finance, 1994, 4: 289-304

[3] Barth A, Schwab C, Zollinger N. Multi-level Monte Carlo finite element method for elliptic PDEs with stochastic coefficients. Num. Math., 2011, 119(1): 123-161

[4] Bates D S. Jumps and stochastic volatility: Exchange rate processes implicit Deutsche mark options. Review Financial Stud., 1996, 9(1): 69-107

[5] Black F, Scholes M. The pricing of options and coorporate liabilities. J. Political Econ., 1973, 81: 637-659

[6] Boyle P. Options: A Monte Carlo approach. Journal of Financial Economics, 1977, 4: 323-338

[7] Carr P, Madan D. Option valuation using the fast Fourier transform. Journal of Computational Finance, 1999, 2(4): 61-73

[8] Dupire B. Pricing with a smile. Risk, 1994 18-20

[9] George P L, Borouchaki H. Delaunay triangulation and meshing, 1998

[10] Glasserman P. Monte-Carlo Methods in Financial Engineering, Stochastic Modeling and Applied Probability, Vol. 53, 2004, New York: Springer-Verlag

[11] Hecht, F., Pironneau, O., Le Yaric, A., et al., freefem++ documentation. http://www.freefem.org

[12] Heston S. A closed form solution for options with stochastic volatility with application to bond and currency options. Review with Financial Studies, 1993, 6(2): 327-343

[13] Lewis A. A simple option formula for general jump-diffusion and other exponential Lévy processes, 2001

[14] Loeper G, Pironneau O. A mixed PDE/Monte-Carlo method for stochastic volatility models. C. R. Acad. Sci. Paris, Ser. I, 2009, 347: 559-563

[15] Merton R C. Option pricing when underlying stock returns are discontinuous. J. Financ. Econ., 1976, 3: 125-144

[16] Pironneau, O., Dupire Identities for Complex Options, Comptes Rendus de l’Académie des Sciences, to appear.

[17] Li, X. S., Demmel, J. W. and Gilbert, J. R., The superLU library. http://crd.lbl.gov/xiaoye/SuperLU

[18] Wilmott P, Howison S, Dewynne J. The Mathematics of Financial Derivatives, 1995, Cambridge: Cambridge University Press