[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Bug-gsl] NaN elements returned by gsl_eigen_hermv()
From: |
Beng Heng, Ng |
Subject: |
[Bug-gsl] NaN elements returned by gsl_eigen_hermv() |
Date: |
Sun, 07 Mar 2010 00:33:38 -0500 |
User-agent: |
Mozilla/5.0 (Windows; U; Windows NT 6.0; en-US; rv:1.9.1.7) Gecko/20100111 Thunderbird/3.0.1 |
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] NaN elements returned by gsl_eigen_hermv(),
Beng Heng, Ng <=