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



Bits error versus a



Bits error versus b
Results
Initial program 0.0
Final simplification0.0
herbie shell --seed 2020152
(FPCore (a b)
:name "(sqrt (/ (+ (- a b) 1) (+ a b)))"
:precision binary64
(sqrt (/ (+ (- a b) 1.0) (+ a b))))