##
**Efficient numerical methods for pricing American options under stochastic volatility.**
*(English)*
Zbl 1152.91516

Summary: Five numerical methods for pricing American put options under Heston’s stochastic volatility model are described and compared. The option prices are obtained as the solution of a two-dimensional parabolic partial differential inequality. A finite difference discretization on nonuniform grids leading to linear complementarity problems with M-matrices is proposed. The projected SOR, a projected multigrid method, an operator splitting method, a penalty method, and a componentwise splitting method are considered. The last one is a direct method while all other methods are iterative. The resulting systems of linear equations in the operator splitting method and in the penalty method are solved using a multigrid method. The projected multigrid method and the componentwise splitting method lead to a sequence of linear complementarity problems with one-dimensional differential operators that are solved using the Brennan and Schwartz algorithm. The numerical experiments compare the accuracy and speed of the considered methods. The accuracies of all methods appear to be similar. Thus, the additional approximations made in the operator splitting method, in the penalty method, and in the componentwise splitting method do not increase the error essentially. The componentwise splitting method is the fastest one. All multigrid-based methods have similar rapid grid independent convergence rates. They are about two or three times slower that the componentwise splitting method. On the coarsest grid, the speed of the projected SOR is comparable with the multigrid methods while on finer grids it is several times slower.

### MSC:

91G60 | Numerical methods (including Monte Carlo methods) |

65M99 | Numerical methods for partial differential equations, initial value and time-dependent initial-boundary value problems |

65M06 | Finite difference methods for initial value and initial-boundary value problems involving PDEs |

91G20 | Derivative securities (option pricing, hedging, etc.) |

### Keywords:

American option pricing; finite difference method; linear complementarity problem; multigrid method; operator splitting method; penalty method; stochastic volatility model
PDF
BibTeX
XML
Cite

\textit{S. Ikonen} and \textit{J. Toivanen}, Numer. Methods Partial Differ. Equations 24, No. 1, 104--126 (2008; Zbl 1152.91516)

Full Text:
DOI

### References:

[1] | Black, J Polit Econ 81 pp 637– (1973) |

[2] | Dupire, Risk 7 pp 18– (1994) |

[3] | and , Financial modelling with jump processes, (Chapman & Hall/CRC financial mathematics series), Chapman & Hall/CRC, Boca Raton, FL, 2004. |

[4] | Merton, J Financ Econ 3 pp 125– (1976) |

[5] | Andersen, Rev Deriv Res 4 pp 231– (2000) |

[6] | , and , Derivatives in financial markets with stochastic volatility, Cambridge University Press, Cambridge, 2000. · Zbl 0954.91025 |

[7] | Heston, Rev Financ Stud 6 pp 327– (1993) |

[8] | Duffie, Econometrica 68 pp 1343– (2000) |

[9] | Brennan, J Finance 32 pp 449– (1977) |

[10] | Ikonen, Appl Math Lett 17 pp 809– (2004) |

[11] | and , Operator splitting methods for pricing American options with stochastic volatility, Tech. Rep. B11/2004, Department of Mathematical Information Technology, University of Jyväskylä, Jyväskylä, Finland, 2004. |

[12] | Ikonen, Int J Theor Appl Finance |

[13] | Rannacher, Numer Math 43 pp 309– (1984) |

[14] | Cryer, SIAM J Contr 9 pp 385– (1971) |

[15] | Tools for computational finance, Universitext, Springer-Verlag, Berlin, 2002. · Zbl 1044.91534 |

[16] | and , Pricing financial instruments: the finite difference method, Wiley, Chichester, 2000. |

[17] | Derivatives, Wiley, Chichester, 1998. |

[18] | Brandt, SIAM J Sci Stat Comput 4 pp 655– (1983) |

[19] | Clarke, Appl Math Finance 6 pp 177– (1999) |

[20] | Oosterlee, Electron Trans Numer Anal 15 pp 165– (2003) |

[21] | Reisinger, Comput Vis Sci 7 pp 189– (2004) · Zbl 1111.91019 |

[22] | Lions, SIAM J Numer Anal 16 pp 964– (1979) |

[23] | Zvan, J Comput Appl Math 91 pp 199– (1998) |

[24] | Forsyth, SIAM J Sci Comput 23 pp 2095– (2002) |

[25] | and , Computational methods for option pricing (Frontiers in applied mathematics), Vol. 30, SIAM, Philadelphia, PA, 2005. |

[26] | Ito, J Math Pure Appl 85 pp 415– (2006) · Zbl 1086.49004 |

[27] | and , Lagrange multiplier approach with optimized finite difference stencils for pricing American options under stochastic volatility, Tech. Rep. B6/2006, Department of Mathematical Information Technology, University of Jyväskylä, Jyväskylä, Finland, 2006. |

[28] | Nielsen, J Comput Finance 5 pp 69– (2002) |

[29] | Villeneuve, Math Oper Res 27 pp 121– (2002) |

[30] | Huang, J Comput Finance 2 pp 31– (1998) |

[31] | Pooley, J Comput Finance 6 pp 25– (2003) |

[32] | M-matrices in numerical analysis, (Teubner texts in mathematics), Vol. 115, BSB B. G. Teubner Verlagsgesellschaft, Leipzig, 1989. |

[33] | Zvan, J Comput Finance 7 pp 37– (2003) |

[34] | Kangro, SIAM J Numer Anal 38 pp 1357– (2000) |

[35] | and , The finite difference method in partial differential equations, Wiley, Chichester, 1980. |

[36] | and , Applied stochastic control of jump diffusions, Universitext, Springer-Verlag, Berlin, 2005. · Zbl 1074.93009 |

[37] | Giles, J Comput Finance 9 pp 89– (2006) |

[38] | , , , , , , , , and , Templates for the solution of linear systems: building blocks for iterative methods, 2nd ed., SIAM, Philadelphia, PA, 1994. |

[39] | Iterative solution of large linear systems, Academic Press, New York, 1971. |

[40] | , and , A multigrid tutorial, 2nd ed., SIAM, Philadelphia, PA, 2000. · Zbl 0958.65128 |

[41] | Multigrid methods and applications, (Springer series in computational mathematics), Vol. 4, Springer-Verlag, Berlin, 1985. |

[42] | , and , Multigrid, Academic Press, San Diego, CA, 2001. |

[43] | An introduction to multigrid methods, Wiley, Chichester, 1992. · Zbl 0760.65092 |

[44] | and , The multigrid solution of two-factor American put options, Tech. Rep. 96-16, Oxford Computing Laboratory, Oxford, 1996. |

[45] | and , Multi-grid methods for two phase Stefan problem, Tech. Rep. 171, TU Berlin, 1987. |

[46] | Jaillet, Acta Appl Math 21 pp 263– (1990) |

[47] | and , Weak and variational methods for moving boundary problems, (Research notes in mathematics), Vol. 59, Pitman, Boston, MA, 1982. |

[48] | Hintermüller, SIAM J Optim 13 pp 865– (2003) |

[49] | Optimization and nonsmooth analysis, Wiley, New York, 1983. · Zbl 0582.49001 |

[50] | Strang, SIAM J Numer Anal 5 pp 506– (1968) |

[51] | Finite element methods for incompressible viscous flow, Handbook of numerical analysis, Vol. IX, North-Holland, Amsterdam, 2003. |

[52] | Lötstedt, Comput Math Appl |

[53] | Cryer, ACM Trans Math Software 9 pp 199– (1983) |

This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.