x + \left(y - z\right) \cdot \frac{t - x}{a - z}\begin{array}{l}
\mathbf{if}\;x + \left(y - z\right) \cdot \frac{t - x}{a - z} \le -9.30598028412173837 \cdot 10^{-248} \lor \neg \left(x + \left(y - z\right) \cdot \frac{t - x}{a - z} \le 0.0\right):\\
\;\;\;\;x + \frac{y - z}{\sqrt[3]{a - z} \cdot \sqrt[3]{a - z}} \cdot \frac{t - x}{\sqrt[3]{a - z}}\\
\mathbf{else}:\\
\;\;\;\;y \cdot \left(\frac{x}{z} - \frac{t}{z}\right) + t\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r127971 = x;
double r127972 = y;
double r127973 = z;
double r127974 = r127972 - r127973;
double r127975 = t;
double r127976 = r127975 - r127971;
double r127977 = a;
double r127978 = r127977 - r127973;
double r127979 = r127976 / r127978;
double r127980 = r127974 * r127979;
double r127981 = r127971 + r127980;
return r127981;
}
double f(double x, double y, double z, double t, double a) {
double r127982 = x;
double r127983 = y;
double r127984 = z;
double r127985 = r127983 - r127984;
double r127986 = t;
double r127987 = r127986 - r127982;
double r127988 = a;
double r127989 = r127988 - r127984;
double r127990 = r127987 / r127989;
double r127991 = r127985 * r127990;
double r127992 = r127982 + r127991;
double r127993 = -9.305980284121738e-248;
bool r127994 = r127992 <= r127993;
double r127995 = 0.0;
bool r127996 = r127992 <= r127995;
double r127997 = !r127996;
bool r127998 = r127994 || r127997;
double r127999 = cbrt(r127989);
double r128000 = r127999 * r127999;
double r128001 = r127985 / r128000;
double r128002 = r127987 / r127999;
double r128003 = r128001 * r128002;
double r128004 = r127982 + r128003;
double r128005 = r127982 / r127984;
double r128006 = r127986 / r127984;
double r128007 = r128005 - r128006;
double r128008 = r127983 * r128007;
double r128009 = r128008 + r127986;
double r128010 = r127998 ? r128004 : r128009;
return r128010;
}



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)))) < -9.305980284121738e-248 or 0.0 < (+ x (* (- y z) (/ (- t x) (- a z)))) Initial program 7.1
rmApplied add-cube-cbrt7.8
Applied *-un-lft-identity7.8
Applied times-frac7.8
Applied associate-*r*4.9
Simplified4.9
if -9.305980284121738e-248 < (+ x (* (- y z) (/ (- t x) (- a z)))) < 0.0Initial program 58.8
rmApplied add-cube-cbrt58.6
Applied *-un-lft-identity58.6
Applied times-frac58.6
Applied associate-*r*56.5
Simplified56.5
rmApplied add-cube-cbrt56.4
Taylor expanded around inf 27.0
Simplified23.0
Final simplification7.5
herbie shell --seed 2020046
(FPCore (x y z t a)
:name "Numeric.Signal:interpolate from hsignal-0.2.7.1"
:precision binary64
(+ x (* (- y z) (/ (- t x) (- a z)))))