[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: Patrick Alken
Subject: [Bug-gsl] [bug #47646] gsl_ran_beta returns NaN for small arguments
Date: Wed, 1 Nov 2017 05:14:53 -0400 (EDT)
User-agent: Mozilla/5.0 (X11; Linux x86_64) AppleWebKit/537.36 (KHTML, like Gecko) Chrome/61.0.3163.100 Safari/537.36

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

[From Yu]

   I was thinking if it's worthwhile to allow the use of close-form cdf for
some functions. Currently, the test routine calculates the expected
probability in each bin using numerical integration of pdf. In this particular
case, it blows up. Using cdf in stead of integration can avoid the problem,
but I am not sure if this is appropriate. There are other functions with
close-form cdf as well, and I am wondering whether it can be used to replace
integration as well.
  In a modified test routine, I added a flag to indicate whether the input
function ``pdf" is actually a cdf, and if that is the case, the probability in
each bin would be simply cdf[x+dx] - cdf[x]. Please see the relevant lines
One obvious issue is that in this function, the input ``pdf" could mean pdf or
cdf depending on the value of the flag. This may be confusing.
   Without much experience in numerical computation, maybe I am proposing
something crazy.



Reply to this item at:


  Message sent via/by Savannah

reply via email to

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