x + \frac{\left(y - z\right) \cdot t}{a - z}x + \left(\left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right) \cdot \frac{y - z}{\left(\sqrt[3]{\sqrt[3]{a - z} \cdot \sqrt[3]{a - z}} \cdot \sqrt[3]{\sqrt[3]{a - z}}\right) \cdot \sqrt[3]{a - z}}\right) \cdot \frac{\sqrt[3]{t}}{\sqrt[3]{a - z}}double f(double x, double y, double z, double t, double a) {
double r386043 = x;
double r386044 = y;
double r386045 = z;
double r386046 = r386044 - r386045;
double r386047 = t;
double r386048 = r386046 * r386047;
double r386049 = a;
double r386050 = r386049 - r386045;
double r386051 = r386048 / r386050;
double r386052 = r386043 + r386051;
return r386052;
}
double f(double x, double y, double z, double t, double a) {
double r386053 = x;
double r386054 = t;
double r386055 = cbrt(r386054);
double r386056 = r386055 * r386055;
double r386057 = y;
double r386058 = z;
double r386059 = r386057 - r386058;
double r386060 = a;
double r386061 = r386060 - r386058;
double r386062 = cbrt(r386061);
double r386063 = r386062 * r386062;
double r386064 = cbrt(r386063);
double r386065 = cbrt(r386062);
double r386066 = r386064 * r386065;
double r386067 = r386066 * r386062;
double r386068 = r386059 / r386067;
double r386069 = r386056 * r386068;
double r386070 = r386055 / r386062;
double r386071 = r386069 * r386070;
double r386072 = r386053 + r386071;
return r386072;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 10.9 |
|---|---|
| Target | 0.5 |
| Herbie | 1.4 |
Initial program 10.9
rmApplied add-cube-cbrt11.3
Applied times-frac1.7
rmApplied *-un-lft-identity1.7
Applied cbrt-prod1.7
Applied add-cube-cbrt1.8
Applied times-frac1.8
Applied associate-*r*1.4
Simplified1.4
rmApplied add-cube-cbrt1.4
Applied cbrt-prod1.4
Final simplification1.4
herbie shell --seed 2019304
(FPCore (x y z t a)
:name "Graphics.Rendering.Plot.Render.Plot.Axis:renderAxisTick from plot-0.2.3.4, A"
:precision binary64
:herbie-target
(if (< t -1.0682974490174067e-39) (+ x (* (/ (- y z) (- a z)) t)) (if (< t 3.9110949887586375e-141) (+ x (/ (* (- y z) t) (- a z))) (+ x (* (/ (- y z) (- a z)) t))))
(+ x (/ (* (- y z) t) (- a z))))