\frac{c - \sqrt{c \cdot d}}{c - d}\frac{c - \sqrt{c \cdot d}}{c - d}double code(double c, double d) {
return ((double) (((double) (c - ((double) sqrt(((double) (c * d)))))) / ((double) (c - d))));
}
double code(double c, double d) {
return ((double) (((double) (c - ((double) sqrt(((double) (c * d)))))) / ((double) (c - d))));
}



Bits error versus c



Bits error versus d
Results
Initial program 12.2
Final simplification12.2
herbie shell --seed 2020152
(FPCore (c d)
:name "(/ (- c (sqrt (* c d))) (- c d))"
:precision binary64
(/ (- c (sqrt (* c d))) (- c d)))