2 \cdot \sqrt{\left(x \cdot y + x \cdot z\right) + y \cdot z}\sqrt{y \cdot z + x \cdot \left(y + z\right)} \cdot 2double f(double x, double y, double z) {
double r518517 = 2.0;
double r518518 = x;
double r518519 = y;
double r518520 = r518518 * r518519;
double r518521 = z;
double r518522 = r518518 * r518521;
double r518523 = r518520 + r518522;
double r518524 = r518519 * r518521;
double r518525 = r518523 + r518524;
double r518526 = sqrt(r518525);
double r518527 = r518517 * r518526;
return r518527;
}
double f(double x, double y, double z) {
double r518528 = y;
double r518529 = z;
double r518530 = r518528 * r518529;
double r518531 = x;
double r518532 = r518528 + r518529;
double r518533 = r518531 * r518532;
double r518534 = r518530 + r518533;
double r518535 = sqrt(r518534);
double r518536 = 2.0;
double r518537 = r518535 * r518536;
return r518537;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 20.3 |
|---|---|
| Target | 19.1 |
| Herbie | 20.3 |
Initial program 20.3
Simplified20.3
Final simplification20.3
herbie shell --seed 2019199
(FPCore (x y z)
:name "Diagrams.TwoD.Apollonian:descartes from diagrams-contrib-1.3.0.5"
:herbie-target
(if (< z 7.636950090573675e+176) (* 2.0 (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.0))
(* 2.0 (sqrt (+ (+ (* x y) (* x z)) (* y z)))))