Circulant preconditioners for pricing options (Q2431151)

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.