x + \frac{y \cdot \left(z - t\right)}{z - a}\frac{\sqrt[3]{z - t} \cdot \sqrt[3]{z - t}}{\sqrt[3]{z - a} \cdot \sqrt[3]{z - a}} \cdot \left(\frac{\sqrt[3]{z - t}}{\sqrt[3]{z - a}} \cdot y\right) + xdouble f(double x, double y, double z, double t, double a) {
double r676802 = x;
double r676803 = y;
double r676804 = z;
double r676805 = t;
double r676806 = r676804 - r676805;
double r676807 = r676803 * r676806;
double r676808 = a;
double r676809 = r676804 - r676808;
double r676810 = r676807 / r676809;
double r676811 = r676802 + r676810;
return r676811;
}
double f(double x, double y, double z, double t, double a) {
double r676812 = z;
double r676813 = t;
double r676814 = r676812 - r676813;
double r676815 = cbrt(r676814);
double r676816 = r676815 * r676815;
double r676817 = a;
double r676818 = r676812 - r676817;
double r676819 = cbrt(r676818);
double r676820 = r676819 * r676819;
double r676821 = r676816 / r676820;
double r676822 = r676815 / r676819;
double r676823 = y;
double r676824 = r676822 * r676823;
double r676825 = r676821 * r676824;
double r676826 = x;
double r676827 = r676825 + r676826;
return r676827;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 10.6 |
|---|---|
| Target | 1.2 |
| Herbie | 0.5 |
Initial program 10.6
Simplified2.9
rmApplied clear-num3.1
rmApplied fma-udef3.1
Simplified2.9
rmApplied associate-/r/1.3
rmApplied add-cube-cbrt1.9
Applied add-cube-cbrt1.7
Applied times-frac1.7
Applied associate-*l*0.5
Final simplification0.5
herbie shell --seed 2019354 +o rules:numerics
(FPCore (x y z t a)
:name "Graphics.Rendering.Plot.Render.Plot.Axis:renderAxisTicks from plot-0.2.3.4, A"
:precision binary64
:herbie-target
(+ x (/ y (/ (- z a) (- z t))))
(+ x (/ (* y (- z t)) (- z a))))