\sqrt{0.5 \cdot \left(1 + \frac{x}{\sqrt{a + x \cdot x}}\right)}\sqrt{0.5 \cdot \left(1 + \frac{x}{\sqrt{a + x \cdot x}}\right)}double code(double x, double a) {
return ((double) sqrt(((double) (0.5 * ((double) (1.0 + ((double) (x / ((double) sqrt(((double) (a + ((double) (x * x))))))))))))));
}
double code(double x, double a) {
return ((double) sqrt(((double) (0.5 * ((double) (1.0 + ((double) (x / ((double) sqrt(((double) (a + ((double) (x * x))))))))))))));
}



Bits error versus x



Bits error versus a
Results
Initial program 22.8
Final simplification22.8
herbie shell --seed 2020152
(FPCore (x a)
:name "(sqrt (* 0.5 (+ 1 (/ x (sqrt (+ a (* x x)))))))"
:precision binary64
(sqrt (* 0.5 (+ 1.0 (/ x (sqrt (+ a (* x x))))))))