x + y \cdot \frac{z - t}{z - a}\frac{\sqrt[3]{z - t} \cdot \sqrt[3]{z - t}}{\sqrt[3]{z - a} \cdot \sqrt[3]{z - a}} \cdot \frac{y}{\frac{\sqrt[3]{z - a}}{\sqrt[3]{z - t}}} + xdouble f(double x, double y, double z, double t, double a) {
double r663334 = x;
double r663335 = y;
double r663336 = z;
double r663337 = t;
double r663338 = r663336 - r663337;
double r663339 = a;
double r663340 = r663336 - r663339;
double r663341 = r663338 / r663340;
double r663342 = r663335 * r663341;
double r663343 = r663334 + r663342;
return r663343;
}
double f(double x, double y, double z, double t, double a) {
double r663344 = z;
double r663345 = t;
double r663346 = r663344 - r663345;
double r663347 = cbrt(r663346);
double r663348 = r663347 * r663347;
double r663349 = a;
double r663350 = r663344 - r663349;
double r663351 = cbrt(r663350);
double r663352 = r663351 * r663351;
double r663353 = r663348 / r663352;
double r663354 = y;
double r663355 = r663351 / r663347;
double r663356 = r663354 / r663355;
double r663357 = r663353 * r663356;
double r663358 = x;
double r663359 = r663357 + r663358;
return r663359;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 1.5 |
|---|---|
| Target | 1.4 |
| Herbie | 0.6 |
Initial program 1.5
Simplified1.5
rmApplied clear-num1.6
rmApplied fma-udef1.6
Simplified1.4
rmApplied add-cube-cbrt1.9
Applied add-cube-cbrt1.8
Applied times-frac1.8
Applied *-un-lft-identity1.8
Applied times-frac0.6
Simplified0.6
Final simplification0.6
herbie shell --seed 2020060 +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)))))