x + \frac{e^{y \cdot \log \left(\frac{y}{z + y}\right)}}{y}x + \frac{{\left(e^{y}\right)}^{\left(\log \left(\frac{y}{y + z}\right)\right)}}{y}double f(double x, double y, double z) {
double r18283490 = x;
double r18283491 = y;
double r18283492 = z;
double r18283493 = r18283492 + r18283491;
double r18283494 = r18283491 / r18283493;
double r18283495 = log(r18283494);
double r18283496 = r18283491 * r18283495;
double r18283497 = exp(r18283496);
double r18283498 = r18283497 / r18283491;
double r18283499 = r18283490 + r18283498;
return r18283499;
}
double f(double x, double y, double z) {
double r18283500 = x;
double r18283501 = y;
double r18283502 = exp(r18283501);
double r18283503 = z;
double r18283504 = r18283501 + r18283503;
double r18283505 = r18283501 / r18283504;
double r18283506 = log(r18283505);
double r18283507 = pow(r18283502, r18283506);
double r18283508 = r18283507 / r18283501;
double r18283509 = r18283500 + r18283508;
return r18283509;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 5.8 |
|---|---|
| Target | 1.2 |
| Herbie | 1.3 |
Initial program 5.8
rmApplied add-log-exp35.8
Applied exp-to-pow1.3
Final simplification1.3
herbie shell --seed 2019165 +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)))