\left(x \cdot \log y + z \cdot \log \left(1.0 - y\right)\right) - t
\left(z \cdot \left(\left(\log 1.0 - 1.0 \cdot y\right) - \frac{\frac{\frac{1}{2}}{\frac{1.0}{y}}}{\frac{1.0}{y}}\right) + \left(\left(\log \left(\sqrt[3]{y}\right) + \log \left(\sqrt[3]{y}\right)\right) \cdot x + x \cdot \log \left({y}^{\frac{1}{3}}\right)\right)\right) - tdouble f(double x, double y, double z, double t) {
double r24098597 = x;
double r24098598 = y;
double r24098599 = log(r24098598);
double r24098600 = r24098597 * r24098599;
double r24098601 = z;
double r24098602 = 1.0;
double r24098603 = r24098602 - r24098598;
double r24098604 = log(r24098603);
double r24098605 = r24098601 * r24098604;
double r24098606 = r24098600 + r24098605;
double r24098607 = t;
double r24098608 = r24098606 - r24098607;
return r24098608;
}
double f(double x, double y, double z, double t) {
double r24098609 = z;
double r24098610 = 1.0;
double r24098611 = log(r24098610);
double r24098612 = y;
double r24098613 = r24098610 * r24098612;
double r24098614 = r24098611 - r24098613;
double r24098615 = 0.5;
double r24098616 = r24098610 / r24098612;
double r24098617 = r24098615 / r24098616;
double r24098618 = r24098617 / r24098616;
double r24098619 = r24098614 - r24098618;
double r24098620 = r24098609 * r24098619;
double r24098621 = cbrt(r24098612);
double r24098622 = log(r24098621);
double r24098623 = r24098622 + r24098622;
double r24098624 = x;
double r24098625 = r24098623 * r24098624;
double r24098626 = 0.3333333333333333;
double r24098627 = pow(r24098612, r24098626);
double r24098628 = log(r24098627);
double r24098629 = r24098624 * r24098628;
double r24098630 = r24098625 + r24098629;
double r24098631 = r24098620 + r24098630;
double r24098632 = t;
double r24098633 = r24098631 - r24098632;
return r24098633;
}




Bits error versus x




Bits error versus y




Bits error versus z




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