x + \frac{e^{y \cdot \log \left(\frac{y}{z + y}\right)}}{y}\frac{e^{y \cdot \log \left(\frac{\sqrt[3]{y}}{\sqrt[3]{y + z}}\right) + y \cdot \left(\log \left(\frac{\sqrt[3]{y}}{\sqrt[3]{y + z}}\right) \cdot 2\right)}}{y} + xdouble f(double x, double y, double z) {
double r325573 = x;
double r325574 = y;
double r325575 = z;
double r325576 = r325575 + r325574;
double r325577 = r325574 / r325576;
double r325578 = log(r325577);
double r325579 = r325574 * r325578;
double r325580 = exp(r325579);
double r325581 = r325580 / r325574;
double r325582 = r325573 + r325581;
return r325582;
}
double f(double x, double y, double z) {
double r325583 = y;
double r325584 = cbrt(r325583);
double r325585 = z;
double r325586 = r325583 + r325585;
double r325587 = cbrt(r325586);
double r325588 = r325584 / r325587;
double r325589 = log(r325588);
double r325590 = r325583 * r325589;
double r325591 = 2.0;
double r325592 = r325589 * r325591;
double r325593 = r325583 * r325592;
double r325594 = r325590 + r325593;
double r325595 = exp(r325594);
double r325596 = r325595 / r325583;
double r325597 = x;
double r325598 = r325596 + r325597;
return r325598;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 5.9 |
|---|---|
| Target | 1.1 |
| Herbie | 1.1 |
Initial program 5.9
rmApplied add-cube-cbrt19.1
Applied add-cube-cbrt5.9
Applied times-frac5.9
Applied log-prod2.3
Applied distribute-lft-in2.3
Simplified1.1
Final simplification1.1
herbie shell --seed 2019194 +o rules:numerics
(FPCore (x y z)
:name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2, G"
:herbie-target
(if (< (/ y (+ z y)) 7.1154157597908e-315) (+ x (/ (exp (/ -1.0 z)) y)) (+ x (/ (exp (log (pow (/ y (+ y z)) y))) y)))
(+ x (/ (exp (* y (log (/ y (+ z y))))) y)))