x + \frac{e^{y \cdot \log \left(\frac{y}{z + y}\right)}}{y}x + \frac{e^{\left(2 \cdot \log \left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)\right) \cdot y + \log \left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right) \cdot y}}{y}double f(double x, double y, double z) {
double r290554 = x;
double r290555 = y;
double r290556 = z;
double r290557 = r290556 + r290555;
double r290558 = r290555 / r290557;
double r290559 = log(r290558);
double r290560 = r290555 * r290559;
double r290561 = exp(r290560);
double r290562 = r290561 / r290555;
double r290563 = r290554 + r290562;
return r290563;
}
double f(double x, double y, double z) {
double r290564 = x;
double r290565 = 2.0;
double r290566 = y;
double r290567 = cbrt(r290566);
double r290568 = z;
double r290569 = r290568 + r290566;
double r290570 = cbrt(r290569);
double r290571 = r290567 / r290570;
double r290572 = log(r290571);
double r290573 = r290565 * r290572;
double r290574 = r290573 * r290566;
double r290575 = r290572 * r290566;
double r290576 = r290574 + r290575;
double r290577 = exp(r290576);
double r290578 = r290577 / r290566;
double r290579 = r290564 + r290578;
return r290579;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 6.0 |
|---|---|
| Target | 1.1 |
| Herbie | 1.1 |
Initial program 6.0
rmApplied add-cube-cbrt19.4
Applied add-cube-cbrt6.1
Applied times-frac6.1
Applied log-prod2.1
Applied distribute-lft-in2.1
Simplified1.1
Simplified1.1
Final simplification1.1
herbie shell --seed 2019303
(FPCore (x y z)
:name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2, G"
:precision binary64
:herbie-target
(if (< (/ y (+ z y)) 7.1154157598e-315) (+ x (/ (exp (/ -1 z)) y)) (+ x (/ (exp (log (pow (/ y (+ y z)) y))) y)))
(+ x (/ (exp (* y (log (/ y (+ z y))))) y)))