\sqrt{0.5 \cdot \left(1 + \frac{q}{\sqrt{p + q \cdot q}}\right)}\sqrt{0.5 \cdot \left(1 + \frac{q}{\sqrt{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) (p + ((double) (q * q))))))))))))));
}
double code(double q, double p) {
return ((double) sqrt(((double) (0.5 * ((double) (1.0 + ((double) (q / ((double) sqrt(((double) (p + ((double) (q * q))))))))))))));
}



Bits error versus q



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