\frac{x - y}{z - y} \cdot t\begin{array}{l}
\mathbf{if}\;\frac{x - y}{z - y} \cdot t = -\infty:\\
\;\;\;\;\frac{1}{z - y} \cdot \left(t \cdot \left(x - y\right)\right)\\
\mathbf{elif}\;\frac{x - y}{z - y} \cdot t \le -2.640233754742455 \cdot 10^{-299}:\\
\;\;\;\;\frac{x - y}{z - y} \cdot t\\
\mathbf{elif}\;\frac{x - y}{z - y} \cdot t \le 0.0:\\
\;\;\;\;\frac{1}{z - y} \cdot \left(t \cdot \left(x - y\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{t}{\frac{z}{x - y} - \frac{y}{x - y}}\\
\end{array}double f(double x, double y, double z, double t) {
double r21267123 = x;
double r21267124 = y;
double r21267125 = r21267123 - r21267124;
double r21267126 = z;
double r21267127 = r21267126 - r21267124;
double r21267128 = r21267125 / r21267127;
double r21267129 = t;
double r21267130 = r21267128 * r21267129;
return r21267130;
}
double f(double x, double y, double z, double t) {
double r21267131 = x;
double r21267132 = y;
double r21267133 = r21267131 - r21267132;
double r21267134 = z;
double r21267135 = r21267134 - r21267132;
double r21267136 = r21267133 / r21267135;
double r21267137 = t;
double r21267138 = r21267136 * r21267137;
double r21267139 = -inf.0;
bool r21267140 = r21267138 <= r21267139;
double r21267141 = 1.0;
double r21267142 = r21267141 / r21267135;
double r21267143 = r21267137 * r21267133;
double r21267144 = r21267142 * r21267143;
double r21267145 = -2.640233754742455e-299;
bool r21267146 = r21267138 <= r21267145;
double r21267147 = 0.0;
bool r21267148 = r21267138 <= r21267147;
double r21267149 = r21267134 / r21267133;
double r21267150 = r21267132 / r21267133;
double r21267151 = r21267149 - r21267150;
double r21267152 = r21267137 / r21267151;
double r21267153 = r21267148 ? r21267144 : r21267152;
double r21267154 = r21267146 ? r21267138 : r21267153;
double r21267155 = r21267140 ? r21267144 : r21267154;
return r21267155;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 2.3 |
|---|---|
| Target | 2.3 |
| Herbie | 1.0 |
if (* (/ (- x y) (- z y)) t) < -inf.0 or -2.640233754742455e-299 < (* (/ (- x y) (- z y)) t) < 0.0Initial program 12.1
rmApplied clear-num12.1
rmApplied associate-*l/11.7
Simplified11.7
rmApplied div-inv11.7
Applied *-un-lft-identity11.7
Applied times-frac0.5
Simplified0.5
if -inf.0 < (* (/ (- x y) (- z y)) t) < -2.640233754742455e-299Initial program 0.3
if 0.0 < (* (/ (- x y) (- z y)) t) Initial program 1.8
rmApplied clear-num2.0
rmApplied associate-*l/1.8
Simplified1.8
rmApplied div-sub1.8
Final simplification1.0
herbie shell --seed 2019168
(FPCore (x y z t)
:name "Numeric.Signal.Multichannel:$cput from hsignal-0.2.7.1"
:herbie-target
(/ t (/ (- z y) (- x y)))
(* (/ (- x y) (- z y)) t))