\sqrt{0.5 \cdot \left(1 - \frac{q}{\sqrt{\left(4 \cdot p\right) \cdot p + q \cdot q}}\right)}\sqrt{0.5 \cdot \left(1 - \frac{q}{\sqrt{\left(4 \cdot p\right) \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) (((double) (4.0 * 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) (((double) (4.0 * p)) * p)) + ((double) (q * q))))))))))))));
}



Bits error versus q



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