2 \cdot \sqrt{\left(x \cdot y + x \cdot z\right) + y \cdot z}2 \cdot \sqrt{x \cdot y + z \cdot \left(x + y\right)}double f(double x, double y, double z) {
double r386290 = 2.0;
double r386291 = x;
double r386292 = y;
double r386293 = r386291 * r386292;
double r386294 = z;
double r386295 = r386291 * r386294;
double r386296 = r386293 + r386295;
double r386297 = r386292 * r386294;
double r386298 = r386296 + r386297;
double r386299 = sqrt(r386298);
double r386300 = r386290 * r386299;
return r386300;
}
double f(double x, double y, double z) {
double r386301 = 2.0;
double r386302 = x;
double r386303 = y;
double r386304 = r386302 * r386303;
double r386305 = z;
double r386306 = r386302 + r386303;
double r386307 = r386305 * r386306;
double r386308 = r386304 + r386307;
double r386309 = sqrt(r386308);
double r386310 = r386301 * r386309;
return r386310;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 19.8 |
|---|---|
| Target | 19.2 |
| Herbie | 19.8 |
Initial program 19.8
rmApplied associate-+l+19.8
Simplified19.8
Final simplification19.8
herbie shell --seed 2019323 +o rules:numerics
(FPCore (x y z)
:name "Diagrams.TwoD.Apollonian:descartes from diagrams-contrib-1.3.0.5"
:precision binary64
:herbie-target
(if (< z 7.636950090573675e+176) (* 2 (sqrt (+ (* (+ x y) z) (* x y)))) (* (* (+ (* 0.25 (* (* (pow y -0.75) (* (pow z -0.75) x)) (+ y z))) (* (pow z 0.25) (pow y 0.25))) (+ (* 0.25 (* (* (pow y -0.75) (* (pow z -0.75) x)) (+ y z))) (* (pow z 0.25) (pow y 0.25)))) 2))
(* 2 (sqrt (+ (+ (* x y) (* x z)) (* y z)))))