x + \frac{y \cdot \left(z - t\right)}{a}\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(z - t\right)\right) + xdouble f(double x, double y, double z, double t, double a) {
double r279507 = x;
double r279508 = y;
double r279509 = z;
double r279510 = t;
double r279511 = r279509 - r279510;
double r279512 = r279508 * r279511;
double r279513 = a;
double r279514 = r279512 / r279513;
double r279515 = r279507 + r279514;
return r279515;
}
double f(double x, double y, double z, double t, double a) {
double r279516 = y;
double r279517 = cbrt(r279516);
double r279518 = r279517 * r279517;
double r279519 = a;
double r279520 = cbrt(r279519);
double r279521 = r279520 * r279520;
double r279522 = r279518 / r279521;
double r279523 = r279517 / r279520;
double r279524 = z;
double r279525 = t;
double r279526 = r279524 - r279525;
double r279527 = r279523 * r279526;
double r279528 = r279522 * r279527;
double r279529 = x;
double r279530 = r279528 + r279529;
return r279530;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 6.1 |
|---|---|
| Target | 0.7 |
| Herbie | 1.0 |
Initial program 6.1
Simplified2.9
rmApplied fma-udef2.9
rmApplied add-cube-cbrt3.4
Applied add-cube-cbrt3.6
Applied times-frac3.6
Applied associate-*l*1.0
Final simplification1.0
herbie shell --seed 2020100 +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)))