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



Bits error versus a



Bits error versus rand
Results
Initial program 0.1
rmApplied sqrt-prod0.1
Applied *-un-lft-identity0.1
Applied times-frac0.2
Applied associate-*l*0.2
rmApplied distribute-lft-in0.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2020091
(FPCore (a rand)
:name "Octave 3.8, oct_fill_randg"
:precision binary64
(* (- a (/ 1 3)) (+ 1 (* (/ 1 (sqrt (* 9 (- a (/ 1 3))))) rand))))