DIDP (v0.2): A MATLAB Package to Compute the Dominant Poles of Large-Scale Descriptor Systems

These Matlab routines are meant to compute a prescribed number of dominant poles of a large-scale descriptor system.

The routines are based on an interpolatory subspace framework. At every iteration, first the dominant poles of a small projected system are computed. Then the projection subspaces are expanded so that the transfer functions of the projected system after subspace expansion and original system, as well as their derivatives, match at these computed dominant poles.