\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 {y}^{\frac{1}{3}}}\right) + x \cdot \log \left(\sqrt[3]{\sqrt[3]{y}}\right)\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 r359584 = x;
double r359585 = y;
double r359586 = log(r359585);
double r359587 = r359584 * r359586;
double r359588 = z;
double r359589 = 1.0;
double r359590 = r359589 - r359585;
double r359591 = log(r359590);
double r359592 = r359588 * r359591;
double r359593 = r359587 + r359592;
double r359594 = t;
double r359595 = r359593 - r359594;
return r359595;
}
double f(double x, double y, double z, double t) {
double r359596 = x;
double r359597 = 2.0;
double r359598 = y;
double r359599 = cbrt(r359598);
double r359600 = log(r359599);
double r359601 = r359597 * r359600;
double r359602 = r359596 * r359601;
double r359603 = 0.3333333333333333;
double r359604 = pow(r359598, r359603);
double r359605 = r359599 * r359604;
double r359606 = cbrt(r359605);
double r359607 = log(r359606);
double r359608 = r359596 * r359607;
double r359609 = cbrt(r359599);
double r359610 = log(r359609);
double r359611 = r359596 * r359610;
double r359612 = r359608 + r359611;
double r359613 = r359602 + r359612;
double r359614 = z;
double r359615 = 1.0;
double r359616 = log(r359615);
double r359617 = r359614 * r359616;
double r359618 = r359614 * r359598;
double r359619 = r359615 * r359618;
double r359620 = 0.5;
double r359621 = pow(r359598, r359597);
double r359622 = r359614 * r359621;
double r359623 = pow(r359615, r359597);
double r359624 = r359622 / r359623;
double r359625 = r359620 * r359624;
double r359626 = r359619 + r359625;
double r359627 = r359617 - r359626;
double r359628 = r359613 + r359627;
double r359629 = t;
double r359630 = r359628 - r359629;
return r359630;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 9.3 |
|---|---|
| Target | 0.3 |
| Herbie | 0.4 |
Initial program 9.3
Taylor expanded around 0 0.4
rmApplied add-cube-cbrt0.4
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
rmApplied pow1/30.4
Final simplification0.4
herbie shell --seed 2020003
(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))