Hello all,

In section 11.5 - Hermitian Matrices, it states that a Hermitian matrix is the complex analog of a real symmetric matrix .... Complex versions of the routines jacobi, tred2, and tqli are quite analogous to their real counterparts.

I have used the eigen_sym.h routine from section 11.1 to find the eigenvalues and eigenvectors of a real symmetric matrix.

My questions are:

How difficult is it to modify the jacobi routine for the complex analog?

Has anyone ever used the routines from the listed references at then end of section 11.5? If so how difficult is it to implement them?

Thanks
Matt

I'm not too sure about how to modify the Jacobi method for Hermitians, but EISPACK has routines for the similarity transformation of a Hermitian matrix to a real symmetric tridiagonal matrix. You might want to take a look to get an idea on how to implement the algorithm in C++.

Jan M.