x + \frac{y \cdot \left(z - x\right)}{t}\begin{array}{l}
\mathbf{if}\;t \le -1.37836980272270247 \cdot 10^{107} \lor \neg \left(t \le 9.89123025060424723 \cdot 10^{-138}\right):\\
\;\;\;\;x + \frac{\sqrt[3]{y} \cdot \sqrt[3]{y}}{\sqrt[3]{t}} \cdot \left(\frac{\sqrt[3]{y}}{\sqrt[3]{t}} \cdot \frac{z - x}{\sqrt[3]{t}}\right)\\
\mathbf{else}:\\
\;\;\;\;x + \left(y \cdot \left(z - x\right)\right) \cdot \frac{1}{t}\\
\end{array}double code(double x, double y, double z, double t) {
return ((double) (x + (((double) (y * ((double) (z - x)))) / t)));
}
double code(double x, double y, double z, double t) {
double VAR;
if (((t <= -1.3783698027227025e+107) || !(t <= 9.891230250604247e-138))) {
VAR = ((double) (x + ((double) ((((double) (((double) cbrt(y)) * ((double) cbrt(y)))) / ((double) cbrt(t))) * ((double) ((((double) cbrt(y)) / ((double) cbrt(t))) * (((double) (z - x)) / ((double) cbrt(t)))))))));
} else {
VAR = ((double) (x + ((double) (((double) (y * ((double) (z - x)))) * (1.0 / t)))));
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 6.7 |
|---|---|
| Target | 2.1 |
| Herbie | 1.9 |
if t < -1.37836980272270247e107 or 9.89123025060424723e-138 < t Initial program 8.7
rmApplied add-cube-cbrt9.1
Applied times-frac1.6
rmApplied add-cube-cbrt1.7
Applied times-frac1.7
Applied associate-*l*1.1
if -1.37836980272270247e107 < t < 9.89123025060424723e-138Initial program 3.1
rmApplied div-inv3.1
Final simplification1.9
herbie shell --seed 2020182
(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)) (* (neg z) (/ y t))))
(+ x (/ (* y (- z x)) t)))