This paper presents a model reduction method for large-scale linear systems that is based on a Lanczos-type approach. A variant of the nonsymmetric Lanczos method, rational Lanczos, is shown to yield a rational interpolant (multi-point Padi approximant) for the large-scale system. An exact expression for the error in the interpolant is derived. Examples are utilized to demonstrate that the rational Lanczos method provides opportunities for significant improvements in the rate of convergence over single-point Lanczos approaches.
Gallivan, K., Grimme, E., & Van Dooren, P. (1996). A rational Lanczos algorithm for model reduction. Numerical Algorithms, 12(1-2), 33-63. https://doi.org/10.1007/BF02141740 (Original work published 1996)