[Bug-gsl] NaN elements returned by gsl_eigen_hermv()

[Beng Heng, Ng]

Hi,

I believe there is a bug in gsl_eigen_hermv for at least one Hermitian 
matrix that results in NaN values in the majority of the elements of 
evec. There are some 0 elements, though just a few of them.

There was a fix in GSL v1.13 for issues with NAN values. However, I'm 
already using GSL v1.13.

Since it is a 264x264 matrix, I've attached a sample program that reads 
in the matrix from a file. I compiled it as follows:
    g++ -lgsl -lgslcblas hermvbug.cpp -o hermvbug.exe

It will read in the matrix, check that the input matrix is indeed 
Hermitian, call gsl_eigen_hermv, and finally write the output complex 
eigenvectors, evec, to "bugoutput.txt".

I'm using GSL v1.13, gcc v4.3.2, Linux kernel 2.6.26-2-amd64. The 
processor is Intel Q9300.

I'll be glad to provide more information, and will definitely appreciate 
any feedback regarding this issue.

Thanks!

Cheers,
Beng Heng

hermvbug.tgz
Description: Binary data