\left(x \cdot \log y + z \cdot \log \left(1 - y\right)\right) - t
\left(x \cdot \log \left(\sqrt[3]{y} \cdot \sqrt[3]{y}\right) + \left(x \cdot \log \left(\sqrt[3]{y}\right) + z \cdot \left(\log 1 - \left(1 \cdot y + \frac{1}{2} \cdot \frac{{y}^{2}}{{1}^{2}}\right)\right)\right)\right) - tdouble f(double x, double y, double z, double t) {
double r526460 = x;
double r526461 = y;
double r526462 = log(r526461);
double r526463 = r526460 * r526462;
double r526464 = z;
double r526465 = 1.0;
double r526466 = r526465 - r526461;
double r526467 = log(r526466);
double r526468 = r526464 * r526467;
double r526469 = r526463 + r526468;
double r526470 = t;
double r526471 = r526469 - r526470;
return r526471;
}
double f(double x, double y, double z, double t) {
double r526472 = x;
double r526473 = y;
double r526474 = cbrt(r526473);
double r526475 = r526474 * r526474;
double r526476 = log(r526475);
double r526477 = r526472 * r526476;
double r526478 = log(r526474);
double r526479 = r526472 * r526478;
double r526480 = z;
double r526481 = 1.0;
double r526482 = log(r526481);
double r526483 = r526481 * r526473;
double r526484 = 0.5;
double r526485 = 2.0;
double r526486 = pow(r526473, r526485);
double r526487 = pow(r526481, r526485);
double r526488 = r526486 / r526487;
double r526489 = r526484 * r526488;
double r526490 = r526483 + r526489;
double r526491 = r526482 - r526490;
double r526492 = r526480 * r526491;
double r526493 = r526479 + r526492;
double r526494 = r526477 + r526493;
double r526495 = t;
double r526496 = r526494 - r526495;
return r526496;
}




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.5 |
Initial program 9.1
Taylor expanded around 0 0.4
rmApplied add-cube-cbrt0.4
Applied log-prod0.5
Applied distribute-lft-in0.5
Applied associate-+l+0.5
Final simplification0.5
herbie shell --seed 2020046
(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))