x + \frac{y \cdot \left(z - x\right)}{t}\begin{array}{l}
\mathbf{if}\;x \le -4.1794824165367897 \cdot 10^{-126}:\\
\;\;\;\;x + \frac{y}{t} \cdot \left(z - x\right)\\
\mathbf{elif}\;x \le -5.76861790673687925 \cdot 10^{-285}:\\
\;\;\;\;x + \frac{1}{t \cdot \frac{1}{y \cdot \left(z - x\right)}}\\
\mathbf{elif}\;x \le 9.5149408658066749 \cdot 10^{-56}:\\
\;\;\;\;x + \frac{y}{\sqrt[3]{t} \cdot \sqrt[3]{t}} \cdot \frac{z - x}{\sqrt[3]{t}}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y}{t} \cdot \left(z - x\right)\\
\end{array}double code(double x, double y, double z, double t) {
return ((double) (x + ((double) (((double) (y * ((double) (z - x)))) / t))));
}
double code(double x, double y, double z, double t) {
double VAR;
if ((x <= -4.17948241653679e-126)) {
VAR = ((double) (x + ((double) (((double) (y / t)) * ((double) (z - x))))));
} else {
double VAR_1;
if ((x <= -5.768617906736879e-285)) {
VAR_1 = ((double) (x + ((double) (1.0 / ((double) (t * ((double) (1.0 / ((double) (y * ((double) (z - x))))))))))));
} else {
double VAR_2;
if ((x <= 9.514940865806675e-56)) {
VAR_2 = ((double) (x + ((double) (((double) (y / ((double) (((double) cbrt(t)) * ((double) cbrt(t)))))) * ((double) (((double) (z - x)) / ((double) cbrt(t))))))));
} else {
VAR_2 = ((double) (x + ((double) (((double) (y / t)) * ((double) (z - x))))));
}
VAR_1 = VAR_2;
}
VAR = VAR_1;
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 6.4 |
|---|---|
| Target | 2.2 |
| Herbie | 1.8 |
if x < -4.17948241653679e-126 or 9.514940865806675e-56 < x Initial program 7.5
rmApplied clear-num7.5
rmApplied associate-/r*0.4
rmApplied div-inv0.5
Applied add-cube-cbrt0.5
Applied times-frac0.6
Simplified0.5
Simplified0.5
if -4.17948241653679e-126 < x < -5.768617906736879e-285Initial program 4.3
rmApplied clear-num4.3
rmApplied div-inv4.6
if -5.768617906736879e-285 < x < 9.514940865806675e-56Initial program 4.9
rmApplied add-cube-cbrt5.6
Applied times-frac3.4
Final simplification1.8
herbie shell --seed 2020121
(FPCore (x y z t)
:name "Optimisation.CirclePacking:place from circle-packing-0.1.0.4, D"
:precision binary64
:herbie-target
(- x (+ (* x (/ y t)) (* (- z) (/ y t))))
(+ x (/ (* y (- z x)) t)))