Circulant preconditioners for pricing options. (English) Zbl 1233.91317
In this article, the authors investigate a non-symmetric Toeplitz system that arises when discretizing and solving the pricing PIDE of European call options in Merton’s jump-diffusion model. If the jump distribution in Merton’s jump-diffusion model has mean \(\mu_J = 0\), then the Toeplitz system is symmetric and has been investigated in other papers. The authors look at the more general case \(\mu_J\neq 0\) and solve the resulting non-symmetric system using Strang’ circulant preconditioner in th enormalized pre-conditioned conjugate gradient method. The main result is that the spectrum is clustered around one and the smallest eigenvalue is uniformly bounded away from zero, leading to superlinear convergence of the conjugate gradient method. This is also highlighted using some numerical examples.
Reviewer: Reinhold Kainhofer (Wien)
MSC:
91G60 | Numerical methods (including Monte Carlo methods) |
91G20 | Derivative securities (option pricing, hedging, etc.) |
65F10 | Iterative numerical methods for linear systems |
47B35 | Toeplitz operators, Hankel operators, Wiener-Hopf operators |
65M06 | Finite difference methods for initial value and initial-boundary value problems involving PDEs |
Keywords:
Strang’s circulant preconditioner; nonsymmetric Toeplitz system; European call option; partial integro-differential equation; normalized preconditioned system; family of generating functionsReferences:
[1] | Almendral, A.; Oosterlee, C., Numerical valuation of options with jumps in the underlying, Appl. Numer. Math., 53, 1-18 (2005) · Zbl 1117.91028 |
[2] | Black, F.; Scholes, M., The pricing of options and corporate liabilities, J. Political Economy, 81, 637-654 (1973) · Zbl 1092.91524 |
[3] | Chan, R.; Jin, X., An Introduction to Iterative Toeplitz Solvers (2007), SIAM: SIAM Philadelphia · Zbl 1146.65028 |
[4] | Chan, R.; Ng, M., Conjugate gradient methods for Toeplitz systems, SIAM Rev., 38, 427-482 (1996) · Zbl 0863.65013 |
[5] | Chan, T., An optimal circulant preconditioner for Toeplitz systems, SIAM J. Sci. Statist. Comput., 9, 766-771 (1988) · Zbl 0646.65042 |
[6] | Cont, R.; Voltchkova, E., A finite-difference scheme for option pricing in jump-diffusion and exponential Lévy models, SIAM J. Numer. Anal., 43, 1596-1626 (2005) · Zbl 1101.47059 |
[7] | Grenander, U.; Szegö, G., Toeplitz Forms and Their Applications (1984), AMS Chelsea Publishing: AMS Chelsea Publishing Rhode Island · Zbl 0611.47018 |
[8] | Horn, R.; Johnson, C., Matrix Analysis (1985), Cambridge University Press: Cambridge University Press Cambridge · Zbl 0576.15001 |
[9] | Jin, X., Developments and Applications of Block Toeplitz Iterative Solvers (2002), Science Press/Kluwer Academic Publishers: Science Press/Kluwer Academic Publishers Beijing/Dordrecht |
[10] | Jin, X., Preconditioning Techniques for Toeplitz Systems (2010), Higher Education Press: Higher Education Press Beijing |
[11] | Matache, A.; Schwab, C.; Wihler, T., Fast numerical solution of parabolic integro-differential equations with applications in finance, SIAM J. Sci. Comput., 27, 369-393 (2005) · Zbl 1098.65123 |
[12] | Merton, R., Option pricing when underlying stock returns are discontinuous, J. Financial Econom., 3, 125-144 (1976) · Zbl 1131.91344 |
[13] | Sachs, E.; Strauss, A., Efficient solution of a partial integro-differential equation in finance, Appl. Numer. Math., 58, 1687-1703 (2008) · Zbl 1155.65109 |
[14] | Strang, G., A proposal for Toeplitz matrix calculations, Stud. Appl. Math., 74, 171-176 (1986) · Zbl 0621.65025 |
[15] | Tyrtyshnikov, E., Optimal and super-optimal circulant preconditioners, SIAM J. Matrix Anal. Appl., 13, 459-473 (1992) · Zbl 0774.65024 |
[16] | Y. Zhang, S. Vong, X. Jin, A family of generating functions with its applications, Appl. Numer. Math., submitted for publication.; Y. Zhang, S. Vong, X. Jin, A family of generating functions with its applications, Appl. Numer. Math., submitted for publication. |
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. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.