[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Bug-gsl] [bug #47646] gsl_ran_beta returns NaN for small arguments

From: Yu Liu
Subject: [Bug-gsl] [bug #47646] gsl_ran_beta returns NaN for small arguments
Date: Fri, 15 Jul 2016 21:46:49 +0000 (UTC)
User-agent: Mozilla/5.0 (Windows NT 6.1; WOW64; rv:47.0) Gecko/20100101 Firefox/47.0

Follow-up Comment #2, bug #47646 (project gsl):

Following up based on numpy issues
(please let me know if this kind of information should be sent over the
mailing list.)

numpy fixed this problem around March 2015
It was related to a Dirichlett distribution bug that is still open.

As for the code, the while loop is there to handle the very RARE cases of `u`
and `v` being zero, because `rk_double` could give a zero. But with
`gsl_ran_uniform_pos` this loop is redundant. 

In terms of test, `numpy` code only check if there is `nan` generated. See
numpy #5858 above.


Reply to this item at:


  Message sent via/by Savannah

reply via email to

[Prev in Thread] Current Thread [Next in Thread]