x + \frac{\left(y - z\right) \cdot t}{a - z}x + \left(\sqrt[3]{\frac{y - z}{\sqrt[3]{a - z} \cdot \sqrt[3]{a - z}}} \cdot \sqrt[3]{\frac{y - z}{\sqrt[3]{a - z} \cdot \sqrt[3]{a - z}}}\right) \cdot \left(\left(\sqrt[3]{\frac{y - z}{\sqrt[3]{a - z} \cdot \sqrt[3]{a - z}}} \cdot \frac{\sqrt[3]{t} \cdot \sqrt[3]{t}}{\sqrt[3]{\sqrt[3]{a - z} \cdot \sqrt[3]{a - z}}}\right) \cdot \frac{\sqrt[3]{t}}{\sqrt[3]{\sqrt[3]{a - z}}}\right)double f(double x, double y, double z, double t, double a) {
double r719224 = x;
double r719225 = y;
double r719226 = z;
double r719227 = r719225 - r719226;
double r719228 = t;
double r719229 = r719227 * r719228;
double r719230 = a;
double r719231 = r719230 - r719226;
double r719232 = r719229 / r719231;
double r719233 = r719224 + r719232;
return r719233;
}
double f(double x, double y, double z, double t, double a) {
double r719234 = x;
double r719235 = y;
double r719236 = z;
double r719237 = r719235 - r719236;
double r719238 = a;
double r719239 = r719238 - r719236;
double r719240 = cbrt(r719239);
double r719241 = r719240 * r719240;
double r719242 = r719237 / r719241;
double r719243 = cbrt(r719242);
double r719244 = r719243 * r719243;
double r719245 = t;
double r719246 = cbrt(r719245);
double r719247 = r719246 * r719246;
double r719248 = cbrt(r719241);
double r719249 = r719247 / r719248;
double r719250 = r719243 * r719249;
double r719251 = cbrt(r719240);
double r719252 = r719246 / r719251;
double r719253 = r719250 * r719252;
double r719254 = r719244 * r719253;
double r719255 = r719234 + r719254;
return r719255;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 10.8 |
|---|---|
| Target | 0.6 |
| Herbie | 1.9 |
Initial program 10.8
rmApplied add-cube-cbrt11.2
Applied times-frac2.0
rmApplied add-cube-cbrt2.1
Applied associate-*l*2.1
rmApplied add-cube-cbrt2.1
Applied cbrt-prod2.1
Applied add-cube-cbrt2.2
Applied times-frac2.2
Applied associate-*r*1.9
Final simplification1.9
herbie shell --seed 2020057
(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))))