[Bug-gsl] bug in gsl_ran_binomial (GSL 1.14)

[Sam Mason]

Hi,

I've just started using GSL and have been suffering with my program
locking up unexpectedly.  I've just figured out that it locks up when
running:

  #include <gsl/gsl_randist.h>
  int main()
  {
    gsl_rng_env_setup();
    gsl_rng *r = gsl_rng_alloc (gsl_rng_default);
    printf ("%i\n", gsl_ran_binomial (r, 1.0, 1u<<31));
    return 0;
  }

I.e. passing 2^31 to gsl_ran_binomial causes an infinite loop.

I'm obviously in "undefined" territory here, but I'd expect to get
garbage back (as I do when I pass 2^31+1) rather than locking up.

-- 
  Sam  http://samason.me.uk/