2 \cdot \sqrt{\left(x \cdot y + x \cdot z\right) + y \cdot z}2 \cdot \sqrt{\left(x \cdot y + x \cdot z\right) + y \cdot z}double f(double x, double y, double z) {
double r674258 = 2.0;
double r674259 = x;
double r674260 = y;
double r674261 = r674259 * r674260;
double r674262 = z;
double r674263 = r674259 * r674262;
double r674264 = r674261 + r674263;
double r674265 = r674260 * r674262;
double r674266 = r674264 + r674265;
double r674267 = sqrt(r674266);
double r674268 = r674258 * r674267;
return r674268;
}
double f(double x, double y, double z) {
double r674269 = 2.0;
double r674270 = x;
double r674271 = y;
double r674272 = r674270 * r674271;
double r674273 = z;
double r674274 = r674270 * r674273;
double r674275 = r674272 + r674274;
double r674276 = r674271 * r674273;
double r674277 = r674275 + r674276;
double r674278 = sqrt(r674277);
double r674279 = r674269 * r674278;
return r674279;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 19.8 |
|---|---|
| Target | 18.9 |
| Herbie | 19.8 |
Initial program 19.8
Final simplification19.8
herbie shell --seed 2020046
(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)))))