x + \left(y - z\right) \cdot \frac{t - x}{a - z}\begin{array}{l}
\mathbf{if}\;x + \frac{t - x}{a - z} \cdot \left(y - z\right) \le -1.923994883449214 \cdot 10^{-303}:\\
\;\;\;\;x + \frac{t - x}{a - z} \cdot \left(y - z\right)\\
\mathbf{elif}\;x + \frac{t - x}{a - z} \cdot \left(y - z\right) \le 0.0:\\
\;\;\;\;\left(t - \frac{t}{\frac{z}{y}}\right) + \frac{x}{\frac{z}{y}}\\
\mathbf{else}:\\
\;\;\;\;\sqrt[3]{\frac{t - x}{\sqrt[3]{a - z}}} \cdot \left(\left(\sqrt[3]{\frac{t - x}{\sqrt[3]{a - z}}} \cdot \sqrt[3]{\frac{t - x}{\sqrt[3]{a - z}}}\right) \cdot \frac{y - z}{\sqrt[3]{a - z} \cdot \sqrt[3]{a - z}}\right) + x\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r6752983 = x;
double r6752984 = y;
double r6752985 = z;
double r6752986 = r6752984 - r6752985;
double r6752987 = t;
double r6752988 = r6752987 - r6752983;
double r6752989 = a;
double r6752990 = r6752989 - r6752985;
double r6752991 = r6752988 / r6752990;
double r6752992 = r6752986 * r6752991;
double r6752993 = r6752983 + r6752992;
return r6752993;
}
double f(double x, double y, double z, double t, double a) {
double r6752994 = x;
double r6752995 = t;
double r6752996 = r6752995 - r6752994;
double r6752997 = a;
double r6752998 = z;
double r6752999 = r6752997 - r6752998;
double r6753000 = r6752996 / r6752999;
double r6753001 = y;
double r6753002 = r6753001 - r6752998;
double r6753003 = r6753000 * r6753002;
double r6753004 = r6752994 + r6753003;
double r6753005 = -1.923994883449214e-303;
bool r6753006 = r6753004 <= r6753005;
double r6753007 = 0.0;
bool r6753008 = r6753004 <= r6753007;
double r6753009 = r6752998 / r6753001;
double r6753010 = r6752995 / r6753009;
double r6753011 = r6752995 - r6753010;
double r6753012 = r6752994 / r6753009;
double r6753013 = r6753011 + r6753012;
double r6753014 = cbrt(r6752999);
double r6753015 = r6752996 / r6753014;
double r6753016 = cbrt(r6753015);
double r6753017 = r6753016 * r6753016;
double r6753018 = r6753014 * r6753014;
double r6753019 = r6753002 / r6753018;
double r6753020 = r6753017 * r6753019;
double r6753021 = r6753016 * r6753020;
double r6753022 = r6753021 + r6752994;
double r6753023 = r6753008 ? r6753013 : r6753022;
double r6753024 = r6753006 ? r6753004 : r6753023;
return r6753024;
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus t



Bits error versus a
Results
if (+ x (* (- y z) (/ (- t x) (- a z)))) < -1.923994883449214e-303Initial program 6.7
if -1.923994883449214e-303 < (+ x (* (- y z) (/ (- t x) (- a z)))) < 0.0Initial program 61.7
rmApplied add-cube-cbrt61.4
Applied *-un-lft-identity61.4
Applied times-frac61.3
Applied associate-*r*61.2
Simplified61.3
rmApplied add-cube-cbrt61.2
Applied associate-*r*61.2
Taylor expanded around inf 25.5
Simplified19.6
if 0.0 < (+ x (* (- y z) (/ (- t x) (- a z)))) Initial program 7.7
rmApplied add-cube-cbrt8.4
Applied *-un-lft-identity8.4
Applied times-frac8.4
Applied associate-*r*5.5
Simplified5.4
rmApplied add-cube-cbrt5.7
Applied associate-*r*5.7
Final simplification7.9
herbie shell --seed 2019164
(FPCore (x y z t a)
:name "Numeric.Signal:interpolate from hsignal-0.2.7.1"
(+ x (* (- y z) (/ (- t x) (- a z)))))