\left(x \cdot \log y + z \cdot \log \left(1 - y\right)\right) - t
\left(\left(x \cdot \log \left(\sqrt[3]{y}\right) + \left(\left(\frac{y}{1} \cdot y\right) \cdot \frac{\frac{-1}{2}}{1} + \left(\log 1 - y \cdot 1\right)\right) \cdot z\right) + \log \left({y}^{\frac{2}{3}}\right) \cdot x\right) - tdouble f(double x, double y, double z, double t) {
double r379271 = x;
double r379272 = y;
double r379273 = log(r379272);
double r379274 = r379271 * r379273;
double r379275 = z;
double r379276 = 1.0;
double r379277 = r379276 - r379272;
double r379278 = log(r379277);
double r379279 = r379275 * r379278;
double r379280 = r379274 + r379279;
double r379281 = t;
double r379282 = r379280 - r379281;
return r379282;
}
double f(double x, double y, double z, double t) {
double r379283 = x;
double r379284 = y;
double r379285 = cbrt(r379284);
double r379286 = log(r379285);
double r379287 = r379283 * r379286;
double r379288 = 1.0;
double r379289 = r379284 / r379288;
double r379290 = r379289 * r379284;
double r379291 = -0.5;
double r379292 = r379291 / r379288;
double r379293 = r379290 * r379292;
double r379294 = log(r379288);
double r379295 = r379284 * r379288;
double r379296 = r379294 - r379295;
double r379297 = r379293 + r379296;
double r379298 = z;
double r379299 = r379297 * r379298;
double r379300 = r379287 + r379299;
double r379301 = 0.6666666666666666;
double r379302 = pow(r379284, r379301);
double r379303 = log(r379302);
double r379304 = r379303 * r379283;
double r379305 = r379300 + r379304;
double r379306 = t;
double r379307 = r379305 - r379306;
return r379307;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 9.6 |
|---|---|
| Target | 0.2 |
| Herbie | 0.4 |
Initial program 9.6
Taylor expanded around 0 0.3
Simplified0.3
rmApplied add-cube-cbrt0.3
Applied log-prod0.4
Applied distribute-rgt-in0.4
Applied associate-+l+0.4
Simplified0.4
rmApplied *-un-lft-identity0.4
Applied cbrt-prod0.4
Applied *-un-lft-identity0.4
Applied cbrt-prod0.4
Applied swap-sqr0.4
Simplified0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2019174
(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) (* (/ 0.3333333333333333 (* 1.0 (* 1.0 1.0))) (* y (* y y))))) (- t (* x (log y))))
(- (+ (* x (log y)) (* z (log (- 1.0 y)))) t))