\frac{x - y}{z - y} \cdot t\begin{array}{l}
\mathbf{if}\;\frac{x - y}{z - y} \leq 1.9403572670315207 \cdot 10^{+93}:\\
\;\;\;\;\left(\frac{x}{z - y} - \frac{y}{z - y}\right) \cdot t\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(x - y\right) \cdot t}{z - y}\\
\end{array}(FPCore (x y z t) :precision binary64 (* (/ (- x y) (- z y)) t))
(FPCore (x y z t) :precision binary64 (if (<= (/ (- x y) (- z y)) 1.9403572670315207e+93) (* (- (/ x (- z y)) (/ y (- z y))) t) (/ (* (- x y) t) (- z y))))
double code(double x, double y, double z, double t) {
return ((x - y) / (z - y)) * t;
}
double code(double x, double y, double z, double t) {
double tmp;
if (((x - y) / (z - y)) <= 1.9403572670315207e+93) {
tmp = ((x / (z - y)) - (y / (z - y))) * t;
} else {
tmp = ((x - y) * t) / (z - y);
}
return tmp;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 2.1 |
|---|---|
| Target | 2.1 |
| Herbie | 1.9 |
if (/.f64 (-.f64 x y) (-.f64 z y)) < 1.94035726703152075e93Initial program 1.7
rmApplied div-sub_binary64_130411.7
if 1.94035726703152075e93 < (/.f64 (-.f64 x y) (-.f64 z y)) Initial program 7.7
rmApplied associate-*l/_binary64_129794.8
Final simplification1.9
herbie shell --seed 2020352
(FPCore (x y z t)
:name "Numeric.Signal.Multichannel:$cput from hsignal-0.2.7.1"
:precision binary64
:herbie-target
(/ t (/ (- z y) (- x y)))
(* (/ (- x y) (- z y)) t))