double code(double x, double y, double z, double t, double a) {
return ((double) (x + ((double) (((double) (y * ((double) (z - t)))) / a))));
}
double code(double x, double y, double z, double t, double a) {
double VAR;
if ((y <= 8.997616379224001e-57)) {
VAR = ((double) (x + ((double) (((double) (y / a)) * ((double) (z - t))))));
} else {
VAR = ((double) (x + ((double) (y * ((double) (((double) (z - t)) / a))))));
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 5.9 |
|---|---|
| Target | 0.7 |
| Herbie | 2.3 |
if y < 8.9976163792240008e-57Initial program 4.1
rmApplied associate-/l*6.8
rmApplied associate-/r/2.6
if 8.9976163792240008e-57 < y Initial program 11.4
Simplified1.3
Final simplification2.3
herbie shell --seed 2020190
(FPCore (x y z t a)
:name "Optimisation.CirclePacking:place from circle-packing-0.1.0.4, E"
:precision binary64
:herbie-target
(if (< y -1.0761266216389975e-10) (+ x (/ 1.0 (/ (/ a (- z t)) y))) (if (< y 2.894426862792089e-49) (+ x (/ (* y (- z t)) a)) (+ x (/ y (/ a (- z t))))))
(+ x (/ (* y (- z t)) a)))