Benner, P., Byers, R., & Barth, E. (2000). Algorithm 800: Fortran 77 Subroutines for Computing the Eigenvalues of Hamiltonian Matrices I: The Square-Reduced Method. ACM Transactions on Mathematical Software, 26(1), 49-77. doi:10.1145/347837.347852.