A fast solver for the complex symmetric eigenproblem

Many numerical methods for studying chemical reaction problems require the computation of the eigenvalues of very large complex symetric matrices. Recently, a new algorithm for this problem has been proposed by Bar-On and Ryaboy [3]. This algorithm is similar in concept and complexity to the Hermitian eigensolver and is based on application of complex orthogonal transformations to preserve symmetry and recovery transformations to preserve stability. We demonstrate the performance of the proposed algorithm on several high performance computers from Digital, SGI, and Cray . The results show that the new algorithm is much faster than the general eigensolver, the present method used for solving these problems.