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.179379256500943923543558174672518347375 \cdot 10^{-293}:\\
\;\;\;\;x + \left(y - z\right) \cdot \frac{t - x}{a - z}\\
\mathbf{elif}\;x + \left(y - z\right) \cdot \frac{t - x}{a - z} \le 4.22759869788422279163979134487375511757 \cdot 10^{-300}:\\
\;\;\;\;y \cdot \left(\frac{x}{z} - \frac{t}{z}\right) + t\\
\mathbf{else}:\\
\;\;\;\;x + \left(\left(\sqrt[3]{\left(y - z\right) \cdot \frac{\sqrt[3]{t - x} \cdot \sqrt[3]{t - x}}{\sqrt[3]{a - z} \cdot \sqrt[3]{a - z}}} \cdot \sqrt[3]{\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 \sqrt[3]{\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}}\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r119701 = x;
double r119702 = y;
double r119703 = z;
double r119704 = r119702 - r119703;
double r119705 = t;
double r119706 = r119705 - r119701;
double r119707 = a;
double r119708 = r119707 - r119703;
double r119709 = r119706 / r119708;
double r119710 = r119704 * r119709;
double r119711 = r119701 + r119710;
return r119711;
}
double f(double x, double y, double z, double t, double a) {
double r119712 = x;
double r119713 = y;
double r119714 = z;
double r119715 = r119713 - r119714;
double r119716 = t;
double r119717 = r119716 - r119712;
double r119718 = a;
double r119719 = r119718 - r119714;
double r119720 = r119717 / r119719;
double r119721 = r119715 * r119720;
double r119722 = r119712 + r119721;
double r119723 = -1.1793792565009439e-293;
bool r119724 = r119722 <= r119723;
double r119725 = 4.227598697884223e-300;
bool r119726 = r119722 <= r119725;
double r119727 = r119712 / r119714;
double r119728 = r119716 / r119714;
double r119729 = r119727 - r119728;
double r119730 = r119713 * r119729;
double r119731 = r119730 + r119716;
double r119732 = cbrt(r119717);
double r119733 = r119732 * r119732;
double r119734 = cbrt(r119719);
double r119735 = r119734 * r119734;
double r119736 = r119733 / r119735;
double r119737 = r119715 * r119736;
double r119738 = cbrt(r119737);
double r119739 = r119738 * r119738;
double r119740 = r119739 * r119738;
double r119741 = r119732 / r119734;
double r119742 = r119740 * r119741;
double r119743 = r119712 + r119742;
double r119744 = r119726 ? r119731 : r119743;
double r119745 = r119724 ? r119722 : r119744;
return r119745;
}



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.1793792565009439e-293Initial program 7.5
if -1.1793792565009439e-293 < (+ x (* (- y z) (/ (- t x) (- a z)))) < 4.227598697884223e-300Initial program 61.2
rmApplied add-cube-cbrt61.0
Applied add-cube-cbrt60.9
Applied times-frac60.9
Applied associate-*r*59.7
rmApplied add-cube-cbrt59.7
Taylor expanded around inf 25.7
Simplified21.1
if 4.227598697884223e-300 < (+ x (* (- y z) (/ (- t x) (- a z)))) Initial program 7.2
rmApplied add-cube-cbrt7.9
Applied add-cube-cbrt8.1
Applied times-frac8.1
Applied associate-*r*4.5
rmApplied add-cube-cbrt4.7
Final simplification8.1
herbie shell --seed 2019298
(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)))))