x - \frac{y \cdot \left(z - t\right)}{a}x - \frac{\sqrt[3]{y} \cdot \sqrt[3]{y}}{\sqrt[3]{\sqrt[3]{a} \cdot \sqrt[3]{a}} \cdot \sqrt[3]{a}} \cdot \frac{\frac{\sqrt[3]{y}}{\sqrt[3]{a}}}{\frac{\sqrt[3]{\sqrt[3]{a}}}{z - t}}double f(double x, double y, double z, double t, double a) {
double r21768695 = x;
double r21768696 = y;
double r21768697 = z;
double r21768698 = t;
double r21768699 = r21768697 - r21768698;
double r21768700 = r21768696 * r21768699;
double r21768701 = a;
double r21768702 = r21768700 / r21768701;
double r21768703 = r21768695 - r21768702;
return r21768703;
}
double f(double x, double y, double z, double t, double a) {
double r21768704 = x;
double r21768705 = y;
double r21768706 = cbrt(r21768705);
double r21768707 = r21768706 * r21768706;
double r21768708 = a;
double r21768709 = cbrt(r21768708);
double r21768710 = r21768709 * r21768709;
double r21768711 = cbrt(r21768710);
double r21768712 = r21768711 * r21768709;
double r21768713 = r21768707 / r21768712;
double r21768714 = r21768706 / r21768709;
double r21768715 = cbrt(r21768709);
double r21768716 = z;
double r21768717 = t;
double r21768718 = r21768716 - r21768717;
double r21768719 = r21768715 / r21768718;
double r21768720 = r21768714 / r21768719;
double r21768721 = r21768713 * r21768720;
double r21768722 = r21768704 - r21768721;
return r21768722;
}




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.2 |
Initial program 6.1
rmApplied associate-/l*5.5
rmApplied *-un-lft-identity5.5
Applied add-cube-cbrt6.0
Applied times-frac6.0
Applied associate-/r*2.8
Simplified2.8
rmApplied *-un-lft-identity2.8
Applied add-cube-cbrt2.9
Applied cbrt-prod2.9
Applied times-frac2.9
Applied add-cube-cbrt3.0
Applied times-frac3.0
Applied times-frac1.2
Simplified1.2
Final simplification1.2
herbie shell --seed 2019171
(FPCore (x y z t a)
:name "Optimisation.CirclePacking:place from circle-packing-0.1.0.4, F"
: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)))