\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) + x \cdot \log \left({\left(\frac{1}{y}\right)}^{\frac{-1}{3}}\right)\right) + \left(z \cdot \log 1 - \left(1 \cdot \left(z \cdot y\right) + \frac{1}{2} \cdot \frac{z \cdot {y}^{2}}{{1}^{2}}\right)\right)\right) - tdouble f(double x, double y, double z, double t) {
double r380717 = x;
double r380718 = y;
double r380719 = log(r380718);
double r380720 = r380717 * r380719;
double r380721 = z;
double r380722 = 1.0;
double r380723 = r380722 - r380718;
double r380724 = log(r380723);
double r380725 = r380721 * r380724;
double r380726 = r380720 + r380725;
double r380727 = t;
double r380728 = r380726 - r380727;
return r380728;
}
double f(double x, double y, double z, double t) {
double r380729 = x;
double r380730 = 2.0;
double r380731 = y;
double r380732 = cbrt(r380731);
double r380733 = log(r380732);
double r380734 = r380730 * r380733;
double r380735 = r380729 * r380734;
double r380736 = 1.0;
double r380737 = r380736 / r380731;
double r380738 = -0.3333333333333333;
double r380739 = pow(r380737, r380738);
double r380740 = log(r380739);
double r380741 = r380729 * r380740;
double r380742 = r380735 + r380741;
double r380743 = z;
double r380744 = 1.0;
double r380745 = log(r380744);
double r380746 = r380743 * r380745;
double r380747 = r380743 * r380731;
double r380748 = r380744 * r380747;
double r380749 = 0.5;
double r380750 = pow(r380731, r380730);
double r380751 = r380743 * r380750;
double r380752 = pow(r380744, r380730);
double r380753 = r380751 / r380752;
double r380754 = r380749 * r380753;
double r380755 = r380748 + r380754;
double r380756 = r380746 - r380755;
double r380757 = r380742 + r380756;
double r380758 = t;
double r380759 = r380757 - r380758;
return r380759;
}




Bits error versus x




Bits error versus y




Bits error versus z




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