[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] gsl_ran_beta can return 0 or 1
From: |
Brian Gough |
Subject: |
Re: [Bug-gsl] gsl_ran_beta can return 0 or 1 |
Date: |
Wed, 17 Nov 2004 19:41:11 +0000 |
Toby Johnson writes:
> I'm not sure whether this is a bug.
>
> gsl_ran_beta() can return values that are 0 or 1. This can cause
> problems if I then want to look up the pdf for the value just
> generated by calling gsl_ran_beta_pdf, because the pdf can be infinite
> at 0 or 1.
Thanks for the bug report, I have logged it in the BUGS file. The
problem appears to occur in gsl_ran_gamma when the parameter "a" is
small.
--
Brian Gough
Network Theory Ltd,
Commercial support for GSL --- http://www.network-theory.co.uk/gsl/