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 r653316 = 2.0;
double r653317 = x;
double r653318 = y;
double r653319 = r653317 * r653318;
double r653320 = z;
double r653321 = r653317 * r653320;
double r653322 = r653319 + r653321;
double r653323 = r653318 * r653320;
double r653324 = r653322 + r653323;
double r653325 = sqrt(r653324);
double r653326 = r653316 * r653325;
return r653326;
}
double f(double x, double y, double z) {
double r653327 = 2.0;
double r653328 = x;
double r653329 = y;
double r653330 = r653328 * r653329;
double r653331 = z;
double r653332 = r653328 * r653331;
double r653333 = r653330 + r653332;
double r653334 = r653329 * r653331;
double r653335 = r653333 + r653334;
double r653336 = sqrt(r653335);
double r653337 = r653327 * r653336;
return r653337;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 19.4 |
|---|---|
| Target | 18.8 |
| Herbie | 19.4 |
Initial program 19.4
Final simplification19.4
herbie shell --seed 2020062 +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)))))