\left(\frac{1}{6} \cdot {\left(-2 \cdot \log u1\right)}^{0.5}\right) \cdot \cos \left(\left(2 \cdot \pi\right) \cdot u2\right) + 0.5\frac{{\left(-2 \cdot \log u1\right)}^{0.5}}{\frac{6}{1}} \cdot \cos \left(\left(2 \cdot \pi\right) \cdot u2\right) + 0.5double code(double u1, double u2) {
return ((double) (((double) (((double) (((double) (1.0 / 6.0)) * ((double) pow(((double) (-2.0 * ((double) log(u1)))), 0.5)))) * ((double) cos(((double) (((double) (2.0 * ((double) M_PI))) * u2)))))) + 0.5));
}
double code(double u1, double u2) {
return ((double) (((double) (((double) (((double) pow(((double) (-2.0 * ((double) log(u1)))), 0.5)) / ((double) (6.0 / 1.0)))) * ((double) cos(((double) (((double) (2.0 * ((double) M_PI))) * u2)))))) + 0.5));
}



Bits error versus u1



Bits error versus u2
Results
Initial program 0.4
rmApplied clear-num0.4
Applied associate-*l/0.3
Simplified0.3
Final simplification0.3
herbie shell --seed 2020114
(FPCore (u1 u2)
:name "normal distribution"
:precision binary64
:pre (and (<= 0.0 u1 1) (<= 0.0 u2 1))
(+ (* (* (/ 1 6) (pow (* -2 (log u1)) 0.5)) (cos (* (* 2 PI) u2))) 0.5))