x - \frac{y \cdot \left(z - t\right)}{a}x - \frac{y \cdot \left(z - t\right)}{a}double f(double x, double y, double z, double t, double a) {
double r245567 = x;
double r245568 = y;
double r245569 = z;
double r245570 = t;
double r245571 = r245569 - r245570;
double r245572 = r245568 * r245571;
double r245573 = a;
double r245574 = r245572 / r245573;
double r245575 = r245567 - r245574;
return r245575;
}
double f(double x, double y, double z, double t, double a) {
double r245576 = x;
double r245577 = y;
double r245578 = z;
double r245579 = t;
double r245580 = r245578 - r245579;
double r245581 = r245577 * r245580;
double r245582 = a;
double r245583 = r245581 / r245582;
double r245584 = r245576 - r245583;
return r245584;
}




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 | 6.1 |
if (/ (* y (- z t)) a) < -2.5229463530234854e+190 or 1.3317847127302016e-49 < (/ (* y (- z t)) a) Initial program 14.0
rmApplied clear-num14.1
rmApplied sub-neg14.1
Simplified3.0
if -2.5229463530234854e+190 < (/ (* y (- z t)) a) < 1.3317847127302016e-49Initial program 0.4
rmApplied div-inv0.5
Final simplification6.1
herbie shell --seed 2019303
(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.07612662163899753e-10) (- x (/ 1 (/ (/ a (- z t)) y))) (if (< y 2.8944268627920891e-49) (- x (/ (* y (- z t)) a)) (- x (/ y (/ a (- z t))))))
(- x (/ (* y (- z t)) a)))