[Fwd: Bug in gsl_ran_chisq()]

[Pascal Vaudrevange <pascal.vaudrevange@........>]

Hi,

sorry, I forgot: the version number is 1.9.3.

Pascal

-------- Original Message --------
Subject: Bug in gsl_ran_chisq()
Date: Wed, 27 Feb 2008 15:01:48 -0500
From: Pascal Vaudrevange <pascal.vaudrevange@........>
To: bug-gsl@gnu.org

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:

pascal@octonion:~/projects/random/minimal_bug > gcc minimal_bug.c -lgsl
-lgslcblas
pascal@octonion:~/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