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 r107946 = x;
double r107947 = y;
double r107948 = z;
double r107949 = r107947 - r107948;
double r107950 = t;
double r107951 = r107950 - r107946;
double r107952 = a;
double r107953 = r107952 - r107948;
double r107954 = r107951 / r107953;
double r107955 = r107949 * r107954;
double r107956 = r107946 + r107955;
return r107956;
}
double f(double x, double y, double z, double t, double a) {
double r107957 = x;
double r107958 = y;
double r107959 = z;
double r107960 = r107958 - r107959;
double r107961 = t;
double r107962 = r107961 - r107957;
double r107963 = a;
double r107964 = r107963 - r107959;
double r107965 = r107962 / r107964;
double r107966 = r107960 * r107965;
double r107967 = r107957 + r107966;
double r107968 = -5.941587734717454e-307;
bool r107969 = r107967 <= r107968;
double r107970 = 0.0;
bool r107971 = r107967 <= r107970;
double r107972 = r107957 / r107959;
double r107973 = r107961 / r107959;
double r107974 = r107972 - r107973;
double r107975 = r107958 * r107974;
double r107976 = r107961 + r107975;
double r107977 = cbrt(r107964);
double r107978 = r107977 * r107977;
double r107979 = r107960 / r107978;
double r107980 = r107962 / r107977;
double r107981 = r107979 * r107980;
double r107982 = r107957 + r107981;
double r107983 = r107971 ? r107976 : r107982;
double r107984 = r107969 ? r107967 : r107983;
return r107984;
}



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)))))