ABSTRACT
In the area of derivatives pricing, the key model for the theoretical evaluation of options is the Black-Scholes partial differential equation. In this paper we present a fourth order accurate discretization scheme in conjunction with Richardson extrapolation method, while for the time integration we consider high order implicit Backward Differences along with an implicit Runge-Kutta method for the numerical solution of the Black-Scholes equation in three space variables. The resulting sparse linear system is solved by the Preconditioned Biconjugate Gradient Stabilized (PBiCG-STAB) method, in conjunction with the Modified Generic Factored Approximate Sparse Inverse (MGenFAspI) scheme, based on approximate inverse sparsity patterns, using reordering schemes. Numerical results indicating the applicability along with discussions concerning the implementation issues of the proposed schemes are presented.
- Achdou, Y., and Pironneau, O. 2005. Computational methods for option pricing. SIAM. Google ScholarDigital Library
- Amestoy, P., Davis, T.A., and Duff, I.S. 1996. An approximate minimum degree ordering algorithm. SIAM Journal on Matrix Analysis and Applications 17(4), 886--905. Google ScholarDigital Library
- Chow, E. 2000. A priori sparsity patterns for parallel sparse approximate inverse preconditioners. SIAM J. Sci. Comput. 21, 1804--1822. Google ScholarDigital Library
- Duffy, D.J. 2006. Finite difference methods in financial engineering: A partial differential equation approach. John Wiley and Sons.Google Scholar
- Filelis-Papadopoulos, C.K., and Gravvanis, G.A. A class of generic factored and multilevel recursive approximate inverse techniques for solving general sparse systems (submitted).Google Scholar
- Grylonakis E-N.G. 2014. On the study and numerical solution of the Black-Scholes equation. Dissertation Thesis, Department of Electrical and Computer Engineering, Democritus University of Thrace.Google Scholar
- Gustafsson, B. 2008. High order difference methods for time dependent PDE. Springer Series in Computational Mathematics, Vol.38Google Scholar
- Haug, E.G. 2007. The complete guide to option pricing formulas. McGraw-Hill.Google Scholar
- Hull, J. 2009. Options, futures and other derivatives. Prentice-Hall.Google Scholar
- Jeong, D., Kim, J., and Wee, I-S. 2009. An accurate and efficient numerical method for Black-Scholes equations. Commun. Korean Math. Soc. 24(4), 617--628Google ScholarCross Ref
- Leentvaar, C.C.W, and Oosterlee, C.W. 2008. On coordinate transformation and grid stretching for sparse grid pricing of basket options. J. Comp. Appl. Math. 222, 193--209. Google ScholarDigital Library
- Persson, J., and von Sydow, L. 2007. Pricing European multi-asset options using a space-time adaptive FD-method. Comput. Vis. Sci. 10, 173--183. Google ScholarDigital Library
- Saad, Y., Schultz, M.H. 1986. GMRES: A generalized minimal residual algorithm for solving nonsymmetric linear systems. SIAM J. Sci. Comp. 7, 856--869. Google ScholarDigital Library
- Süli, E., and Mayers, D.F. 2003. An introduction to numerical analysis. Cambridge University Press.Google Scholar
- Trottenberg, U., Osterlee, C.W., and Schuller, A. 2000. Multigrid. Academic Press. Google ScholarDigital Library
- Van der Vorst, H.A. 1992. Bi-CGSTAB: A fast and smoothly converging variant of Bi-CG for the solution of nonsymmetric linear systems. SIAM J. Sci. and Stat. Comput. 13(2), 631--644. Google ScholarDigital Library
- Wilmott, P., Dewynne, J., and Howison, S. 1994. Option pricing: Mathematical models and computation. Oxford Financial Press.Google Scholar
- Zhang, P.G. 1998. Exotic Options: A guide to second-generation options. World Scientific, 2nd edition.Google ScholarCross Ref
Index Terms
- Higher Order Finite Difference Scheme for solving 3D Black-Scholes equation based on Generic Factored Approximate Sparse Inverse Preconditioning using Reordering Schemes
Recommendations
Robust Approximate Inverse Preconditioning for the Conjugate Gradient Method
We present a variant of the AINV factorized sparse approximate inverse algorithm which is applicable to any symmetric positive definite matrix. The new preconditioner is breakdown-free and, when used in conjunction with the conjugate gradient method, ...
Efficient and accurate finite difference schemes for solving one-dimensional Burgers' equation
In this paper, two efficient fourth-order compact finite difference algorithms have been developed to solve the one-dimensional Burgers' equation: ut+u ux=ε uxx. The methods are based on the Hopf-Cole transformation, Richardson's extrapolation, and ...
Deflation techniques in conjunction with generic factored approximate sparse inverse preconditioning
PCI '15: Proceedings of the 19th Panhellenic Conference on InformaticsA class of Modified Generic Factored Approximate Sparse Inverse matrix, based on sparsity patterns and a column-wise computational approach, has been recently proposed and used for solving general sparse linear systems. For certain classes of problems ...
Comments