\frac{x - y}{z - y} \cdot t\begin{array}{l}
\mathbf{if}\;\frac{x - y}{z - y} \le 2.470328229206232720882843964341106861825 \cdot 10^{-323}:\\
\;\;\;\;\frac{\frac{t}{\sqrt[3]{z - y}} \cdot \frac{x - y}{\sqrt[3]{z - y}}}{\sqrt[3]{z - y}}\\
\mathbf{elif}\;\frac{x - y}{z - y} \le 5.779702072575541191674228138407204633316 \cdot 10^{148}:\\
\;\;\;\;\left(\frac{1}{z - y} \cdot \left(x - y\right)\right) \cdot t\\
\mathbf{else}:\\
\;\;\;\;\left(x - y\right) \cdot \frac{t}{z - y}\\
\end{array}double f(double x, double y, double z, double t) {
double r446367 = x;
double r446368 = y;
double r446369 = r446367 - r446368;
double r446370 = z;
double r446371 = r446370 - r446368;
double r446372 = r446369 / r446371;
double r446373 = t;
double r446374 = r446372 * r446373;
return r446374;
}
double f(double x, double y, double z, double t) {
double r446375 = x;
double r446376 = y;
double r446377 = r446375 - r446376;
double r446378 = z;
double r446379 = r446378 - r446376;
double r446380 = r446377 / r446379;
double r446381 = 2.4703282292062e-323;
bool r446382 = r446380 <= r446381;
double r446383 = t;
double r446384 = cbrt(r446379);
double r446385 = r446383 / r446384;
double r446386 = r446377 / r446384;
double r446387 = r446385 * r446386;
double r446388 = r446387 / r446384;
double r446389 = 5.779702072575541e+148;
bool r446390 = r446380 <= r446389;
double r446391 = 1.0;
double r446392 = r446391 / r446379;
double r446393 = r446392 * r446377;
double r446394 = r446393 * r446383;
double r446395 = r446383 / r446379;
double r446396 = r446377 * r446395;
double r446397 = r446390 ? r446394 : r446396;
double r446398 = r446382 ? r446388 : r446397;
return r446398;
}




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.9 |
if (/ (- x y) (- z y)) < 2.4703282292062e-323Initial program 4.9
rmApplied add-cube-cbrt5.7
Applied *-un-lft-identity5.7
Applied times-frac5.7
rmApplied associate-*r/5.7
Applied associate-*l/4.7
Simplified4.8
if 2.4703282292062e-323 < (/ (- x y) (- z y)) < 5.779702072575541e+148Initial program 0.1
rmApplied div-inv0.3
if 5.779702072575541e+148 < (/ (- x y) (- z y)) Initial program 12.2
rmApplied div-inv12.3
Applied associate-*l*2.4
Simplified2.3
Final simplification1.9
herbie shell --seed 2019194
(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))