\frac{x \cdot y}{\left(z \cdot z\right) \cdot \left(z + 1\right)}\left(\sqrt[3]{\frac{\sqrt[3]{x} \cdot \sqrt[3]{x}}{z} \cdot \left(\frac{\sqrt[3]{x}}{z} \cdot \frac{y}{z + 1}\right)} \cdot \sqrt[3]{\frac{\sqrt[3]{x} \cdot \sqrt[3]{x}}{z} \cdot \left(\frac{\sqrt[3]{x}}{z} \cdot \frac{y}{z + 1}\right)}\right) \cdot \sqrt[3]{\frac{\sqrt[3]{x} \cdot \sqrt[3]{x}}{z} \cdot \left(\frac{\sqrt[3]{x}}{z} \cdot \frac{y}{z + 1}\right)}double f(double x, double y, double z) {
double r354274 = x;
double r354275 = y;
double r354276 = r354274 * r354275;
double r354277 = z;
double r354278 = r354277 * r354277;
double r354279 = 1.0;
double r354280 = r354277 + r354279;
double r354281 = r354278 * r354280;
double r354282 = r354276 / r354281;
return r354282;
}
double f(double x, double y, double z) {
double r354283 = x;
double r354284 = cbrt(r354283);
double r354285 = r354284 * r354284;
double r354286 = z;
double r354287 = r354285 / r354286;
double r354288 = r354284 / r354286;
double r354289 = y;
double r354290 = 1.0;
double r354291 = r354286 + r354290;
double r354292 = r354289 / r354291;
double r354293 = r354288 * r354292;
double r354294 = r354287 * r354293;
double r354295 = cbrt(r354294);
double r354296 = r354295 * r354295;
double r354297 = r354296 * r354295;
return r354297;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 14.6 |
|---|---|
| Target | 4.1 |
| Herbie | 1.2 |
Initial program 14.6
rmApplied times-frac11.0
rmApplied add-cube-cbrt11.3
Applied times-frac6.4
Applied associate-*l*1.1
rmApplied add-cube-cbrt1.2
Final simplification1.2
herbie shell --seed 2020057 +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))))