x + y \cdot \frac{z - t}{a - t}\frac{\frac{y}{\sqrt[3]{a - t} \cdot \sqrt[3]{a - t}}}{\frac{\sqrt[3]{a - t}}{z - t}} + xdouble f(double x, double y, double z, double t, double a) {
double r446489 = x;
double r446490 = y;
double r446491 = z;
double r446492 = t;
double r446493 = r446491 - r446492;
double r446494 = a;
double r446495 = r446494 - r446492;
double r446496 = r446493 / r446495;
double r446497 = r446490 * r446496;
double r446498 = r446489 + r446497;
return r446498;
}
double f(double x, double y, double z, double t, double a) {
double r446499 = y;
double r446500 = a;
double r446501 = t;
double r446502 = r446500 - r446501;
double r446503 = cbrt(r446502);
double r446504 = r446503 * r446503;
double r446505 = r446499 / r446504;
double r446506 = z;
double r446507 = r446506 - r446501;
double r446508 = r446503 / r446507;
double r446509 = r446505 / r446508;
double r446510 = x;
double r446511 = r446509 + r446510;
return r446511;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 1.3 |
|---|---|
| Target | 0.4 |
| Herbie | 2.3 |
Initial program 1.3
Simplified1.3
rmApplied clear-num1.4
rmApplied fma-udef1.4
Simplified1.3
rmApplied *-un-lft-identity1.3
Applied add-cube-cbrt1.8
Applied times-frac1.8
Applied associate-/r*2.3
Simplified2.3
Final simplification2.3
herbie shell --seed 2019209 +o rules:numerics
(FPCore (x y z t a)
:name "Graphics.Rendering.Plot.Render.Plot.Axis:renderAxisLine from plot-0.2.3.4, B"
:precision binary64
:herbie-target
(if (< y -8.50808486055124107e-17) (+ x (* y (/ (- z t) (- a t)))) (if (< y 2.8944268627920891e-49) (+ x (* (* y (- z t)) (/ 1 (- a t)))) (+ x (* y (/ (- z t) (- a t))))))
(+ x (* y (/ (- z t) (- a t)))))