[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
https://github.com/numpy/numpy/pull/5858
It was related to a Dirichlett distribution bug that is still open.
https://github.com/numpy/numpy/issues/5851
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:
<http://savannah.gnu.org/bugs/?47646>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/