\frac{1 - \sqrt{1 - \left(4 \cdot a\right) \cdot b}}{2}\frac{1 - \sqrt{1 - \left(4 \cdot a\right) \cdot b}}{2}double code(double a, double b) {
return ((double) (((double) (1.0 - ((double) sqrt(((double) (1.0 - ((double) (((double) (4.0 * a)) * b)))))))) / 2.0));
}
double code(double a, double b) {
return ((double) (((double) (1.0 - ((double) sqrt(((double) (1.0 - ((double) (((double) (4.0 * a)) * b)))))))) / 2.0));
}



Bits error versus a



Bits error versus b
Results
Initial program 34.6
Final simplification34.6
herbie shell --seed 2020152
(FPCore (a b)
:name "(/ (- 1.0 (sqrt (- 1.0 (* (* 4.0 a) b)))) 2.0)"
:precision binary64
(/ (- 1.0 (sqrt (- 1.0 (* (* 4.0 a) b)))) 2.0))