\frac{x \cdot y}{\left(z \cdot z\right) \cdot \left(z + 1.0\right)}\frac{\frac{y}{z + 1.0}}{z} \cdot \frac{x}{z}double f(double x, double y, double z) {
double r15617049 = x;
double r15617050 = y;
double r15617051 = r15617049 * r15617050;
double r15617052 = z;
double r15617053 = r15617052 * r15617052;
double r15617054 = 1.0;
double r15617055 = r15617052 + r15617054;
double r15617056 = r15617053 * r15617055;
double r15617057 = r15617051 / r15617056;
return r15617057;
}
double f(double x, double y, double z) {
double r15617058 = y;
double r15617059 = z;
double r15617060 = 1.0;
double r15617061 = r15617059 + r15617060;
double r15617062 = r15617058 / r15617061;
double r15617063 = r15617062 / r15617059;
double r15617064 = x;
double r15617065 = r15617064 / r15617059;
double r15617066 = r15617063 * r15617065;
return r15617066;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 14.5 |
|---|---|
| Target | 3.9 |
| Herbie | 3.0 |
Initial program 14.5
rmApplied times-frac11.0
rmApplied associate-/r*5.8
rmApplied div-inv5.9
Applied associate-*l*3.0
Simplified3.0
Final simplification3.0
herbie shell --seed 2019162 +o rules:numerics
(FPCore (x y z)
:name "Statistics.Distribution.Beta:$cvariance from math-functions-0.1.5.2"
:herbie-target
(if (< z 249.6182814532307) (/ (* y (/ x z)) (+ z (* z z))) (/ (* (/ (/ y z) (+ 1 z)) x) z))
(/ (* x y) (* (* z z) (+ z 1.0))))