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.9756495134823586e-122) || !(t <= 8.566572551246772e+57))) {
VAR = ((double) (((double) (x + ((double) (((double) (y / 3.0)) * ((double) (-1.0 / z)))))) + ((double) (t / ((double) (y * ((double) (z * 3.0))))))));
} else {
VAR = ((double) (x + ((double) (((double) (((double) (t / y)) - 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.6 |
if t < -1.9756495134823586e-122 or 8.56657255124677238e57 < t Initial program 1.0
rmApplied *-un-lft-identity1.0
Applied times-frac1.0
if -1.9756495134823586e-122 < t < 8.56657255124677238e57Initial program 6.0
Simplified0.3
Final simplification0.6
herbie shell --seed 2020191
(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))))