bug-gsl
[Top][All Lists]
Advanced

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

Re: [Bug-gsl] cdf_poisson_P


From: Georgios
Subject: Re: [Bug-gsl] cdf_poisson_P
Date: Thu, 18 Oct 2012 12:37:30 +0100
User-agent: Mozilla/5.0 (X11; Linux x86_64; rv:16.0) Gecko/20121011 Thunderbird/16.0.1

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



reply via email to

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