[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] Bug in special functions ellint.c
|
From: |
Brian Gough |
|
Subject: |
Re: [Bug-gsl] Bug in special functions ellint.c |
|
Date: |
Sat, 29 Nov 2008 08:12:33 +0000 |
|
User-agent: |
Wanderlust/2.14.0 (Africa) Emacs/22.2 Mule/5.0 (SAKAKI) |
At Sat, 29 Nov 2008 08:48:06 +1100,
Sascha Merz wrote:
> I still did not get the right result and found the bug. Line 489 is
> const double tb = -log(y) * (b[0] * y*(b[1] + y*b[2]));
> however, according to Abramowitz and Stegun, equation 17.3.34, it should be
> const double tb = -log(y) * (b[0] + y*(b[1] + y*b[2]));
> i. e. replace the second * by a +.
Thanks, this was recently fixed in the repository and will be in the
next release: http://savannah.gnu.org/bugs/?24146
--
Brian Gough
Support freedom by joining the FSF
http://www.fsf.org/news/fall-2008-fundraiser