Hi, I think the line:
j = gsl_rng_uniform_int (r, 16777216); /* sample from 2^24 */ should be j = gsl_rng_uniform_int (r, 16777216-1); /* sample from 2^24 */ as it crashes on 2^24 (exceeds upped limit for gsl_rng_uniform_int. Thanks, Gareth _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
