\left(x - \frac{y}{z \cdot 3}\right) + \frac{t}{\left(z \cdot 3\right) \cdot y}\begin{array}{l}
\mathbf{if}\;t \le -1.11239943139134638 \cdot 10^{-76}:\\
\;\;\;\;\left(x - y \cdot \frac{1}{z \cdot 3}\right) + \frac{t}{\left(z \cdot 3\right) \cdot y}\\
\mathbf{elif}\;t \le 2.5459723469470038 \cdot 10^{-216}:\\
\;\;\;\;\left(x - \frac{y}{z \cdot 3}\right) + \frac{1}{z \cdot 3} \cdot \frac{t}{y}\\
\mathbf{else}:\\
\;\;\;\;\left(x - \frac{y}{z \cdot 3}\right) + \frac{\frac{t}{z}}{3 \cdot y}\\
\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 <= -1.1123994313913464e-76)) {
VAR = ((double) (((double) (x - ((double) (y * ((double) (1.0 / ((double) (z * 3.0)))))))) + ((double) (t / ((double) (((double) (z * 3.0)) * y))))));
} else {
double VAR_1;
if ((t <= 2.5459723469470038e-216)) {
VAR_1 = ((double) (((double) (x - ((double) (y / ((double) (z * 3.0)))))) + ((double) (((double) (1.0 / ((double) (z * 3.0)))) * ((double) (t / y))))));
} else {
VAR_1 = ((double) (((double) (x - ((double) (y / ((double) (z * 3.0)))))) + ((double) (((double) (t / z)) / ((double) (3.0 * y))))));
}
VAR = VAR_1;
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 3.7 |
|---|---|
| Target | 1.8 |
| Herbie | 1.1 |
if t < -1.1123994313913464e-76Initial program 0.8
rmApplied div-inv0.8
if -1.1123994313913464e-76 < t < 2.5459723469470038e-216Initial program 7.0
rmApplied *-un-lft-identity7.0
Applied times-frac0.2
if 2.5459723469470038e-216 < t Initial program 3.1
rmApplied associate-*l*3.1
rmApplied associate-/r*2.0
Final simplification1.1
herbie shell --seed 2020148
(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))))