Numerical computation of deflating subspaces of skew-Hamiltonian/Hamiltonian pencils

Abstract

We discuss the numerical solution of structured generalized eigenvalue problems that arise from linear- quadratic optimal control problems, H infinity optimization, multibody systems, and many other areas of applied mathematics, physics, and chemistry. The classical approach for these problems requires computing invariant and deflating subspaces of matrices and matrix pencils with Hamiltonian and/ or skew- Hamiltonian structure. We extend the recently developed methods for Hamiltonian matrices to the general case of skew- Hamiltonian/ Hamiltonian pencils. The algorithms circumvent problems with skew- Hamiltonian/ Hamiltonian matrix pencils that lack structured Schur forms by embedding them into matrix pencils that always admit a structured Schur form. The rounding error analysis of the resulting algorithms is favorable. For the embedded matrix pencils, the algorithms use structure- preserving unitary matrix computations and are strongly backwards stable, i. e., they compute the exact structured Schur form of a nearby matrix pencil with the same structure.