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) (z - t)) <= -7.801339917478246e+140)) {
VAR = ((double) (x + ((double) (((double) (z - t)) * ((double) (y / a))))));
} else {
double VAR_1;
if ((((double) (z - t)) <= 7.123466189650817e-09)) {
VAR_1 = ((double) (x + ((double) (y * ((double) (((double) (z - t)) / a))))));
} else {
double VAR_2;
if ((((double) (z - t)) <= 3.383832292408153e+254)) {
VAR_2 = ((double) (x + ((double) (((double) (z - t)) * ((double) (y / a))))));
} else {
VAR_2 = ((double) (x + ((double) (((double) (((double) cbrt(y)) * ((double) (((double) cbrt(y)) * ((double) sqrt(((double) (z - t)))))))) * ((double) (((double) sqrt(((double) (z - t)))) * ((double) (((double) cbrt(y)) / a))))))));
}
VAR_1 = VAR_2;
}
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 | 6.4 |
|---|---|
| Target | 0.8 |
| Herbie | 1.9 |
if (- z t) < -7.80133991747824628e140 or 7.12346618965081666e-9 < (- z t) < 3.3838322924081531e254Initial program 9.3
rmApplied associate-/l*8.3
rmApplied associate-/r/1.7
if -7.80133991747824628e140 < (- z t) < 7.12346618965081666e-9Initial program 2.1
Simplified1.9
if 3.3838322924081531e254 < (- z t) Initial program 15.0
rmApplied associate-/l*13.6
rmApplied add-sqr-sqrt13.8
Applied *-un-lft-identity13.8
Applied times-frac13.7
Applied add-cube-cbrt14.1
Applied times-frac6.1
Simplified6.1
Simplified3.5
Final simplification1.9
herbie shell --seed 2020180
(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)))