[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Bug-gsl] bug in gsl_eigen_hermv

From: Brian Gough
Subject: Re: [Bug-gsl] bug in gsl_eigen_hermv
Date: Mon, 12 Nov 2007 12:18:46 +0000
User-agent: Wanderlust/2.14.0 (Africa) Emacs/22.1 Mule/5.0 (SAKAKI)

At Sat, 10 Nov 2007 09:05:46 -0700 (MST),
Matthew B. Hastings wrote:
> I think there is a bug in gsl_eigen_hermv which causes it to give nan if
> there ar many very small entries or entries equal to zero.  I checked the
> bug archives, and I found a similar bug report,
> http://lists.gnu.org/archive/html/bug-gsl/2002-10/msg00002.html , for
> gsl-1.2.  However, I am using gsl-1.10.  The bug does not seem to appear
> on gsl-1.4.  I am using Mac OS X 10.4.9 on a PowerPC G4.

Thanks for your bug report.  

It seems to be an architecture dependent problem, I'm not able to
reproduce the NaNs on x86 or x86_64.

Is is possible for you to identify where the NaNs are coming from? e.g
by trapping any division by zero, overflow and invalid operation
exceptions while running under GDB.  The function gsl_ieee_env_setup()
can be used to enable these.

Brian Gough

Network Theory Ltd,
Publishing Free Software Manuals --- http://www.network-theory.co.uk/

reply via email to

[Prev in Thread] Current Thread [Next in Thread]