Tridiagonal implicit method to evaluate European and American options under infinite activity Lévy models

We propose an efficient implicit method to evaluate European and American options when the underlying asset follows an infinite activity Lévy model. Since the Lévy measure of the infinite activity model has the singularity at the origin, we approximate infinitely many small jumps by samples of a diffusion. The proposed methods to solve partial integro–differential equations for European options and linear complementarity problems for American options via an operator splitting method involve solving linear systems with tridiagonal matrices and so can significantly reduce the computations associated with the discrete integral operators. The numerical experiments verify that the proposed method has the second-order convergence rate under an infinite activity Lévy model.