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) + \left(y \cdot \log \left(\frac{\sqrt[3]{y}}{\sqrt[3]{y + z}}\right) + y \cdot \log \left(\frac{\sqrt[3]{y}}{\sqrt[3]{y + z}}\right)\right)}}{y} + xdouble f(double x, double y, double z) {
double r16129285 = x;
double r16129286 = y;
double r16129287 = z;
double r16129288 = r16129287 + r16129286;
double r16129289 = r16129286 / r16129288;
double r16129290 = log(r16129289);
double r16129291 = r16129286 * r16129290;
double r16129292 = exp(r16129291);
double r16129293 = r16129292 / r16129286;
double r16129294 = r16129285 + r16129293;
return r16129294;
}
double f(double x, double y, double z) {
double r16129295 = y;
double r16129296 = cbrt(r16129295);
double r16129297 = z;
double r16129298 = r16129295 + r16129297;
double r16129299 = cbrt(r16129298);
double r16129300 = r16129296 / r16129299;
double r16129301 = log(r16129300);
double r16129302 = r16129295 * r16129301;
double r16129303 = r16129302 + r16129302;
double r16129304 = r16129302 + r16129303;
double r16129305 = exp(r16129304);
double r16129306 = r16129305 / r16129295;
double r16129307 = x;
double r16129308 = r16129306 + r16129307;
return r16129308;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 5.8 |
|---|---|
| Target | 1.0 |
| Herbie | 1.0 |
Initial program 5.8
rmApplied add-cube-cbrt19.6
Applied add-cube-cbrt5.8
Applied times-frac5.8
Applied log-prod2.0
Applied distribute-lft-in2.0
Simplified1.0
Final simplification1.0
herbie shell --seed 2019168 +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 z)) y)) (+ x (/ (exp (log (pow (/ y (+ y z)) y))) y)))
(+ x (/ (exp (* y (log (/ y (+ z y))))) y)))