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 ((((double) (y * ((double) (z - t)))) <= -3.13336732550441e+195)) {
VAR = ((double) (x + ((double) (y * ((double) (((double) (z - t)) / a))))));
} else {
double VAR_1;
if ((((double) (y * ((double) (z - t)))) <= 9.935411130541129e+225)) {
VAR_1 = ((double) (x + ((double) (((double) (y * ((double) (z - t)))) * ((double) (1.0 / a))))));
} else {
VAR_1 = ((double) (x + ((double) (y / ((double) (a / ((double) (z - t))))))));
}
VAR = VAR_1;
}
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 | 0.5 |
if (* y (- z t)) < -3.13336732550441e195Initial program 26.5
Simplified1.0
if -3.13336732550441e195 < (* y (- z t)) < 9.935411130541129e225Initial program 0.4
rmApplied div-inv0.4
if 9.935411130541129e225 < (* y (- z t)) Initial program 34.9
rmApplied associate-/l*0.8
Final simplification0.5
herbie shell --seed 2020192
(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)))