\left(x \cdot \log y + z \cdot \log \left(1 - y\right)\right) - t
\left(\left(x \cdot \left(2 \cdot \log \left(\sqrt[3]{y}\right)\right) + \log \left({y}^{\frac{1}{3}}\right) \cdot x\right) + \left(z \cdot \left(\log 1 - 1 \cdot y\right) - \frac{1}{2} \cdot \frac{z \cdot {y}^{2}}{{1}^{2}}\right)\right) - tdouble f(double x, double y, double z, double t) {
double r367486 = x;
double r367487 = y;
double r367488 = log(r367487);
double r367489 = r367486 * r367488;
double r367490 = z;
double r367491 = 1.0;
double r367492 = r367491 - r367487;
double r367493 = log(r367492);
double r367494 = r367490 * r367493;
double r367495 = r367489 + r367494;
double r367496 = t;
double r367497 = r367495 - r367496;
return r367497;
}
double f(double x, double y, double z, double t) {
double r367498 = x;
double r367499 = 2.0;
double r367500 = y;
double r367501 = cbrt(r367500);
double r367502 = log(r367501);
double r367503 = r367499 * r367502;
double r367504 = r367498 * r367503;
double r367505 = 0.3333333333333333;
double r367506 = pow(r367500, r367505);
double r367507 = log(r367506);
double r367508 = r367507 * r367498;
double r367509 = r367504 + r367508;
double r367510 = z;
double r367511 = 1.0;
double r367512 = log(r367511);
double r367513 = r367511 * r367500;
double r367514 = r367512 - r367513;
double r367515 = r367510 * r367514;
double r367516 = 0.5;
double r367517 = pow(r367500, r367499);
double r367518 = r367510 * r367517;
double r367519 = pow(r367511, r367499);
double r367520 = r367518 / r367519;
double r367521 = r367516 * r367520;
double r367522 = r367515 - r367521;
double r367523 = r367509 + r367522;
double r367524 = t;
double r367525 = r367523 - r367524;
return r367525;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 9.9 |
|---|---|
| Target | 0.2 |
| Herbie | 0.3 |
Initial program 9.9
Taylor expanded around 0 0.3
Simplified0.3
rmApplied add-cube-cbrt0.3
Applied log-prod0.4
Applied distribute-lft-in0.4
Simplified0.4
Simplified0.4
rmApplied pow1/30.3
Final simplification0.3
herbie shell --seed 2019303
(FPCore (x y z t)
:name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2, B"
:precision binary64
:herbie-target
(- (* (- z) (+ (+ (* 0.5 (* y y)) y) (* (/ 0.333333333333333315 (* 1 (* 1 1))) (* y (* y y))))) (- t (* x (log y))))
(- (+ (* x (log y)) (* z (log (- 1 y)))) t))