x + \left(y - z\right) \cdot \frac{t - x}{a - z}\begin{array}{l}
\mathbf{if}\;z \le -2.7916675455541999 \cdot 10^{190}:\\
\;\;\;\;t + y \cdot \left(\frac{x}{z} - \frac{t}{z}\right)\\
\mathbf{elif}\;z \le -7.29199372430577621 \cdot 10^{117}:\\
\;\;\;\;x + \frac{y - z}{-\left(a - z\right)} \cdot \left(-\left(t - x\right)\right)\\
\mathbf{elif}\;z \le -7.9425005692074541 \cdot 10^{65}:\\
\;\;\;\;t + y \cdot \left(\frac{x}{z} - \frac{t}{z}\right)\\
\mathbf{elif}\;z \le 2.50608572105253218 \cdot 10^{82}:\\
\;\;\;\;x + \left(\left(y - z\right) \cdot \frac{\sqrt[3]{t - x} \cdot \sqrt[3]{t - x}}{\sqrt[3]{a - z} \cdot \sqrt[3]{a - z}}\right) \cdot \frac{\sqrt[3]{t - x}}{\sqrt[3]{a - z}}\\
\mathbf{else}:\\
\;\;\;\;t + y \cdot \left(\frac{x}{z} - \frac{t}{z}\right)\\
\end{array}double code(double x, double y, double z, double t, double a) {
return ((double) (x + ((double) (((double) (y - z)) * ((double) (((double) (t - x)) / ((double) (a - z))))))));
}
double code(double x, double y, double z, double t, double a) {
double VAR;
if ((z <= -2.7916675455542e+190)) {
VAR = ((double) (t + ((double) (y * ((double) (((double) (x / z)) - ((double) (t / z))))))));
} else {
double VAR_1;
if ((z <= -7.291993724305776e+117)) {
VAR_1 = ((double) (x + ((double) (((double) (((double) (y - z)) / ((double) -(((double) (a - z)))))) * ((double) -(((double) (t - x))))))));
} else {
double VAR_2;
if ((z <= -7.942500569207454e+65)) {
VAR_2 = ((double) (t + ((double) (y * ((double) (((double) (x / z)) - ((double) (t / z))))))));
} else {
double VAR_3;
if ((z <= 2.5060857210525322e+82)) {
VAR_3 = ((double) (x + ((double) (((double) (((double) (y - z)) * ((double) (((double) (((double) cbrt(((double) (t - x)))) * ((double) cbrt(((double) (t - x)))))) / ((double) (((double) cbrt(((double) (a - z)))) * ((double) cbrt(((double) (a - z)))))))))) * ((double) (((double) cbrt(((double) (t - x)))) / ((double) cbrt(((double) (a - z))))))))));
} else {
VAR_3 = ((double) (t + ((double) (y * ((double) (((double) (x / z)) - ((double) (t / z))))))));
}
VAR_2 = VAR_3;
}
VAR_1 = VAR_2;
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus t



Bits error versus a
Results
if z < -2.7916675455542e+190 or -7.291993724305776e+117 < z < -7.942500569207454e+65 or 2.5060857210525322e+82 < z Initial program 24.7
rmApplied add-cube-cbrt25.2
Applied add-cube-cbrt25.4
Applied times-frac25.4
Applied associate-*r*21.6
Taylor expanded around inf 26.6
Simplified18.3
if -2.7916675455542e+190 < z < -7.291993724305776e+117Initial program 21.3
rmApplied clear-num21.7
Applied un-div-inv21.6
rmApplied frac-2neg21.6
Applied associate-/r/18.6
if -7.942500569207454e+65 < z < 2.5060857210525322e+82Initial program 8.2
rmApplied add-cube-cbrt8.7
Applied add-cube-cbrt8.8
Applied times-frac8.8
Applied associate-*r*5.4
Final simplification10.9
herbie shell --seed 2020113
(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)))))