double gauss_diff(){ static pfi::math::random::random<mersenne_twister> r; static double last_val = r.next_gaussian(); double cur_val = r.next_gaussian(); double ret = (last_val-cur_val)/sqrt(2.0); last_val=cur_val; return ret; }
double gauss_shifted(){ static pfi::math::random::random<mersenne_twister> r; return (r.next_gaussian(42.0, 24.0)-42.0)/24.0; }
double gauss(){ static pfi::math::random::random<mersenne_twister> r; return r.next_gaussian(); }