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 r692028 = 2.0;
double r692029 = x;
double r692030 = y;
double r692031 = r692029 * r692030;
double r692032 = z;
double r692033 = r692029 * r692032;
double r692034 = r692031 + r692033;
double r692035 = r692030 * r692032;
double r692036 = r692034 + r692035;
double r692037 = sqrt(r692036);
double r692038 = r692028 * r692037;
return r692038;
}
double f(double x, double y, double z) {
double r692039 = 2.0;
double r692040 = x;
double r692041 = y;
double r692042 = r692040 * r692041;
double r692043 = z;
double r692044 = r692040 * r692043;
double r692045 = r692042 + r692044;
double r692046 = r692041 * r692043;
double r692047 = r692045 + r692046;
double r692048 = sqrt(r692047);
double r692049 = r692039 * r692048;
return r692049;
}




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)))))