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




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 5.8 |
|---|---|
| Target | 0.6 |
| Herbie | 0.9 |
Initial program 5.8
Simplified2.2
rmApplied add-cube-cbrt_binary64_48872.7
Applied add-cube-cbrt_binary64_48872.8
Applied times-frac_binary64_48582.8
Applied associate-*l*_binary64_47930.9
Final simplification0.9
herbie shell --seed 2021176
(FPCore (x y z t a)
:name "Optimisation.CirclePacking:place from circle-packing-0.1.0.4, F"
: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)))