x + \frac{e^{y \cdot \log \left(\frac{y}{z + y}\right)}}{y}x + \frac{e^{y \cdot \left(2 \cdot \log \left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right) + \log \left(\frac{\sqrt[3]{y}}{\sqrt[3]{z + y}}\right)\right)}}{y}double f(double x, double y, double z) {
double r364226 = x;
double r364227 = y;
double r364228 = z;
double r364229 = r364228 + r364227;
double r364230 = r364227 / r364229;
double r364231 = log(r364230);
double r364232 = r364227 * r364231;
double r364233 = exp(r364232);
double r364234 = r364233 / r364227;
double r364235 = r364226 + r364234;
return r364235;
}
double f(double x, double y, double z) {
double r364236 = x;
double r364237 = y;
double r364238 = 2.0;
double r364239 = cbrt(r364237);
double r364240 = z;
double r364241 = r364240 + r364237;
double r364242 = cbrt(r364241);
double r364243 = r364239 / r364242;
double r364244 = log(r364243);
double r364245 = r364238 * r364244;
double r364246 = r364245 + r364244;
double r364247 = r364237 * r364246;
double r364248 = exp(r364247);
double r364249 = r364248 / r364237;
double r364250 = r364236 + r364249;
return r364250;
}




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