\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) + \left(x \cdot \log \left(\sqrt[3]{\sqrt[3]{y} \cdot \sqrt[3]{y}}\right) + x \cdot \log \left(\sqrt[3]{\sqrt[3]{y}}\right)\right)\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 r502522 = x;
double r502523 = y;
double r502524 = log(r502523);
double r502525 = r502522 * r502524;
double r502526 = z;
double r502527 = 1.0;
double r502528 = r502527 - r502523;
double r502529 = log(r502528);
double r502530 = r502526 * r502529;
double r502531 = r502525 + r502530;
double r502532 = t;
double r502533 = r502531 - r502532;
return r502533;
}
double f(double x, double y, double z, double t) {
double r502534 = x;
double r502535 = 2.0;
double r502536 = y;
double r502537 = cbrt(r502536);
double r502538 = log(r502537);
double r502539 = r502535 * r502538;
double r502540 = r502534 * r502539;
double r502541 = r502537 * r502537;
double r502542 = cbrt(r502541);
double r502543 = log(r502542);
double r502544 = r502534 * r502543;
double r502545 = cbrt(r502537);
double r502546 = log(r502545);
double r502547 = r502534 * r502546;
double r502548 = r502544 + r502547;
double r502549 = r502540 + r502548;
double r502550 = z;
double r502551 = 1.0;
double r502552 = log(r502551);
double r502553 = r502551 * r502536;
double r502554 = 0.5;
double r502555 = pow(r502536, r502535);
double r502556 = pow(r502551, r502535);
double r502557 = r502555 / r502556;
double r502558 = r502554 * r502557;
double r502559 = r502553 + r502558;
double r502560 = r502552 - r502559;
double r502561 = r502550 * r502560;
double r502562 = r502549 + r502561;
double r502563 = t;
double r502564 = r502562 - r502563;
return r502564;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 9.4 |
|---|---|
| Target | 0.3 |
| Herbie | 0.4 |
Initial program 9.4
Taylor expanded around 0 0.3
rmApplied add-cube-cbrt0.3
Applied log-prod0.4
Applied distribute-lft-in0.4
Simplified0.4
rmApplied add-cube-cbrt0.4
Applied cbrt-prod0.4
Applied log-prod0.4
Applied distribute-lft-in0.4
Final simplification0.4
herbie shell --seed 2020001
(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))