Wavelet Galerkin pricing of American options on Lévy driven assets

The price of an American style contract on assets driven by Lévy processes with infinite jump activity is expressed as solution of a parabolic variational integro-differential inequality (PIDI). A Galerkin discretization in logarithmic price using a wavelet basis is presented with compression of the moment matrix of the jump part of the price process' Dynkin operator. An iterative solver with wavelet preconditioning for the resulting large matrix inequality problems is presented and its efficiency is demonstrated by numerical experiments.