Rational Krylov methods for functions of matrices with applications to fractional partial differential equations

Abstract: In this paper, we propose a new choice of poles to define reliable rational Krylov methods. These methods are used for approximating function of positive definite matrices. In particular, the fractional power and the fractional resolvent are considered because of their importance in the numerical solution of fractional partial differential equations. The results of the numerical experiments we have carried out on some fractional models confirm that the proposed approach is promising.