Abstract

We present a fraction-free approach to the computation of matrix Padé systems. The method relies on determining a modified Schur complement for the coefficient matrices of the linear systems of equations that are associated to matrix Padé approximation problems. By using this modified Schur complement for these matrices we are able to obtain a hybrid fraction-free algorithm for their computation.

The algorithm that is presented is general and requires no extra assumptions on its input. The algorithm is also fast in the sense that it is usually an order of magnitute faster than existing fraction-free methods for the corresponding problem.