x + \frac{y \cdot \left(z - t\right)}{a}x + \frac{1}{\frac{\frac{a}{y}}{z - t}}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 + (1.0 / ((a / y) / (z - t))));
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 6.3 |
|---|---|
| Target | 0.7 |
| Herbie | 2.5 |
Initial program 6.3
rmApplied *-commutative6.3
Applied associate-/l*2.5
rmApplied clear-num2.5
Final simplification2.5
herbie shell --seed 2020071 +o rules:numerics
(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 (/ (/ a (- z t)) y))) (if (< y 2.894426862792089e-49) (+ x (/ (* y (- z t)) a)) (+ x (/ y (/ a (- z t))))))
(+ x (/ (* y (- z t)) a)))