Krylov subspace methods for functions of fractional differential operators

The paper deals with the computation of functions of fractional powers of differential operators. The spectral properties of these operators naturally suggest the use of rational approximations. In this view we analyze the convergence properties of the shift-and-invert Krylov method applied to operator functions arising from the numerical solution of differential equations involving fractional diffusion.