\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 r200386 = x;
double r200387 = y;
double r200388 = r200386 * r200387;
double r200389 = z;
double r200390 = r200389 * r200389;
double r200391 = 1.0;
double r200392 = r200389 + r200391;
double r200393 = r200390 * r200392;
double r200394 = r200388 / r200393;
return r200394;
}
double f(double x, double y, double z) {
double r200395 = x;
double r200396 = z;
double r200397 = r200395 / r200396;
double r200398 = y;
double r200399 = r200397 * r200398;
double r200400 = r200399 / r200396;
double r200401 = 1.0;
double r200402 = r200396 + r200401;
double r200403 = r200400 / r200402;
return r200403;
}




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))))