Diagonalization of complex symmetric matrices: Generalized Householder reflections, iterative deflation and implicit shifts

For simplicity, the “wrapper” program which contains an example application and the HTDQLS routine are distributed in the same file.