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 -1.2753538938243992 \cdot 10^{-306} \lor \neg \left(x + \left(y - z\right) \cdot \frac{t - x}{a - z} \le 0.0\right):\\
\;\;\;\;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{else}:\\
\;\;\;\;\left(\frac{x \cdot y}{z} + t\right) - \frac{t \cdot y}{z}\\
\end{array}double code(double x, double y, double z, double t, double a) {
return (x + ((y - z) * ((t - x) / (a - z))));
}
double code(double x, double y, double z, double t, double a) {
double temp;
if ((((x + ((y - z) * ((t - x) / (a - z)))) <= -1.2753538938243992e-306) || !((x + ((y - z) * ((t - x) / (a - z)))) <= 0.0))) {
temp = (x + (((cbrt((y - z)) * cbrt((y - z))) / cbrt((a - z))) * ((cbrt((y - z)) / cbrt((a - z))) * ((t - x) / cbrt((a - z))))));
} else {
temp = ((((x * y) / z) + t) - ((t * y) / z));
}
return temp;
}



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.2753538938243992e-306 or 0.0 < (+ x (* (- y z) (/ (- t x) (- a z)))) Initial program 8.0
rmApplied add-cube-cbrt8.6
Applied *-un-lft-identity8.6
Applied times-frac8.6
Applied associate-*r*5.3
Simplified5.3
rmApplied add-cube-cbrt5.1
Applied times-frac5.1
Applied associate-*l*4.8
if -1.2753538938243992e-306 < (+ x (* (- y z) (/ (- t x) (- a z)))) < 0.0Initial program 61.8
Taylor expanded around inf 24.8
Final simplification7.4
herbie shell --seed 2020053
(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)))))