\left(a - \frac{1}{3}\right) \cdot \left(1 + \frac{1}{\sqrt{9 \cdot \left(a - \frac{1}{3}\right)}} \cdot rand\right)\left(a - \frac{1}{3}\right) \cdot \left(1 + \frac{-1}{\frac{-\sqrt{9 \cdot \left(a - \frac{1}{3}\right)}}{1 \cdot rand}}\right)double code(double a, double rand) {
return ((a - (1.0 / 3.0)) * (1.0 + ((1.0 / sqrt((9.0 * (a - (1.0 / 3.0))))) * rand)));
}
double code(double a, double rand) {
return ((a - (1.0 / 3.0)) * (1.0 + (-1.0 / (-sqrt((9.0 * (a - (1.0 / 3.0)))) / (1.0 * rand)))));
}



Bits error versus a



Bits error versus rand
Results
Initial program 0.1
rmApplied frac-2neg0.1
Applied associate-*l/0.1
rmApplied neg-mul-10.1
Applied associate-*l*0.1
Applied associate-/l*0.1
Final simplification0.1
herbie shell --seed 2020071 +o rules:numerics
(FPCore (a rand)
:name "Octave 3.8, oct_fill_randg"
:precision binary64
(* (- a (/ 1 3)) (+ 1 (* (/ 1 (sqrt (* 9 (- a (/ 1 3))))) rand))))