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



Bits error versus q



Bits error versus r



Bits error versus p
Results
Initial program 29.9
Final simplification29.9
herbie shell --seed 2020152
(FPCore (q r p)
:name "(sqrt (+ 0.5 (/ (* 0.5 (- q r)) (sqrt (+ (* p p) (* (- q r) (- q r)))))))"
:precision binary64
(sqrt (+ 0.5 (/ (* 0.5 (- q r)) (sqrt (+ (* p p) (* (- q r) (- q r))))))))