\left(a - \frac{1}{3}\right) \cdot \left(1 + \frac{1}{\sqrt{9 \cdot \left(a - \frac{1}{3}\right)}} \cdot rand\right)\frac{1 \cdot \left(a - \frac{1}{3}\right)}{\frac{\sqrt{9 \cdot \left(a - \frac{1}{3}\right)}}{rand}} + \left(a - \frac{1}{3}\right) \cdot 1double 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 (((1.0 * (a - (1.0 / 3.0))) / (sqrt((9.0 * (a - (1.0 / 3.0)))) / rand)) + ((a - (1.0 / 3.0)) * 1.0));
}



Bits error versus a



Bits error versus rand
Results
Initial program 0.1
rmApplied associate-*l/0.1
Simplified0.1
rmApplied +-commutative0.1
Applied distribute-lft-in0.1
rmApplied *-commutative0.1
Applied associate-/l*0.1
Applied associate-*r/0.1
Simplified0.1
Final simplification0.1
herbie shell --seed 2020066
(FPCore (a rand)
:name "Octave 3.8, oct_fill_randg"
:precision binary64
(* (- a (/ 1 3)) (+ 1 (* (/ 1 (sqrt (* 9 (- a (/ 1 3))))) rand))))