x - \frac{y \cdot \left(z - t\right)}{a}x - \frac{\frac{y}{\frac{\sqrt[3]{a}}{\sqrt[3]{z - t}} \cdot \frac{\sqrt[3]{a}}{\sqrt[3]{z - t}}}}{\frac{\sqrt[3]{a}}{\sqrt[3]{z - t}}}double f(double x, double y, double z, double t, double a) {
double r20431572 = x;
double r20431573 = y;
double r20431574 = z;
double r20431575 = t;
double r20431576 = r20431574 - r20431575;
double r20431577 = r20431573 * r20431576;
double r20431578 = a;
double r20431579 = r20431577 / r20431578;
double r20431580 = r20431572 - r20431579;
return r20431580;
}
double f(double x, double y, double z, double t, double a) {
double r20431581 = x;
double r20431582 = y;
double r20431583 = a;
double r20431584 = cbrt(r20431583);
double r20431585 = z;
double r20431586 = t;
double r20431587 = r20431585 - r20431586;
double r20431588 = cbrt(r20431587);
double r20431589 = r20431584 / r20431588;
double r20431590 = r20431589 * r20431589;
double r20431591 = r20431582 / r20431590;
double r20431592 = r20431591 / r20431589;
double r20431593 = r20431581 - r20431592;
return r20431593;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 6.3 |
|---|---|
| Target | 0.6 |
| Herbie | 1.5 |
Initial program 6.3
rmApplied associate-/l*5.8
rmApplied add-cube-cbrt6.2
Applied add-cube-cbrt6.4
Applied times-frac6.4
Applied associate-/r*1.5
Simplified1.5
Final simplification1.5
herbie shell --seed 2019174 +o rules:numerics
(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)))