x + \frac{\left(y - z\right) \cdot \left(t - x\right)}{a - z}\begin{array}{l}
\mathbf{if}\;a \le -3.652286650956768103984227980253051511603 \cdot 10^{-176}:\\
\;\;\;\;x + \frac{\sqrt[3]{y - z} \cdot \sqrt[3]{y - z}}{\sqrt[3]{a - z}} \cdot \left(\frac{\sqrt[3]{y - z}}{\sqrt[3]{a - z}} \cdot \frac{t - x}{\sqrt[3]{a - z}}\right)\\
\mathbf{elif}\;a \le 7.246661161646684881081730449605960355659 \cdot 10^{-111}:\\
\;\;\;\;\left(\frac{x \cdot y}{z} + t\right) - \frac{t \cdot y}{z}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{\sqrt[3]{y - z} \cdot \left(\sqrt[3]{\sqrt[3]{y - z} \cdot \sqrt[3]{y - z}} \cdot \sqrt[3]{\sqrt[3]{y - z}}\right)}{\sqrt[3]{a - z}} \cdot \left(\frac{\sqrt[3]{y - z}}{\sqrt[3]{a - z}} \cdot \frac{t - x}{\sqrt[3]{a - z}}\right)\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r457867 = x;
double r457868 = y;
double r457869 = z;
double r457870 = r457868 - r457869;
double r457871 = t;
double r457872 = r457871 - r457867;
double r457873 = r457870 * r457872;
double r457874 = a;
double r457875 = r457874 - r457869;
double r457876 = r457873 / r457875;
double r457877 = r457867 + r457876;
return r457877;
}
double f(double x, double y, double z, double t, double a) {
double r457878 = a;
double r457879 = -3.652286650956768e-176;
bool r457880 = r457878 <= r457879;
double r457881 = x;
double r457882 = y;
double r457883 = z;
double r457884 = r457882 - r457883;
double r457885 = cbrt(r457884);
double r457886 = r457885 * r457885;
double r457887 = r457878 - r457883;
double r457888 = cbrt(r457887);
double r457889 = r457886 / r457888;
double r457890 = r457885 / r457888;
double r457891 = t;
double r457892 = r457891 - r457881;
double r457893 = r457892 / r457888;
double r457894 = r457890 * r457893;
double r457895 = r457889 * r457894;
double r457896 = r457881 + r457895;
double r457897 = 7.246661161646685e-111;
bool r457898 = r457878 <= r457897;
double r457899 = r457881 * r457882;
double r457900 = r457899 / r457883;
double r457901 = r457900 + r457891;
double r457902 = r457891 * r457882;
double r457903 = r457902 / r457883;
double r457904 = r457901 - r457903;
double r457905 = cbrt(r457886);
double r457906 = cbrt(r457885);
double r457907 = r457905 * r457906;
double r457908 = r457885 * r457907;
double r457909 = r457908 / r457888;
double r457910 = r457909 * r457894;
double r457911 = r457881 + r457910;
double r457912 = r457898 ? r457904 : r457911;
double r457913 = r457880 ? r457896 : r457912;
return r457913;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 24.3 |
|---|---|
| Target | 12.1 |
| Herbie | 10.8 |
if a < -3.652286650956768e-176Initial program 22.8
rmApplied add-cube-cbrt23.2
Applied times-frac10.3
rmApplied add-cube-cbrt10.2
Applied times-frac10.2
Applied associate-*l*10.1
if -3.652286650956768e-176 < a < 7.246661161646685e-111Initial program 29.6
Taylor expanded around inf 14.2
if 7.246661161646685e-111 < a Initial program 22.4
rmApplied add-cube-cbrt22.8
Applied times-frac9.6
rmApplied add-cube-cbrt9.6
Applied times-frac9.6
Applied associate-*l*9.3
rmApplied add-cube-cbrt9.4
Applied cbrt-prod9.4
Final simplification10.8
herbie shell --seed 2019208
(FPCore (x y z t a)
:name "Graphics.Rendering.Chart.Axis.Types:invLinMap from Chart-1.5.3"
:precision binary64
:herbie-target
(if (< z -1.25361310560950359e188) (- t (* (/ y z) (- t x))) (if (< z 4.44670236911381103e64) (+ x (/ (- y z) (/ (- a z) (- t x)))) (- t (* (/ y z) (- t x)))))
(+ x (/ (* (- y z) (- t x)) (- a z))))