Revisiting IRKA: Connections with Pole Placement and Backward Stability

Christopher Beattie,Zlatko Drmač,Serkan Gugercin
DOI: https://doi.org/10.1007/s10013-020-00424-0
2020-08-05
Vietnam Journal of Mathematics
Abstract:The iterative rational Krylov algorithm (<span class="u-sans-serif">IRKA</span>) is a popular approach for producing locally optimal reduced-order <span class="mathjax-tex">\({\mathscr{H}}_{2}\)</span>-approximations to linear time-invariant (LTI) dynamical systems. Overall, <span class="u-sans-serif">IRKA</span> has seen significant practical success in computing high fidelity (locally) optimal reduced models and has been successfully applied in a variety of large-scale settings. Moreover, <span class="u-sans-serif">IRKA</span> has provided a foundation for recent extensions to the systematic model reduction of bilinear and nonlinear dynamical systems. Convergence of the basic <span class="u-sans-serif">IRKA</span> iteration is generally observed to be rapid—but not always; and despite the simplicity of the iteration, its convergence behavior is remarkably complex and not well understood aside from a few special cases. The overall effectiveness and computational robustness of the basic <span class="u-sans-serif">IRKA</span> iteration is surprising since its algorithmic goals are very similar to a pole assignment problem, which can be notoriously ill-conditioned. We investigate this connection here and discuss a variety of nice properties of the <span class="u-sans-serif">IRKA</span> iteration that are revealed when the iteration is framed with respect to a primitive basis. We find that the connection with pole assignment suggests refinements to the basic algorithm that can improve convergence behavior, leading also to new choices for termination criteria that assure backward stability.
English Else
What problem does this paper attempt to address?