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
_______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
