double r = 1.005; // some overall skew just to be evil
for (unsigned i = 0; i < num_samples; ++i) {
- double r_delta = 3e-3 * (rand() / (RAND_MAX+1.0) - .5);
+ double r_delta = 3e-6 * (rand() / (RAND_MAX+1.0) - .5);
r += r_delta;
fwrite(&r, sizeof(double), 1, stdout);