\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)\frac{a - \frac{1.0}{3.0}}{\sqrt{a - \frac{1.0}{3.0}}} \cdot \frac{rand}{3} + \left(a - \frac{1.0}{3.0}\right)double f(double a, double rand) {
double r3040864 = a;
double r3040865 = 1.0;
double r3040866 = 3.0;
double r3040867 = r3040865 / r3040866;
double r3040868 = r3040864 - r3040867;
double r3040869 = 1.0;
double r3040870 = 9.0;
double r3040871 = r3040870 * r3040868;
double r3040872 = sqrt(r3040871);
double r3040873 = r3040869 / r3040872;
double r3040874 = rand;
double r3040875 = r3040873 * r3040874;
double r3040876 = r3040869 + r3040875;
double r3040877 = r3040868 * r3040876;
return r3040877;
}
double f(double a, double rand) {
double r3040878 = a;
double r3040879 = 1.0;
double r3040880 = 3.0;
double r3040881 = r3040879 / r3040880;
double r3040882 = r3040878 - r3040881;
double r3040883 = sqrt(r3040882);
double r3040884 = r3040882 / r3040883;
double r3040885 = rand;
double r3040886 = 3.0;
double r3040887 = r3040885 / r3040886;
double r3040888 = r3040884 * r3040887;
double r3040889 = r3040888 + r3040882;
return r3040889;
}



Bits error versus a



Bits error versus rand
Results
Initial program 0.1
Simplified0.1
rmApplied sqrt-prod0.1
Applied *-un-lft-identity0.1
Applied *-un-lft-identity0.1
Applied distribute-lft-out--0.1
Applied times-frac0.2
Applied associate-*r*0.1
Simplified0.1
Final simplification0.1
herbie shell --seed 2019143
(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))))