x + \frac{y \cdot \left(z - t\right)}{a - t}\frac{\sqrt[3]{z - t} \cdot \sqrt[3]{z - t}}{\sqrt[3]{a - t} \cdot \sqrt[3]{a - t}} \cdot \left(\frac{\sqrt[3]{z - t}}{\sqrt[3]{a - t}} \cdot y\right) + xdouble f(double x, double y, double z, double t, double a) {
double r681562 = x;
double r681563 = y;
double r681564 = z;
double r681565 = t;
double r681566 = r681564 - r681565;
double r681567 = r681563 * r681566;
double r681568 = a;
double r681569 = r681568 - r681565;
double r681570 = r681567 / r681569;
double r681571 = r681562 + r681570;
return r681571;
}
double f(double x, double y, double z, double t, double a) {
double r681572 = z;
double r681573 = t;
double r681574 = r681572 - r681573;
double r681575 = cbrt(r681574);
double r681576 = r681575 * r681575;
double r681577 = a;
double r681578 = r681577 - r681573;
double r681579 = cbrt(r681578);
double r681580 = r681579 * r681579;
double r681581 = r681576 / r681580;
double r681582 = r681575 / r681579;
double r681583 = y;
double r681584 = r681582 * r681583;
double r681585 = r681581 * r681584;
double r681586 = x;
double r681587 = r681585 + r681586;
return r681587;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 11.2 |
|---|---|
| Target | 1.2 |
| Herbie | 0.6 |
Initial program 11.2
Simplified3.2
rmApplied *-un-lft-identity3.2
Applied *-un-lft-identity3.2
Applied times-frac3.2
Simplified3.2
rmApplied clear-num3.4
rmApplied fma-udef3.4
Simplified1.4
rmApplied add-cube-cbrt1.9
Applied add-cube-cbrt1.7
Applied times-frac1.7
Applied associate-*l*0.6
Final simplification0.6
herbie shell --seed 2020020 +o rules:numerics
(FPCore (x y z t a)
:name "Graphics.Rendering.Plot.Render.Plot.Axis:renderAxisTicks from plot-0.2.3.4, B"
:precision binary64
:herbie-target
(+ x (/ y (/ (- a t) (- z t))))
(+ x (/ (* y (- z t)) (- a t))))