Just to clarify: the function I am talking about is *gsl_cdf_poisson_P*,
Not gsl_ran_poisson_pdf that I mistakenly mentioned below. The function
*gsl_ran_poisson_pdf* (unsigned int k, double mu) is fine, and
*gsl_ran_poisson_pdf*(-1,5) returns the expected zero.
The answer to your question is no, I do not see any warning messages. In
the past the compiler produced warning messages for other cases, so the
warning must be on.
On 17/10/12 22:01, Rhys Ulerich wrote:
This time with the list CCed... Sorry I keep screwing up on the CC.
I think there is an issue with function *gsl_ran_poisson_pdf* (unsigned int
k, double mu):
gsl_cdf_poisson_P(-1,5) returns 1.0000
while
gsl_cdf_poisson_P(1,5) returns 0.040428
The second number is fine, but the first one should be either equal to the
second one, or if the sign makes a difference, it should be zero.
If you turn on warnings, does the compiler shout at you for providing
a negative number (-1) as an unsigned argument (unsigned int)?
- Rhys
