Different eigenvectors obtained from gsl_eigen_hermv compared with scipy.linalg.eigh

[Jiasen Guo]

Dear all,
I am having an issue using gsl_eigen_hermv to reproduce a complex hermitian
eigenvalue problem previously solved with scipy.linalg.eigh. I have
included the cpp and python code in the attached zip file. Basically, the
matrix is created in cpp and solved using gsl. the same matrix (*real.csv*
for the real part and *imag.csv* for the imaginary part) is dumped out and
solved in Python again. While I get the same eigenvalues from the two
approaches, the eigenvalues are, however, different. And if I calculated
the module of the complex numbers in the eigenvector, then they are
identical from both methods.

I use g++ --std=c++11 test.cpp -o test.o -lgsl -lgslcblas for compilation.
Please help me understand what is the problem here. Thanks for your kind
help.

Best,
Jiasen

test.zip
Description: Zip compressed data