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



Bits error versus p



Bits error versus q



Bits error versus r
Results
Initial program 27.3
Final simplification27.3
herbie shell --seed 2020153
(FPCore (p q r)
:name "(/ p (sqrt (* 0.5 (+ (+ (* (* 4 p) p) (* (- q r) (- q r))) (* (- q r) q)))))"
:precision binary64
(/ p (sqrt (* 0.5 (+ (+ (* (* 4.0 p) p) (* (- q r) (- q r))) (* (- q r) q))))))