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
