Restricted Denominator (RD) rational approximations to the matrix exponential operator are constructed by interpolation in points related to Krylov subspaces associated to a rational transform of the particular matrix considered. Convergence analysis are provided. Numerical experiments show the effectiveness of the proposed methods in applications involving discretizations of differential operators.
