\frac{x \cdot y}{\left(z \cdot z\right) \cdot \left(z + 1\right)}\frac{\frac{y}{z + 1} \cdot \frac{x}{z}}{z}double f(double x, double y, double z) {
double r199722 = x;
double r199723 = y;
double r199724 = r199722 * r199723;
double r199725 = z;
double r199726 = r199725 * r199725;
double r199727 = 1.0;
double r199728 = r199725 + r199727;
double r199729 = r199726 * r199728;
double r199730 = r199724 / r199729;
return r199730;
}
double f(double x, double y, double z) {
double r199731 = y;
double r199732 = z;
double r199733 = 1.0;
double r199734 = r199732 + r199733;
double r199735 = r199731 / r199734;
double r199736 = x;
double r199737 = r199736 / r199732;
double r199738 = r199735 * r199737;
double r199739 = r199738 / r199732;
return r199739;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 15.5 |
|---|---|
| Target | 4.2 |
| Herbie | 2.6 |
Initial program 15.5
rmApplied times-frac11.4
rmApplied *-un-lft-identity11.4
Applied times-frac6.1
Applied associate-*l*2.7
rmApplied associate-*l/5.8
Applied associate-*r/5.8
Simplified2.6
Final simplification2.6
herbie shell --seed 2019326
(FPCore (x y z)
:name "Statistics.Distribution.Beta:$cvariance from math-functions-0.1.5.2"
:precision binary64
:herbie-target
(if (< z 249.6182814532307) (/ (* y (/ x z)) (+ z (* z z))) (/ (* (/ (/ y z) (+ 1 z)) x) z))
(/ (* x y) (* (* z z) (+ z 1))))