x + \frac{y \cdot \left(z - t\right)}{a}x + \frac{\sqrt[3]{y} \cdot \sqrt[3]{y}}{\frac{\sqrt[3]{a} \cdot \sqrt[3]{a}}{\sqrt[3]{z - t} \cdot \sqrt[3]{z - t}}} \cdot \frac{\sqrt[3]{y}}{\frac{\sqrt[3]{a}}{\sqrt[3]{z - t}}}(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 (- z t)) (cbrt (- z t)))))
(/ (cbrt y) (/ (cbrt a) (cbrt (- 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 + (((cbrt(y) * cbrt(y)) / ((cbrt(a) * cbrt(a)) / (cbrt(z - t) * cbrt(z - t)))) * (cbrt(y) / (cbrt(a) / cbrt(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.2 |
|---|---|
| Target | 0.7 |
| Herbie | 1.5 |
Initial program 6.2
rmApplied associate-/l*_binary64_109355.8
rmApplied add-cube-cbrt_binary64_110256.2
Applied add-cube-cbrt_binary64_110256.3
Applied times-frac_binary64_109966.3
Applied add-cube-cbrt_binary64_110256.4
Applied times-frac_binary64_109961.5
Final simplification1.5
herbie shell --seed 2020356
(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)))