\frac{x \cdot y}{\left(z \cdot z\right) \cdot \left(z + 1\right)}\frac{\frac{\frac{x}{z} \cdot y}{z}}{z + 1}double f(double x, double y, double z) {
double r276052 = x;
double r276053 = y;
double r276054 = r276052 * r276053;
double r276055 = z;
double r276056 = r276055 * r276055;
double r276057 = 1.0;
double r276058 = r276055 + r276057;
double r276059 = r276056 * r276058;
double r276060 = r276054 / r276059;
return r276060;
}
double f(double x, double y, double z) {
double r276061 = x;
double r276062 = z;
double r276063 = r276061 / r276062;
double r276064 = y;
double r276065 = r276063 * r276064;
double r276066 = r276065 / r276062;
double r276067 = 1.0;
double r276068 = r276062 + r276067;
double r276069 = r276066 / r276068;
return r276069;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 15.5 |
|---|---|
| Target | 4.2 |
| Herbie | 3.3 |
Initial program 15.5
rmApplied times-frac11.4
rmApplied *-un-lft-identity11.4
Applied times-frac6.1
Applied associate-*l*2.7
rmApplied associate-*r/3.3
Applied associate-*r/3.3
Simplified3.3
Final simplification3.3
herbie shell --seed 2019326 +o rules:numerics
(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))))