\left(x - \frac{y}{z \cdot 3}\right) + \frac{t}{\left(z \cdot 3\right) \cdot y}\begin{array}{l}
\mathbf{if}\;t \le -2.20900859555046478 \cdot 10^{76} \lor \neg \left(t \le 3.4334795866523095 \cdot 10^{-75}\right):\\
\;\;\;\;\left(x - \frac{y}{z \cdot 3}\right) + 0.333333333333333315 \cdot \frac{t}{z \cdot y}\\
\mathbf{else}:\\
\;\;\;\;\left(x - \frac{\frac{y}{z}}{3}\right) + \frac{\frac{t}{y}}{z \cdot 3}\\
\end{array}double code(double x, double y, double z, double t) {
return ((double) (((double) (x - ((double) (y / ((double) (z * 3.0)))))) + ((double) (t / ((double) (((double) (z * 3.0)) * y))))));
}
double code(double x, double y, double z, double t) {
double VAR;
if (((t <= -2.2090085955504648e+76) || !(t <= 3.4334795866523095e-75))) {
VAR = ((double) (((double) (x - ((double) (y / ((double) (z * 3.0)))))) + ((double) (0.3333333333333333 * ((double) (t / ((double) (z * y))))))));
} else {
VAR = ((double) (((double) (x - ((double) (((double) (y / z)) / 3.0)))) + ((double) (((double) (t / y)) / ((double) (z * 3.0))))));
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 3.6 |
|---|---|
| Target | 1.8 |
| Herbie | 0.7 |
if t < -2.20900859555046478e76 or 3.4334795866523095e-75 < t Initial program 1.0
Taylor expanded around 0 1.0
if -2.20900859555046478e76 < t < 3.4334795866523095e-75Initial program 5.5
rmApplied add-cube-cbrt5.6
Applied times-frac0.9
rmApplied associate-*l/0.6
Simplified0.4
rmApplied associate-/r*0.4
Final simplification0.7
herbie shell --seed 2020173
(FPCore (x y z t)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, H"
:precision binary64
:herbie-target
(+ (- x (/ y (* z 3.0))) (/ (/ t (* z 3.0)) y))
(+ (- x (/ y (* z 3.0))) (/ t (* (* z 3.0) y))))