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 <[EMAIL PROTECTED]> To: [email protected] 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: [EMAIL PROTECTED]:~/projects/random/minimal_bug > gcc minimal_bug.c -lgsl -lgslcblas [EMAIL PROTECTED]:~/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 _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
