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/