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 r859846 = 2.0;
double r859847 = x;
double r859848 = y;
double r859849 = r859847 * r859848;
double r859850 = z;
double r859851 = r859847 * r859850;
double r859852 = r859849 + r859851;
double r859853 = r859848 * r859850;
double r859854 = r859852 + r859853;
double r859855 = sqrt(r859854);
double r859856 = r859846 * r859855;
return r859856;
}
double f(double x, double y, double z) {
double r859857 = 2.0;
double r859858 = x;
double r859859 = y;
double r859860 = r859858 * r859859;
double r859861 = z;
double r859862 = r859858 * r859861;
double r859863 = r859860 + r859862;
double r859864 = r859859 * r859861;
double r859865 = r859863 + r859864;
double r859866 = sqrt(r859865);
double r859867 = r859857 * r859866;
return r859867;
}




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
Final simplification19.8
herbie shell --seed 2020042 +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)))))