Fast numerical solution of parabolic integrodifferential equations with applications in finance

摘要

We numerically solve parabolic problems u(t) + Au = 0 in (0, T) x Omega, T < infinity, where Omega subset of R is a bounded interval and A is a strongly elliptic integrodi. erential operator of order rho is an element of [ 0, 2]. A discontinuous Galerkin (dG) discretization in time and a wavelet discretization in space are used. The densely populated matrices in the corresponding linear systems of equations are replaced by sparse ones using appropriate wavelet compression techniques. The linear systems in each time step are solved by an incomplete GMRES iteration. Under these conditions, we show that the complexity of our algorithm is linear ( up to some logarithmic terms) in the number of spatial degrees of freedom and present error estimates. Applications to purely discontinuous Levy processes arising in finance are given.