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



Bits error versus p



Bits error versus x



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