In this paper, the numerical evaluation of matrix functions expressed in partial fraction form is addressed.
The shift-and-invert Krylov method is analyzed, with special attention to error estimates. Such estimates
give insights into the selection of the shift parameter and lead to a simple and effective restart
procedure. Applications to the class of Mittag–Leffler functions are presented.
