x + \frac{y \cdot \left(z - t\right)}{a}x + \frac{\sqrt[3]{y} \cdot \sqrt[3]{y}}{\sqrt[3]{a}} \cdot \left(\frac{z - t}{\sqrt[3]{a}} \cdot \frac{\sqrt[3]{y}}{\sqrt[3]{a}}\right)(FPCore (x y z t a) :precision binary64 (+ x (/ (* y (- z t)) a)))
(FPCore (x y z t a) :precision binary64 (+ x (* (/ (* (cbrt y) (cbrt y)) (cbrt a)) (* (/ (- z t) (cbrt a)) (/ (cbrt y) (cbrt 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 + (((cbrt(y) * cbrt(y)) / cbrt(a)) * (((z - t) / cbrt(a)) * (cbrt(y) / cbrt(a))));
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 6.0 |
|---|---|
| Target | 0.7 |
| Herbie | 2.5 |
Initial program 6.0
rmApplied add-cube-cbrt_binary64_121546.5
Applied times-frac_binary64_121253.3
rmApplied add-cube-cbrt_binary64_121543.4
Applied times-frac_binary64_121253.5
Applied associate-*l*_binary64_120602.5
Simplified2.5
Final simplification2.5
herbie shell --seed 2020289
(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)))