\frac{x}{y} \cdot \left(z - t\right) + t\begin{array}{l}
\mathbf{if}\;y \le -6.87411635940552805 \cdot 10^{59}:\\
\;\;\;\;\frac{\sqrt[3]{x} \cdot \sqrt[3]{x}}{1} \cdot \left(\frac{\sqrt[3]{x}}{y} \cdot \left(z - t\right)\right) + t\\
\mathbf{elif}\;y \le 1.2846808360778829 \cdot 10^{-79}:\\
\;\;\;\;\left(\frac{x \cdot z}{y} - \frac{t \cdot x}{y}\right) + t\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y} \cdot \left(z - t\right) + t\\
\end{array}double code(double x, double y, double z, double t) {
return ((double) (((double) (((double) (x / y)) * ((double) (z - t)))) + t));
}
double code(double x, double y, double z, double t) {
double VAR;
if ((y <= -6.874116359405528e+59)) {
VAR = ((double) (((double) (((double) (((double) (((double) cbrt(x)) * ((double) cbrt(x)))) / 1.0)) * ((double) (((double) (((double) cbrt(x)) / y)) * ((double) (z - t)))))) + t));
} else {
double VAR_1;
if ((y <= 1.2846808360778829e-79)) {
VAR_1 = ((double) (((double) (((double) (((double) (x * z)) / y)) - ((double) (((double) (t * x)) / y)))) + t));
} else {
VAR_1 = ((double) (((double) (((double) (x / y)) * ((double) (z - t)))) + t));
}
VAR = VAR_1;
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 2.2 |
|---|---|
| Target | 2.4 |
| Herbie | 1.6 |
if y < -6.87411635940552805e59Initial program 1.3
rmApplied *-un-lft-identity1.3
Applied add-cube-cbrt1.7
Applied times-frac1.7
Applied associate-*l*1.1
if -6.87411635940552805e59 < y < 1.2846808360778829e-79Initial program 3.3
rmApplied add-cube-cbrt4.1
Applied associate-*l*4.1
Taylor expanded around 0 2.3
if 1.2846808360778829e-79 < y Initial program 1.1
Final simplification1.6
herbie shell --seed 2020162
(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))