b - \sqrt{b \cdot b - a}b - \sqrt{b \cdot b - a}double code(double b, double a) {
return ((double) (b - ((double) sqrt(((double) (((double) (b * b)) - a))))));
}
double code(double b, double a) {
return ((double) (b - ((double) sqrt(((double) (((double) (b * b)) - a))))));
}



Bits error versus b



Bits error versus a
Results
Initial program 29.5
Final simplification29.5
herbie shell --seed 2020152
(FPCore (b a)
:name "(- b (sqrt (- (* b b) a)))"
:precision binary64
(- b (sqrt (- (* b b) a))))