[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] Bug in normalized incomplete Gamma function/upper tail Gam
Re: [Bug-gsl] Bug in normalized incomplete Gamma function/upper tail Gamma dist'n probs.
Wed, 7 Sep 2005 11:47:37 +0100
Brian Williams writes:
> Version number of GSL: 1.6
> Hardware & OS: Dell PowerEdge 6850, Intel Xeon processor at 3.0GHz/8MB
> Cache, Linux 2.4.21-27.EL
> Compiler: gcc 3.2.3 20030502 (Red Hat Linux 3.2.3-52)
> gcc -o test.out -Wall -lgsl -lgslcblas -lm test.c
> Description: For some inputs, answers are returned that are incorrect.
> Here is a test program (test.c):
Thanks for the bug report. This will be fixed in the next release,
which should be available in the next week or so. The cause was an
instability in a continued fraction, you can find the corrected code
in cvs at http://sources.redhat.com/gsl/devel.html (specfunc/gamma_inc.c).
Network Theory Ltd,
Publishing Free Software Manuals --- http://www.network-theory.co.uk/