x + \frac{y \cdot \left(z - t\right)}{a}x + \left(z - t\right) \cdot \frac{y}{a}(FPCore (x y z t a) :precision binary64 (+ x (/ (* y (- z t)) a)))
(FPCore (x y z t a) :precision binary64 (+ x (* (- z t) (/ y a))))
double code(double x, double y, double z, double t, double a) {
return x + ((y * (z - t)) / a);
}
double code(double x, double y, double z, double t, double a) {
return x + ((z - t) * (y / a));
}




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.7 |
| Herbie | 2.4 |
Initial program 6.4
rmApplied clear-num_binary64_116716.5
Simplified6.5
rmApplied *-un-lft-identity_binary64_116726.5
Applied times-frac_binary64_116782.5
Applied add-sqr-sqrt_binary64_116942.5
Applied times-frac_binary64_116782.6
Simplified2.6
Simplified2.4
Final simplification2.4
herbie shell --seed 2021050
(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)))