A numerical method for computing the Hamiltonian Schur form (Q861657)

\textit{C. Paige} and \textit{C. Van Loan} [Linear Algebra Appl. 41, 11--32 (1981)] posed the open question to derive an \(O(n^3)\) method that is numerically strongly backward stable, to compute the real Hamiltonian Schur form of a Hamiltonian matrix. This problem is solved in this paper. The authors present a numerical method based on a symplectic URV-decomposition. The numerical performance of the new method indicates that if no eigenvalues of the Hamiltonian matrix are close to the imaginary axis then the method is numerically strongly backward stable.