s + \frac{x - {s}^{2}}{2 \cdot s}s + \frac{x - {s}^{2}}{2 \cdot s}double code(double s, double x) {
return ((double) (s + ((double) (((double) (x - ((double) pow(s, 2.0)))) / ((double) (2.0 * s))))));
}
double code(double s, double x) {
return ((double) (s + ((double) (((double) (x - ((double) pow(s, 2.0)))) / ((double) (2.0 * s))))));
}



Bits error versus s



Bits error versus x
Results
Initial program 18.4
Final simplification18.4
herbie shell --seed 2020152
(FPCore (s x)
:name "(+ s (/ (- x (pow s 2)) (* 2 s)))"
:precision binary64
(+ s (/ (- x (pow s 2.0)) (* 2.0 s))))