\frac{x \cdot y}{\left(z \cdot z\right) \cdot \left(z + 1\right)}\frac{\sqrt[3]{x} \cdot \sqrt[3]{x}}{z} \cdot \left(\left(\frac{\sqrt[3]{x}}{z} \cdot \left(\sqrt[3]{y} \cdot \sqrt[3]{y}\right)\right) \cdot \frac{\sqrt[3]{y}}{z + 1}\right)double f(double x, double y, double z) {
double r977630 = x;
double r977631 = y;
double r977632 = r977630 * r977631;
double r977633 = z;
double r977634 = r977633 * r977633;
double r977635 = 1.0;
double r977636 = r977633 + r977635;
double r977637 = r977634 * r977636;
double r977638 = r977632 / r977637;
return r977638;
}
double f(double x, double y, double z) {
double r977639 = x;
double r977640 = cbrt(r977639);
double r977641 = r977640 * r977640;
double r977642 = z;
double r977643 = r977641 / r977642;
double r977644 = r977640 / r977642;
double r977645 = y;
double r977646 = cbrt(r977645);
double r977647 = r977646 * r977646;
double r977648 = r977644 * r977647;
double r977649 = 1.0;
double r977650 = r977642 + r977649;
double r977651 = r977646 / r977650;
double r977652 = r977648 * r977651;
double r977653 = r977643 * r977652;
return r977653;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 15.7 |
|---|---|
| Target | 4.4 |
| Herbie | 1.4 |
Initial program 15.7
rmApplied times-frac11.8
rmApplied add-cube-cbrt12.1
Applied times-frac7.0
Applied associate-*l*1.3
rmApplied *-un-lft-identity1.3
Applied add-cube-cbrt1.5
Applied times-frac1.5
Applied associate-*r*1.4
Simplified1.4
Final simplification1.4
herbie shell --seed 2020047
(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))))