x + \left(y - z\right) \cdot \frac{t - x}{a - z}\begin{array}{l}
\mathbf{if}\;a \le -1.317914642947996499854808142021174692446 \cdot 10^{-107}:\\
\;\;\;\;x + \frac{y - z}{\sqrt[3]{a - z} \cdot \sqrt[3]{a - z}} \cdot \frac{t - x}{\sqrt[3]{a - z}}\\
\mathbf{elif}\;a \le 1.510853063573210616754360888454383623335 \cdot 10^{-160}:\\
\;\;\;\;\left(\frac{x \cdot y}{z} + t\right) - \frac{t \cdot y}{z}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{\frac{y - z}{\sqrt[3]{a - z} \cdot \sqrt[3]{a - z}}}{\sqrt[3]{\sqrt[3]{a - z} \cdot \sqrt[3]{a - z}}} \cdot \frac{t - x}{\sqrt[3]{\sqrt[3]{a - z}}}\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r112237 = x;
double r112238 = y;
double r112239 = z;
double r112240 = r112238 - r112239;
double r112241 = t;
double r112242 = r112241 - r112237;
double r112243 = a;
double r112244 = r112243 - r112239;
double r112245 = r112242 / r112244;
double r112246 = r112240 * r112245;
double r112247 = r112237 + r112246;
return r112247;
}
double f(double x, double y, double z, double t, double a) {
double r112248 = a;
double r112249 = -1.3179146429479965e-107;
bool r112250 = r112248 <= r112249;
double r112251 = x;
double r112252 = y;
double r112253 = z;
double r112254 = r112252 - r112253;
double r112255 = r112248 - r112253;
double r112256 = cbrt(r112255);
double r112257 = r112256 * r112256;
double r112258 = r112254 / r112257;
double r112259 = t;
double r112260 = r112259 - r112251;
double r112261 = r112260 / r112256;
double r112262 = r112258 * r112261;
double r112263 = r112251 + r112262;
double r112264 = 1.5108530635732106e-160;
bool r112265 = r112248 <= r112264;
double r112266 = r112251 * r112252;
double r112267 = r112266 / r112253;
double r112268 = r112267 + r112259;
double r112269 = r112259 * r112252;
double r112270 = r112269 / r112253;
double r112271 = r112268 - r112270;
double r112272 = cbrt(r112257);
double r112273 = r112258 / r112272;
double r112274 = cbrt(r112256);
double r112275 = r112260 / r112274;
double r112276 = r112273 * r112275;
double r112277 = r112251 + r112276;
double r112278 = r112265 ? r112271 : r112277;
double r112279 = r112250 ? r112263 : r112278;
return r112279;
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus t



Bits error versus a
Results
if a < -1.3179146429479965e-107Initial program 11.9
rmApplied add-cube-cbrt12.4
Applied *-un-lft-identity12.4
Applied times-frac12.4
Applied associate-*r*10.4
Simplified10.4
if -1.3179146429479965e-107 < a < 1.5108530635732106e-160Initial program 25.8
Taylor expanded around inf 14.7
if 1.5108530635732106e-160 < a Initial program 11.2
rmApplied add-cube-cbrt11.7
Applied *-un-lft-identity11.7
Applied times-frac11.7
Applied associate-*r*9.6
Simplified9.6
rmApplied add-cube-cbrt9.7
Applied cbrt-prod9.7
Applied *-un-lft-identity9.7
Applied times-frac9.7
Applied associate-*r*9.5
Simplified9.5
Final simplification11.1
herbie shell --seed 2019326
(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)))))