\left(\left(\left(x - 0.5\right) \cdot \log x - x\right) + 0.9189385332046700050057097541866824030876\right) + \frac{\left(\left(y + 7.936500793651000149400709382518925849581 \cdot 10^{-4}\right) \cdot z - 0.002777777777777800001512975569539776188321\right) \cdot z + 0.08333333333333299564049667651488562114537}{x}\left(\left(x - 0.5\right) \cdot \log \left(\sqrt[3]{x} \cdot {\left(\frac{1}{x}\right)}^{\frac{-1}{3}}\right) + \left(\left(\log \left(\sqrt[3]{x}\right) \cdot \left(x - 0.5\right) - x\right) + 0.9189385332046700050057097541866824030876\right)\right) + \frac{\left(\left(y + 7.936500793651000149400709382518925849581 \cdot 10^{-4}\right) \cdot z - 0.002777777777777800001512975569539776188321\right) \cdot z + 0.08333333333333299564049667651488562114537}{x}double f(double x, double y, double z) {
double r424355 = x;
double r424356 = 0.5;
double r424357 = r424355 - r424356;
double r424358 = log(r424355);
double r424359 = r424357 * r424358;
double r424360 = r424359 - r424355;
double r424361 = 0.91893853320467;
double r424362 = r424360 + r424361;
double r424363 = y;
double r424364 = 0.0007936500793651;
double r424365 = r424363 + r424364;
double r424366 = z;
double r424367 = r424365 * r424366;
double r424368 = 0.0027777777777778;
double r424369 = r424367 - r424368;
double r424370 = r424369 * r424366;
double r424371 = 0.083333333333333;
double r424372 = r424370 + r424371;
double r424373 = r424372 / r424355;
double r424374 = r424362 + r424373;
return r424374;
}
double f(double x, double y, double z) {
double r424375 = x;
double r424376 = 0.5;
double r424377 = r424375 - r424376;
double r424378 = cbrt(r424375);
double r424379 = 1.0;
double r424380 = r424379 / r424375;
double r424381 = -0.3333333333333333;
double r424382 = pow(r424380, r424381);
double r424383 = r424378 * r424382;
double r424384 = log(r424383);
double r424385 = r424377 * r424384;
double r424386 = log(r424378);
double r424387 = r424386 * r424377;
double r424388 = r424387 - r424375;
double r424389 = 0.91893853320467;
double r424390 = r424388 + r424389;
double r424391 = r424385 + r424390;
double r424392 = y;
double r424393 = 0.0007936500793651;
double r424394 = r424392 + r424393;
double r424395 = z;
double r424396 = r424394 * r424395;
double r424397 = 0.0027777777777778;
double r424398 = r424396 - r424397;
double r424399 = r424398 * r424395;
double r424400 = 0.083333333333333;
double r424401 = r424399 + r424400;
double r424402 = r424401 / r424375;
double r424403 = r424391 + r424402;
return r424403;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 6.1 |
|---|---|
| Target | 1.2 |
| Herbie | 6.1 |
Initial program 6.1
rmApplied add-cube-cbrt6.1
Applied log-prod6.1
Applied distribute-lft-in6.1
Applied associate--l+6.1
Applied associate-+l+6.1
Simplified6.1
Taylor expanded around inf 6.1
Final simplification6.1
herbie shell --seed 2019322
(FPCore (x y z)
:name "Numeric.SpecFunctions:$slogFactorial from math-functions-0.1.5.2, B"
:precision binary64
:herbie-target
(+ (+ (+ (* (- x 0.5) (log x)) (- 0.91893853320467001 x)) (/ 0.0833333333333329956 x)) (* (/ z x) (- (* z (+ y 7.93650079365100015e-4)) 0.0027777777777778)))
(+ (+ (- (* (- x 0.5) (log x)) x) 0.91893853320467001) (/ (+ (* (- (* (+ y 7.93650079365100015e-4) z) 0.0027777777777778) z) 0.0833333333333329956) x)))