Follow-up Comment #2, bug #43326 (group gsl):

First order patch in this case means:
1- test cases not provided, although something like the sample code below
proves the result
2- solution does not work in conjunction with histograms, as this will yield a
xmin larger than xmax error; however, given that the case is degenerate, more
than likely other error messages will be presented too.

#include <stdio.h>
#include <gsl/gsl_rng.h>
#include <gsl/gsl_randist.h>

int
main (void)
{
  for (int k = 3; k >= 0; k--)
    for (double mu_ = 1; mu_ >= 0; mu_ -= 0.1)
      printf("%d\t%.1f\t%f\n", k, mu_, gsl_ran_poisson_pdf(k, mu_));

  return 0;
}



    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/bugs/?43326>

_______________________________________________
Message sent via Savannah
https://savannah.gnu.org/

Attachment: signature.asc
Description: PGP signature

Reply via email to