\sqrt{0.5 \cdot \left(1 + \frac{x}{\sqrt{10^{5} + x \cdot x}}\right)}\sqrt{0.5 \cdot \left(1 + \frac{x}{\sqrt{10^{5} + x \cdot x}}\right)}double code(double x) {
return ((double) sqrt(((double) (0.5 * ((double) (1.0 + ((double) (x / ((double) sqrt(((double) (100000.0 + ((double) (x * x))))))))))))));
}
double code(double x) {
return ((double) sqrt(((double) (0.5 * ((double) (1.0 + ((double) (x / ((double) sqrt(((double) (100000.0 + ((double) (x * x))))))))))))));
}



Bits error versus x
Results
Initial program 21.6
Final simplification21.6
herbie shell --seed 2020153
(FPCore (x)
:name "(sqrt (* 0.5 (+ 1 (/ x (sqrt (+ 100000.0 (* x x)))))))"
:precision binary64
(sqrt (* 0.5 (+ 1.0 (/ x (sqrt (+ 100000.0 (* x x))))))))