shgeqz(3) on Linux
The text 5c2e17f77e24, as Debian 4.0 (Etch), Debian 3.1 (Sarge) shipped it.
implement a single-/double-shift version of the QZ method for finding the generalized eigenvalues w(j)=(ALPHAR(j) + i*ALPHAI(j))/BETAR(j) of the equation det( A - w(i) B ) = 0 In addition, the pair A,B may be reduced to generalized Schur form