\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 r838231 = x;
double r838232 = y;
double r838233 = r838231 * r838232;
double r838234 = z;
double r838235 = r838234 * r838234;
double r838236 = 1.0;
double r838237 = r838234 + r838236;
double r838238 = r838235 * r838237;
double r838239 = r838233 / r838238;
return r838239;
}
double f(double x, double y, double z) {
double r838240 = x;
double r838241 = cbrt(r838240);
double r838242 = z;
double r838243 = cbrt(r838242);
double r838244 = r838243 * r838243;
double r838245 = r838241 / r838244;
double r838246 = r838241 / r838243;
double r838247 = r838241 / r838242;
double r838248 = y;
double r838249 = 1.0;
double r838250 = r838242 + r838249;
double r838251 = r838248 / r838250;
double r838252 = r838247 * r838251;
double r838253 = r838246 * r838252;
double r838254 = r838245 * r838253;
return r838254;
}




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