[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] Erroneous values from gsl_ran_chisq_pdf
From: |
Brian Gough |
Subject: |
Re: [Bug-gsl] Erroneous values from gsl_ran_chisq_pdf |
Date: |
Sat, 20 Nov 2010 12:05:38 +0000 |
User-agent: |
Wanderlust/2.15.6 (Almost Unreal) Emacs/23.2 Mule/6.0 (HANACHIRUSATO) |
At Wed, 10 Nov 2010 14:37:33 +0100,
Teemu Ikonen wrote:
> The current implementation of the chi-squared distribution pdf in GSL
> returns 0.0 for inputs with x<=0.0.
>
> The domain of this function is [0.0, inf), (see e.g. Abramovitz-Stegun
> here: http://people.math.sfu.ca/~cbm/aands/page_940.htm ),
> so wrong values at x=0.0 are returned for nu=1.0
> (gsl_ran_chisq_pdf(0.0, 1.0) should be inf) and nu=2.0 (should be
> 0.5).
>
> Attached is a patch which fixes this.
Thanks for the patch, I have committed it to the repository.
--
Brian Gough