| #include <math.h> | |
| #include "gsl/gsl_rng.h" | |
| int main() | |
| { | |
| // return code | |
| int retval = 1; | |
| // create a generator | |
| gsl_rng* generator; | |
| generator = gsl_rng_alloc(gsl_rng_mt19937); | |
| // Read a value. | |
| double const Result = gsl_rng_uniform(generator); | |
| // Check value | |
| double const expectedResult(0.999741748906672); | |
| if (fabs(expectedResult - Result) < 1.0e-6) | |
| retval = 0; | |
| // free allocated memory | |
| gsl_rng_free(generator); | |
| return retval; | |
| } |