\sqrt{x - 1} \cdot \sqrt{x}\left(x - \frac{0.125}{x}\right) - 0.5double f(double x) {
double r332940 = x;
double r332941 = 1.0;
double r332942 = r332940 - r332941;
double r332943 = sqrt(r332942);
double r332944 = sqrt(r332940);
double r332945 = r332943 * r332944;
return r332945;
}
double f(double x) {
double r332946 = x;
double r332947 = 0.125;
double r332948 = r332947 / r332946;
double r332949 = r332946 - r332948;
double r332950 = 0.5;
double r332951 = r332949 - r332950;
return r332951;
}



Bits error versus x
Results
Initial program 0.5
Taylor expanded around inf 0.5
Simplified0.5
Final simplification0.5
herbie shell --seed 2019179 +o rules:numerics
(FPCore (x)
:name "sqrt times"
(* (sqrt (- x 1.0)) (sqrt x)))