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


_______________________________________________
Bug-gsl mailing list
[email protected]
http://lists.gnu.org/mailman/listinfo/bug-gsl

Reply via email to