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 -5.94158773471745396 \cdot 10^{-307}:\\
\;\;\;\;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 0.0:\\
\;\;\;\;t + y \cdot \left(\frac{x}{z} - \frac{t}{z}\right)\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y - z}{\sqrt[3]{a - z} \cdot \sqrt[3]{a - z}} \cdot \frac{t - x}{\sqrt[3]{a - z}}\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r85506 = x;
double r85507 = y;
double r85508 = z;
double r85509 = r85507 - r85508;
double r85510 = t;
double r85511 = r85510 - r85506;
double r85512 = a;
double r85513 = r85512 - r85508;
double r85514 = r85511 / r85513;
double r85515 = r85509 * r85514;
double r85516 = r85506 + r85515;
return r85516;
}
double f(double x, double y, double z, double t, double a) {
double r85517 = x;
double r85518 = y;
double r85519 = z;
double r85520 = r85518 - r85519;
double r85521 = t;
double r85522 = r85521 - r85517;
double r85523 = a;
double r85524 = r85523 - r85519;
double r85525 = r85522 / r85524;
double r85526 = r85520 * r85525;
double r85527 = r85517 + r85526;
double r85528 = -5.941587734717454e-307;
bool r85529 = r85527 <= r85528;
double r85530 = 0.0;
bool r85531 = r85527 <= r85530;
double r85532 = r85517 / r85519;
double r85533 = r85521 / r85519;
double r85534 = r85532 - r85533;
double r85535 = r85518 * r85534;
double r85536 = r85521 + r85535;
double r85537 = cbrt(r85524);
double r85538 = r85537 * r85537;
double r85539 = r85520 / r85538;
double r85540 = r85522 / r85537;
double r85541 = r85539 * r85540;
double r85542 = r85517 + r85541;
double r85543 = r85531 ? r85536 : r85542;
double r85544 = r85529 ? r85527 : r85543;
return r85544;
}



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)))) < -5.941587734717454e-307Initial program 6.9
if -5.941587734717454e-307 < (+ x (* (- y z) (/ (- t x) (- a z)))) < 0.0Initial program 61.8
rmApplied add-cube-cbrt61.4
Taylor expanded around inf 24.7
Simplified18.9
if 0.0 < (+ x (* (- y z) (/ (- t x) (- a z)))) Initial program 7.7
rmApplied add-cube-cbrt8.4
Applied *-un-lft-identity8.4
Applied times-frac8.4
Applied associate-*r*5.0
Simplified5.0
Final simplification7.6
herbie shell --seed 2020047
(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)))))