[Top][All Lists]

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

Re: [Help-gsl] Robustness/reliability of GSL

From: Brian Gough
Subject: Re: [Help-gsl] Robustness/reliability of GSL
Date: Tue, 08 Aug 2006 18:15:41 +0100
User-agent: Wanderlust/2.14.0 (Africa) Emacs/21.3 Mule/5.0 (SAKAKI)

At Sun, 6 Aug 2006 10:10:41 +0300,
David Klein wrote:
> I have used slatec and Lapack in the past, and I am thiniking of
> switching  to GSL for the thread safety. A major concern of mine is
> the robustness of the code. Numerical software is notoriously tricky
> to get right (e.g.I found a bug in cephes gamma routine where it was
> assumed that  x<1 implies that  x+1<2. As far as I know, there is
> exactly one double precision number for which this is wrong.) A big
> advantage of slatec and lapack is that they have been extensively
> tested and used. Does anyone have any feedback to provide on GSL?


Some general testing was done at CERN, see the "Reports and Papers"
section at

There is a GSL testsuite, which is fairly extensive, use "make check
GSL_TEST_VERBOSE=1" to see the full output.  Extensions are welcome.

If you want to see what bugs have been reported, the bug mailing list
archives are at

best regards, 

Brian Gough
(GSL Maintainer)

Network Theory Ltd,
Commercial support for GSL ---

reply via email to

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