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 r645060 = 2.0;
double r645061 = x;
double r645062 = y;
double r645063 = r645061 * r645062;
double r645064 = z;
double r645065 = r645061 * r645064;
double r645066 = r645063 + r645065;
double r645067 = r645062 * r645064;
double r645068 = r645066 + r645067;
double r645069 = sqrt(r645068);
double r645070 = r645060 * r645069;
return r645070;
}
double f(double x, double y, double z) {
double r645071 = 2.0;
double r645072 = x;
double r645073 = y;
double r645074 = r645072 * r645073;
double r645075 = z;
double r645076 = r645072 * r645075;
double r645077 = r645074 + r645076;
double r645078 = r645073 * r645075;
double r645079 = r645077 + r645078;
double r645080 = sqrt(r645079);
double r645081 = r645071 * r645080;
return r645081;
}




Bits error versus x




Bits error versus y




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