\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(rand \cdot \frac{\frac{1}{3}}{\sqrt{a - \frac{1.0}{3.0}}}\right) \cdot \left(a - \frac{1.0}{3.0}\right) + \left(a - \frac{1.0}{3.0}\right)double f(double a, double rand) {
double r3730924 = a;
double r3730925 = 1.0;
double r3730926 = 3.0;
double r3730927 = r3730925 / r3730926;
double r3730928 = r3730924 - r3730927;
double r3730929 = 1.0;
double r3730930 = 9.0;
double r3730931 = r3730930 * r3730928;
double r3730932 = sqrt(r3730931);
double r3730933 = r3730929 / r3730932;
double r3730934 = rand;
double r3730935 = r3730933 * r3730934;
double r3730936 = r3730929 + r3730935;
double r3730937 = r3730928 * r3730936;
return r3730937;
}
double f(double a, double rand) {
double r3730938 = rand;
double r3730939 = 0.3333333333333333;
double r3730940 = a;
double r3730941 = 1.0;
double r3730942 = 3.0;
double r3730943 = r3730941 / r3730942;
double r3730944 = r3730940 - r3730943;
double r3730945 = sqrt(r3730944);
double r3730946 = r3730939 / r3730945;
double r3730947 = r3730938 * r3730946;
double r3730948 = r3730947 * r3730944;
double r3730949 = r3730948 + r3730944;
return r3730949;
}



Bits error versus a



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