[Bug-gsl] [bug #47193] gsl_ran_poisson_pdf with mu=0

[Patrick Alken]

URL:
  <http://savannah.gnu.org/bugs/?47193>

                 Summary: gsl_ran_poisson_pdf with mu=0
                 Project: GNU Scientific Library
            Submitted by: psa
            Submitted on: Thu 18 Feb 2016 04:39:48 PM GMT
                Category: None
                Severity: 3 - Normal
        Operating System: 
                  Status: None
             Assigned to: None
             Open/Closed: Open
                 Release: 
         Discussion Lock: Any

    _______________________________________________________

Details:

from trojkan =at= gmail =dot= com

Hi,

I noticed that gsl_ran_poisson_pdf(k, mu) error out if mu=0.0. IMHO a Poisson
distribution with zero mean is a perfectly valid distribition (although it is
not very random), so I find it unexpected that the function doesn't support
it. I suggest to add some code taking care of the special case. Something
like:

if (!mu) {
  if (k)
    return 0.0
  return 1.0
}





    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?47193>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/