x - \frac{y \cdot \left(z - t\right)}{a}x - \frac{1}{\frac{\frac{\sqrt[3]{a} \cdot \sqrt[3]{a}}{\sqrt[3]{y} \cdot \sqrt[3]{y}}}{\frac{z - t}{\frac{\sqrt[3]{a}}{\sqrt[3]{y}}}}}double f(double x, double y, double z, double t, double a) {
double r365719 = x;
double r365720 = y;
double r365721 = z;
double r365722 = t;
double r365723 = r365721 - r365722;
double r365724 = r365720 * r365723;
double r365725 = a;
double r365726 = r365724 / r365725;
double r365727 = r365719 - r365726;
return r365727;
}
double f(double x, double y, double z, double t, double a) {
double r365728 = x;
double r365729 = 1.0;
double r365730 = a;
double r365731 = cbrt(r365730);
double r365732 = r365731 * r365731;
double r365733 = y;
double r365734 = cbrt(r365733);
double r365735 = r365734 * r365734;
double r365736 = r365732 / r365735;
double r365737 = z;
double r365738 = t;
double r365739 = r365737 - r365738;
double r365740 = r365731 / r365734;
double r365741 = r365739 / r365740;
double r365742 = r365736 / r365741;
double r365743 = r365729 / r365742;
double r365744 = r365728 - r365743;
return r365744;
}




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.6 |
| Herbie | 0.9 |
Initial program 6.0
rmApplied sub-neg6.0
Applied distribute-lft-in6.0
rmApplied clear-num6.1
Simplified2.6
rmApplied add-cube-cbrt3.0
Applied add-cube-cbrt3.2
Applied times-frac3.2
Applied associate-/l*0.9
Final simplification0.9
herbie shell --seed 2020049
(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 (/ (/ a (- z t)) y))) (if (< y 2.894426862792089e-49) (- x (/ (* y (- z t)) a)) (- x (/ y (/ a (- z t))))))
(- x (/ (* y (- z t)) a)))