\frac{x}{y} \cdot \left(z - t\right) + t\begin{array}{l}
\mathbf{if}\;z \le -2.54272324424329845 \cdot 10^{-9} \lor \neg \left(z \le 1.09387645143986444 \cdot 10^{45}\right):\\
\;\;\;\;\left(\left(\frac{x}{y} \cdot \left(\sqrt[3]{z} \cdot \sqrt[3]{z}\right)\right) \cdot \sqrt[3]{z} + \frac{x}{y} \cdot \left(-t\right)\right) + t\\
\mathbf{else}:\\
\;\;\;\;x \cdot \frac{z - t}{y} + t\\
\end{array}double f(double x, double y, double z, double t) {
double r435404 = x;
double r435405 = y;
double r435406 = r435404 / r435405;
double r435407 = z;
double r435408 = t;
double r435409 = r435407 - r435408;
double r435410 = r435406 * r435409;
double r435411 = r435410 + r435408;
return r435411;
}
double f(double x, double y, double z, double t) {
double r435412 = z;
double r435413 = -2.5427232442432984e-09;
bool r435414 = r435412 <= r435413;
double r435415 = 1.0938764514398644e+45;
bool r435416 = r435412 <= r435415;
double r435417 = !r435416;
bool r435418 = r435414 || r435417;
double r435419 = x;
double r435420 = y;
double r435421 = r435419 / r435420;
double r435422 = cbrt(r435412);
double r435423 = r435422 * r435422;
double r435424 = r435421 * r435423;
double r435425 = r435424 * r435422;
double r435426 = t;
double r435427 = -r435426;
double r435428 = r435421 * r435427;
double r435429 = r435425 + r435428;
double r435430 = r435429 + r435426;
double r435431 = r435412 - r435426;
double r435432 = r435431 / r435420;
double r435433 = r435419 * r435432;
double r435434 = r435433 + r435426;
double r435435 = r435418 ? r435430 : r435434;
return r435435;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 1.9 |
|---|---|
| Target | 2.2 |
| Herbie | 3.4 |
if z < -2.5427232442432984e-09 or 1.0938764514398644e+45 < z Initial program 1.5
rmApplied sub-neg1.5
Applied distribute-lft-in1.5
rmApplied add-cube-cbrt2.0
Applied associate-*r*2.0
if -2.5427232442432984e-09 < z < 1.0938764514398644e+45Initial program 2.2
rmApplied div-inv2.3
Applied associate-*l*4.4
Simplified4.4
Final simplification3.4
herbie shell --seed 2020081
(FPCore (x y z t)
:name "Numeric.Signal.Multichannel:$cget from hsignal-0.2.7.1"
:precision binary64
:herbie-target
(if (< z 2.759456554562692e-282) (+ (* (/ x y) (- z t)) t) (if (< z 2.326994450874436e-110) (+ (* x (/ (- z t) y)) t) (+ (* (/ x y) (- z t)) t)))
(+ (* (/ x y) (- z t)) t))