\frac{x \cdot y}{\left(z \cdot z\right) \cdot \left(z + 1\right)}\frac{\frac{x}{z} \cdot \frac{y}{z + 1}}{z}double f(double x, double y, double z) {
double r177499 = x;
double r177500 = y;
double r177501 = r177499 * r177500;
double r177502 = z;
double r177503 = r177502 * r177502;
double r177504 = 1.0;
double r177505 = r177502 + r177504;
double r177506 = r177503 * r177505;
double r177507 = r177501 / r177506;
return r177507;
}
double f(double x, double y, double z) {
double r177508 = x;
double r177509 = z;
double r177510 = r177508 / r177509;
double r177511 = y;
double r177512 = 1.0;
double r177513 = r177509 + r177512;
double r177514 = r177511 / r177513;
double r177515 = r177510 * r177514;
double r177516 = r177515 / r177509;
return r177516;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 14.8 |
|---|---|
| Target | 4.1 |
| Herbie | 2.4 |
Initial program 14.8
rmApplied times-frac10.9
rmApplied *-un-lft-identity10.9
Applied times-frac6.0
Applied associate-*l*2.5
rmApplied *-un-lft-identity2.5
Applied add-sqr-sqrt2.5
Applied times-frac2.5
Applied associate-*l*2.5
Simplified2.4
Final simplification2.4
herbie shell --seed 2019304 +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.618281453230708) (/ (* y (/ x z)) (+ z (* z z))) (/ (* (/ (/ y z) (+ 1 z)) x) z))
(/ (* x y) (* (* z z) (+ z 1))))