\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({y}^{\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 r527451 = x;
double r527452 = y;
double r527453 = log(r527452);
double r527454 = r527451 * r527453;
double r527455 = z;
double r527456 = 1.0;
double r527457 = r527456 - r527452;
double r527458 = log(r527457);
double r527459 = r527455 * r527458;
double r527460 = r527454 + r527459;
double r527461 = t;
double r527462 = r527460 - r527461;
return r527462;
}
double f(double x, double y, double z, double t) {
double r527463 = x;
double r527464 = 2.0;
double r527465 = y;
double r527466 = cbrt(r527465);
double r527467 = log(r527466);
double r527468 = r527464 * r527467;
double r527469 = r527463 * r527468;
double r527470 = 0.3333333333333333;
double r527471 = pow(r527465, r527470);
double r527472 = log(r527471);
double r527473 = r527463 * r527472;
double r527474 = r527469 + r527473;
double r527475 = z;
double r527476 = 1.0;
double r527477 = log(r527476);
double r527478 = r527475 * r527477;
double r527479 = r527475 * r527465;
double r527480 = r527476 * r527479;
double r527481 = 0.5;
double r527482 = pow(r527465, r527464);
double r527483 = r527475 * r527482;
double r527484 = pow(r527476, r527464);
double r527485 = r527483 / r527484;
double r527486 = r527481 * r527485;
double r527487 = r527480 + r527486;
double r527488 = r527478 - r527487;
double r527489 = r527474 + r527488;
double r527490 = t;
double r527491 = r527489 - r527490;
return r527491;
}




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.4
rmApplied add-cube-cbrt0.4
Applied log-prod0.4
Applied distribute-lft-in0.4
Simplified0.4
rmApplied pow1/30.4
Final simplification0.4
herbie shell --seed 2020060
(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))