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.3
Final simplification29.3
herbie shell --seed 2020153
(FPCore (b a)
:name "(- b (sqrt (- (* b b) a)))"
:precision binary64
(- b (sqrt (- (* b b) a))))