\sqrt{\frac{\left(x - y\right) + 1}{x + y}}\sqrt{\frac{\left(x - y\right) + 1}{x + y}}double code(double x, double y) {
return ((double) sqrt(((double) (((double) (((double) (x - y)) + 1.0)) / ((double) (x + y))))));
}
double code(double x, double y) {
return ((double) sqrt(((double) (((double) (((double) (x - y)) + 1.0)) / ((double) (x + y))))));
}



Bits error versus x



Bits error versus y
Results
Initial program 0.1
Final simplification0.1
herbie shell --seed 2020153
(FPCore (x y)
:name "(sqrt (/ (+ (- x y) 1) (+ x y)))"
:precision binary64
(sqrt (/ (+ (- x y) 1.0) (+ x y))))