[Bug-gsl] Bug in gsl_ran_chisq()

[Pascal Vaudrevange]

Hi,

when using the random number generator gsl_rng_ranmar, a call to the 
distribution gsl_ran_chisq() results in an error. Using the generator 
gsl_rng_default works alright.

Sample code:
--- snip ---

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

int main(int argc, char * argv[]){
  gsl_rng * ran_set;
  const gsl_rng_type *gsltype = gsl_rng_ranmar;
  int seed=0;
  double x;

  gsl_rng_env_setup();
  ran_set=gsl_rng_alloc(gsltype);

  gsl_rng_set(ran_set, seed);
  x=gsl_ran_chisq(ran_set, 13.0);
  return 0;
}
--- snip ---

Output:

address@hidden:~/projects/random/minimal_bug > gcc minimal_bug.c -lgsl 
-lgslcblas
address@hidden:~/projects/random/minimal_bug > ./a.out
gsl: ../gsl/gsl_rng.h:201: ERROR: invalid n, either 0 or exceeds maximum 
value of generator
Default GSL error handler invoked.
Aborted (core dumped)


Thanks,

Pascal