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.58469193808142249 \cdot 10^{-283} \lor \neg \left(x + \left(y - z\right) \cdot \frac{t - x}{a - z} \le 7.4730959852087275 \cdot 10^{-283}\right):\\
\;\;\;\;\mathsf{fma}\left(\sqrt[3]{x} \cdot \sqrt[3]{x}, \sqrt[3]{\sqrt[3]{x} \cdot \sqrt[3]{x}} \cdot \sqrt[3]{\sqrt[3]{x}}, \frac{y - z}{\sqrt[3]{a - z} \cdot \sqrt[3]{a - z}} \cdot \frac{t - x}{\sqrt[3]{a - z}}\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{x}{z}, y, t - \frac{t \cdot y}{z}\right)\\
\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.5846919380814225e-283) || !((x + ((y - z) * ((t - x) / (a - z)))) <= 7.4730959852087275e-283))) {
temp = fma((cbrt(x) * cbrt(x)), (cbrt((cbrt(x) * cbrt(x))) * cbrt(cbrt(x))), (((y - z) / (cbrt((a - z)) * cbrt((a - z)))) * ((t - x) / cbrt((a - z)))));
} else {
temp = fma((x / z), y, (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.5846919380814225e-283 or 7.4730959852087275e-283 < (+ x (* (- y z) (/ (- t x) (- a z)))) Initial program 7.6
rmApplied add-cube-cbrt8.2
Applied *-un-lft-identity8.2
Applied times-frac8.2
Applied associate-*r*5.6
Simplified5.6
rmApplied add-cube-cbrt6.0
Applied fma-def6.0
rmApplied add-cube-cbrt6.0
Applied cbrt-prod6.1
if -1.5846919380814225e-283 < (+ x (* (- y z) (/ (- t x) (- a z)))) < 7.4730959852087275e-283Initial program 59.2
Taylor expanded around inf 27.8
Simplified22.8
Final simplification8.5
herbie shell --seed 2020058 +o rules:numerics
(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)))))