x + \left(y - z\right) \cdot \frac{t - x}{a - z}\begin{array}{l}
\mathbf{if}\;\frac{t - x}{a - z} \cdot \left(y - z\right) + x \le -5.373786410763023350407955876767747404474 \cdot 10^{-308}:\\
\;\;\;\;\frac{t - x}{\sqrt[3]{a - z}} \cdot \frac{y - z}{\sqrt[3]{a - z} \cdot \sqrt[3]{a - z}} + x\\
\mathbf{elif}\;\frac{t - x}{a - z} \cdot \left(y - z\right) + x \le 0.0:\\
\;\;\;\;t + y \cdot \left(\frac{x}{z} - \frac{t}{z}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{t - x}{a - z} \cdot \left(y - z\right) + x\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r6765980 = x;
double r6765981 = y;
double r6765982 = z;
double r6765983 = r6765981 - r6765982;
double r6765984 = t;
double r6765985 = r6765984 - r6765980;
double r6765986 = a;
double r6765987 = r6765986 - r6765982;
double r6765988 = r6765985 / r6765987;
double r6765989 = r6765983 * r6765988;
double r6765990 = r6765980 + r6765989;
return r6765990;
}
double f(double x, double y, double z, double t, double a) {
double r6765991 = t;
double r6765992 = x;
double r6765993 = r6765991 - r6765992;
double r6765994 = a;
double r6765995 = z;
double r6765996 = r6765994 - r6765995;
double r6765997 = r6765993 / r6765996;
double r6765998 = y;
double r6765999 = r6765998 - r6765995;
double r6766000 = r6765997 * r6765999;
double r6766001 = r6766000 + r6765992;
double r6766002 = -5.373786410763023e-308;
bool r6766003 = r6766001 <= r6766002;
double r6766004 = cbrt(r6765996);
double r6766005 = r6765993 / r6766004;
double r6766006 = r6766004 * r6766004;
double r6766007 = r6765999 / r6766006;
double r6766008 = r6766005 * r6766007;
double r6766009 = r6766008 + r6765992;
double r6766010 = 0.0;
bool r6766011 = r6766001 <= r6766010;
double r6766012 = r6765992 / r6765995;
double r6766013 = r6765991 / r6765995;
double r6766014 = r6766012 - r6766013;
double r6766015 = r6765998 * r6766014;
double r6766016 = r6765991 + r6766015;
double r6766017 = r6766011 ? r6766016 : r6766001;
double r6766018 = r6766003 ? r6766009 : r6766017;
return r6766018;
}



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)))) < -5.373786410763023e-308Initial program 7.5
rmApplied add-cube-cbrt8.2
Applied *-un-lft-identity8.2
Applied times-frac8.2
Applied associate-*r*5.3
Simplified5.2
if -5.373786410763023e-308 < (+ x (* (- y z) (/ (- t x) (- a z)))) < 0.0Initial program 61.8
rmApplied add-cube-cbrt61.6
Applied add-cube-cbrt61.6
Applied times-frac61.6
Applied associate-*r*61.4
Simplified61.4
Taylor expanded around inf 26.6
Simplified20.7
if 0.0 < (+ x (* (- y z) (/ (- t x) (- a z)))) Initial program 7.6
Final simplification8.3
herbie shell --seed 2019200
(FPCore (x y z t a)
:name "Numeric.Signal:interpolate from hsignal-0.2.7.1"
(+ x (* (- y z) (/ (- t x) (- a z)))))