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



Bits error versus q



Bits error versus p
Results
Initial program 22.6
Final simplification22.6
herbie shell --seed 2020153
(FPCore (q p)
:name "(sqrt (* 0.5 (+ 1 (/ q (sqrt (+ (* p p) (* q q)))))))"
:precision binary64
(sqrt (* 0.5 (+ 1.0 (/ q (sqrt (+ (* p p) (* q q))))))))