\left(x \cdot \log y + z \cdot \log \left(1 - y\right)\right) - t
\left(\left(\left(2 \cdot \log \left(\sqrt[3]{y}\right)\right) \cdot x + \log \left({\left(\frac{1}{y}\right)}^{\frac{-1}{3}}\right) \cdot x\right) + z \cdot \left(\log 1 - \left(1 \cdot y + \frac{1}{2} \cdot \frac{{y}^{2}}{{1}^{2}}\right)\right)\right) - tdouble f(double x, double y, double z, double t) {
double r307475 = x;
double r307476 = y;
double r307477 = log(r307476);
double r307478 = r307475 * r307477;
double r307479 = z;
double r307480 = 1.0;
double r307481 = r307480 - r307476;
double r307482 = log(r307481);
double r307483 = r307479 * r307482;
double r307484 = r307478 + r307483;
double r307485 = t;
double r307486 = r307484 - r307485;
return r307486;
}
double f(double x, double y, double z, double t) {
double r307487 = 2.0;
double r307488 = y;
double r307489 = cbrt(r307488);
double r307490 = log(r307489);
double r307491 = r307487 * r307490;
double r307492 = x;
double r307493 = r307491 * r307492;
double r307494 = 1.0;
double r307495 = r307494 / r307488;
double r307496 = -0.3333333333333333;
double r307497 = pow(r307495, r307496);
double r307498 = log(r307497);
double r307499 = r307498 * r307492;
double r307500 = r307493 + r307499;
double r307501 = z;
double r307502 = 1.0;
double r307503 = log(r307502);
double r307504 = r307502 * r307488;
double r307505 = 0.5;
double r307506 = pow(r307488, r307487);
double r307507 = pow(r307502, r307487);
double r307508 = r307506 / r307507;
double r307509 = r307505 * r307508;
double r307510 = r307504 + r307509;
double r307511 = r307503 - r307510;
double r307512 = r307501 * r307511;
double r307513 = r307500 + r307512;
double r307514 = t;
double r307515 = r307513 - r307514;
return r307515;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 9.7 |
|---|---|
| Target | 0.3 |
| Herbie | 0.4 |
Initial program 9.7
Taylor expanded around 0 0.3
rmApplied add-cube-cbrt0.3
Applied log-prod0.4
Applied distribute-lft-in0.4
Simplified0.4
Simplified0.4
Taylor expanded around inf 0.4
Final simplification0.4
herbie shell --seed 2019323
(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.3333333333333333 (* 1 (* 1 1))) (* y (* y y))))) (- t (* x (log y))))
(- (+ (* x (log y)) (* z (log (- 1 y)))) t))