[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Bug-gsl] [Fwd: Bug in gsl_ran_chisq()]
From: |
Pascal Vaudrevange |
Subject: |
[Bug-gsl] [Fwd: Bug in gsl_ran_chisq()] |
Date: |
Wed, 27 Feb 2008 15:11:36 -0500 |
User-agent: |
Thunderbird 2.0.0.6 (X11/20071022) |
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 <address@hidden>
To: address@hidden
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
[Prev in Thread] |
Current Thread |
[Next in Thread] |
- [Bug-gsl] [Fwd: Bug in gsl_ran_chisq()],
Pascal Vaudrevange <=