x + y \cdot \frac{z - t}{z - a}x + \left(y \cdot \frac{\sqrt[3]{z - t} \cdot \sqrt[3]{z - t}}{\sqrt[3]{z - a} \cdot \sqrt[3]{z - a}}\right) \cdot \frac{\sqrt[3]{z - t}}{\sqrt[3]{z - a}}double f(double x, double y, double z, double t, double a) {
double r413597 = x;
double r413598 = y;
double r413599 = z;
double r413600 = t;
double r413601 = r413599 - r413600;
double r413602 = a;
double r413603 = r413599 - r413602;
double r413604 = r413601 / r413603;
double r413605 = r413598 * r413604;
double r413606 = r413597 + r413605;
return r413606;
}
double f(double x, double y, double z, double t, double a) {
double r413607 = x;
double r413608 = y;
double r413609 = z;
double r413610 = t;
double r413611 = r413609 - r413610;
double r413612 = cbrt(r413611);
double r413613 = r413612 * r413612;
double r413614 = a;
double r413615 = r413609 - r413614;
double r413616 = cbrt(r413615);
double r413617 = r413616 * r413616;
double r413618 = r413613 / r413617;
double r413619 = r413608 * r413618;
double r413620 = r413612 / r413616;
double r413621 = r413619 * r413620;
double r413622 = r413607 + r413621;
return r413622;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 1.7 |
|---|---|
| Target | 1.5 |
| Herbie | 0.5 |
Initial program 1.7
rmApplied add-cube-cbrt2.2
Applied add-cube-cbrt2.0
Applied times-frac2.0
Applied associate-*r*0.5
Final simplification0.5
herbie shell --seed 2019323 +o rules:numerics
(FPCore (x y z t a)
:name "Graphics.Rendering.Plot.Render.Plot.Axis:renderAxisLine from plot-0.2.3.4, A"
:precision binary64
:herbie-target
(+ x (/ y (/ (- z a) (- z t))))
(+ x (* y (/ (- z t) (- z a)))))