\left(a - \frac{1.0}{3.0}\right) \cdot \left(1 + \frac{1}{\sqrt{9 \cdot \left(a - \frac{1.0}{3.0}\right)}} \cdot rand\right)\left(a - \frac{1.0}{3.0}\right) + \frac{rand}{\sqrt{3 \cdot \left(\left(a - \frac{1.0}{3.0}\right) \cdot 3\right)}} \cdot \left(a - \frac{1.0}{3.0}\right)double f(double a, double rand) {
double r3555727 = a;
double r3555728 = 1.0;
double r3555729 = 3.0;
double r3555730 = r3555728 / r3555729;
double r3555731 = r3555727 - r3555730;
double r3555732 = 1.0;
double r3555733 = 9.0;
double r3555734 = r3555733 * r3555731;
double r3555735 = sqrt(r3555734);
double r3555736 = r3555732 / r3555735;
double r3555737 = rand;
double r3555738 = r3555736 * r3555737;
double r3555739 = r3555732 + r3555738;
double r3555740 = r3555731 * r3555739;
return r3555740;
}
double f(double a, double rand) {
double r3555741 = a;
double r3555742 = 1.0;
double r3555743 = 3.0;
double r3555744 = r3555742 / r3555743;
double r3555745 = r3555741 - r3555744;
double r3555746 = rand;
double r3555747 = 3.0;
double r3555748 = r3555745 * r3555747;
double r3555749 = r3555747 * r3555748;
double r3555750 = sqrt(r3555749);
double r3555751 = r3555746 / r3555750;
double r3555752 = r3555751 * r3555745;
double r3555753 = r3555745 + r3555752;
return r3555753;
}



Bits error versus a



Bits error versus rand
Results
Initial program 0.1
rmApplied +-commutative0.1
Applied distribute-lft-in0.1
Simplified0.1
Simplified0.1
rmApplied add-sqr-sqrt0.1
Applied associate-*l*0.1
Final simplification0.1
herbie shell --seed 2019158
(FPCore (a rand)
:name "Octave 3.8, oct_fill_randg"
(* (- a (/ 1.0 3.0)) (+ 1 (* (/ 1 (sqrt (* 9 (- a (/ 1.0 3.0))))) rand))))