x + \frac{y \cdot \left(z - t\right)}{a}y \cdot \left(\frac{z}{a} - \frac{t}{a}\right) + xdouble f(double x, double y, double z, double t, double a) {
double r267227 = x;
double r267228 = y;
double r267229 = z;
double r267230 = t;
double r267231 = r267229 - r267230;
double r267232 = r267228 * r267231;
double r267233 = a;
double r267234 = r267232 / r267233;
double r267235 = r267227 + r267234;
return r267235;
}
double f(double x, double y, double z, double t, double a) {
double r267236 = y;
double r267237 = z;
double r267238 = a;
double r267239 = r267237 / r267238;
double r267240 = t;
double r267241 = r267240 / r267238;
double r267242 = r267239 - r267241;
double r267243 = r267236 * r267242;
double r267244 = x;
double r267245 = r267243 + r267244;
return r267245;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 5.9 |
|---|---|
| Target | 0.8 |
| Herbie | 6.2 |
if y < -3.079662010669181e+61 or 1.737882189922714e-20 < y Initial program 15.8
Taylor expanded around 0 15.8
Simplified3.9
Taylor expanded around 0 15.8
Simplified1.0
if -3.079662010669181e+61 < y < 1.737882189922714e-20Initial program 0.7
Taylor expanded around 0 0.7
Final simplification6.2
herbie shell --seed 2019291
(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.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)))