\frac{x \cdot y}{\left(z \cdot z\right) \cdot \left(z + 1\right)}\frac{\sqrt[3]{x}}{\sqrt[3]{z} \cdot \sqrt[3]{z}} \cdot \left(\frac{\sqrt[3]{x}}{\sqrt[3]{z}} \cdot \left(\frac{\sqrt[3]{x}}{z} \cdot \frac{y}{z + 1}\right)\right)double f(double x, double y, double z) {
double r276725 = x;
double r276726 = y;
double r276727 = r276725 * r276726;
double r276728 = z;
double r276729 = r276728 * r276728;
double r276730 = 1.0;
double r276731 = r276728 + r276730;
double r276732 = r276729 * r276731;
double r276733 = r276727 / r276732;
return r276733;
}
double f(double x, double y, double z) {
double r276734 = x;
double r276735 = cbrt(r276734);
double r276736 = z;
double r276737 = cbrt(r276736);
double r276738 = r276737 * r276737;
double r276739 = r276735 / r276738;
double r276740 = r276735 / r276737;
double r276741 = r276735 / r276736;
double r276742 = y;
double r276743 = 1.0;
double r276744 = r276736 + r276743;
double r276745 = r276742 / r276744;
double r276746 = r276741 * r276745;
double r276747 = r276740 * r276746;
double r276748 = r276739 * r276747;
return r276748;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 14.9 |
|---|---|
| Target | 4.2 |
| Herbie | 1.3 |
Initial program 14.9
rmApplied times-frac10.8
rmApplied add-cube-cbrt11.2
Applied times-frac6.3
Applied associate-*l*1.2
rmApplied add-cube-cbrt1.4
Applied times-frac1.4
Applied associate-*l*1.3
Final simplification1.3
herbie shell --seed 2020065 +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))))