Abstract
In this paper, we consider numerical pricing of European and American options under the Bates model, a model which gives rise to a partial-integro differential equation. This equation is discretized in space using adaptive finite differences while an IMEX scheme is employed in time. The sparse linear systems of equations in each time-step are solved using an LU-decomposition and an operator splitting technique is employed for the linear complementarity problems arising for American options. The integral part of the equation is treated explicitly in time which means that we have to perform matrix-vector multiplications each time-step with a matrix with dense blocks. These multiplications are accomplished through fast Fourier transforms. The great performance of the method is demonstrated through numerical experiments.
Acknowledgments
The work in this paper builds on the sequence of Master theses/project reports [Citation31,Citation32,Citation37,Citation38]. The authors are thankful for the early contributions of these students. The computations were performed on resources provided by the Swedish National Infrastructure for Computing (SNIC) through Uppsala Multidisciplinary Center for Advanced Computational Science (UPPMAX) under Projects snic2014-3-24 and snic2015-6-77. We thank the referees for their suggestions and comments.
Disclosure statement
No potential conflict of interest was reported by the authors.