x + y \cdot \frac{z - t}{a - t}x + \left(y \cdot \frac{\sqrt[3]{z - t} \cdot \sqrt[3]{z - t}}{\sqrt[3]{a - t} \cdot \sqrt[3]{a - t}}\right) \cdot \frac{\sqrt[3]{z - t}}{\sqrt[3]{a - t}}double f(double x, double y, double z, double t, double a) {
double r559917 = x;
double r559918 = y;
double r559919 = z;
double r559920 = t;
double r559921 = r559919 - r559920;
double r559922 = a;
double r559923 = r559922 - r559920;
double r559924 = r559921 / r559923;
double r559925 = r559918 * r559924;
double r559926 = r559917 + r559925;
return r559926;
}
double f(double x, double y, double z, double t, double a) {
double r559927 = x;
double r559928 = y;
double r559929 = z;
double r559930 = t;
double r559931 = r559929 - r559930;
double r559932 = cbrt(r559931);
double r559933 = r559932 * r559932;
double r559934 = a;
double r559935 = r559934 - r559930;
double r559936 = cbrt(r559935);
double r559937 = r559936 * r559936;
double r559938 = r559933 / r559937;
double r559939 = r559928 * r559938;
double r559940 = r559932 / r559936;
double r559941 = r559939 * r559940;
double r559942 = r559927 + r559941;
return r559942;
}




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 | 0.5 |
| Herbie | 0.6 |
Initial program 1.5
rmApplied add-cube-cbrt2.0
Applied add-cube-cbrt1.9
Applied times-frac1.9
Applied associate-*r*0.6
Final simplification0.6
herbie shell --seed 2019356
(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.508084860551241e-17) (+ x (* y (/ (- z t) (- a t)))) (if (< y 2.894426862792089e-49) (+ x (* (* y (- z t)) (/ 1 (- a t)))) (+ x (* y (/ (- z t) (- a t))))))
(+ x (* y (/ (- z t) (- a t)))))