\frac{x - y}{z - y} \cdot t\begin{array}{l}
\mathbf{if}\;\frac{x - y}{z - y} \le -1.477269956526729882709689256541591977713 \cdot 10^{-242}:\\
\;\;\;\;1 \cdot \frac{t}{\frac{z - y}{x - y}}\\
\mathbf{elif}\;\frac{x - y}{z - y} \le -0.0:\\
\;\;\;\;1 \cdot \left(\frac{t}{z - y} \cdot \left(x - y\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{x - y}{z - y} \cdot t\\
\end{array}double f(double x, double y, double z, double t) {
double r563264 = x;
double r563265 = y;
double r563266 = r563264 - r563265;
double r563267 = z;
double r563268 = r563267 - r563265;
double r563269 = r563266 / r563268;
double r563270 = t;
double r563271 = r563269 * r563270;
return r563271;
}
double f(double x, double y, double z, double t) {
double r563272 = x;
double r563273 = y;
double r563274 = r563272 - r563273;
double r563275 = z;
double r563276 = r563275 - r563273;
double r563277 = r563274 / r563276;
double r563278 = -1.4772699565267299e-242;
bool r563279 = r563277 <= r563278;
double r563280 = 1.0;
double r563281 = t;
double r563282 = r563276 / r563274;
double r563283 = r563281 / r563282;
double r563284 = r563280 * r563283;
double r563285 = -0.0;
bool r563286 = r563277 <= r563285;
double r563287 = r563281 / r563276;
double r563288 = r563287 * r563274;
double r563289 = r563280 * r563288;
double r563290 = r563277 * r563281;
double r563291 = r563286 ? r563289 : r563290;
double r563292 = r563279 ? r563284 : r563291;
return r563292;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 2.4 |
|---|---|
| Target | 2.4 |
| Herbie | 1.4 |
if (/ (- x y) (- z y)) < -1.4772699565267299e-242Initial program 2.7
rmApplied clear-num2.8
rmApplied *-un-lft-identity2.8
Applied *-un-lft-identity2.8
Applied times-frac2.8
Applied *-un-lft-identity2.8
Applied times-frac2.8
Applied associate-*l*2.8
Simplified2.5
if -1.4772699565267299e-242 < (/ (- x y) (- z y)) < -0.0Initial program 13.3
rmApplied clear-num14.8
rmApplied *-un-lft-identity14.8
Applied *-un-lft-identity14.8
Applied times-frac14.8
Applied *-un-lft-identity14.8
Applied times-frac14.8
Applied associate-*l*14.8
Simplified14.7
rmApplied associate-/r/0.2
if -0.0 < (/ (- x y) (- z y)) Initial program 1.6
Final simplification1.4
herbie shell --seed 2020001
(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))