x + \left(y - z\right) \cdot \frac{t - x}{a - z}\begin{array}{l}
\mathbf{if}\;a \leq -3.980018115154435 \cdot 10^{-169} \lor \neg \left(a \leq 3.528776701234468 \cdot 10^{-87}\right):\\
\;\;\;\;x + \frac{y - z}{a - z} \cdot \left(t - x\right)\\
\mathbf{else}:\\
\;\;\;\;t + y \cdot \left(\frac{x}{z} - \frac{t}{z}\right)\\
\end{array}(FPCore (x y z t a) :precision binary64 (+ x (* (- y z) (/ (- t x) (- a z)))))
(FPCore (x y z t a) :precision binary64 (if (or (<= a -3.980018115154435e-169) (not (<= a 3.528776701234468e-87))) (+ x (* (/ (- y z) (- a z)) (- t x))) (+ t (* y (- (/ x z) (/ t z))))))
double code(double x, double y, double z, double t, double a) {
return ((double) (x + ((double) (((double) (y - z)) * (((double) (t - x)) / ((double) (a - z)))))));
}
double code(double x, double y, double z, double t, double a) {
double VAR;
if (((a <= -3.980018115154435e-169) || !(a <= 3.528776701234468e-87))) {
VAR = ((double) (x + ((double) ((((double) (y - z)) / ((double) (a - z))) * ((double) (t - x))))));
} else {
VAR = ((double) (t + ((double) (y * ((double) ((x / z) - (t / z)))))));
}
return VAR;
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus t



Bits error versus a
Results
if a < -3.9800181151544353e-169 or 3.528776701234468e-87 < a Initial program 11.3
rmApplied clear-num11.5
rmApplied associate-/r/11.4
Applied associate-*r*9.0
Simplified8.9
if -3.9800181151544353e-169 < a < 3.528776701234468e-87Initial program 24.9
Taylor expanded around inf 14.8
Simplified13.2
Final simplification10.0
herbie shell --seed 2020198
(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)))))