\sqrt{x - 1} \cdot \sqrt{x}\left(x - \frac{0.125}{x}\right) - 0.5double f(double x) {
double r12073 = x;
double r12074 = 1.0;
double r12075 = r12073 - r12074;
double r12076 = sqrt(r12075);
double r12077 = sqrt(r12073);
double r12078 = r12076 * r12077;
return r12078;
}
double f(double x) {
double r12079 = x;
double r12080 = 0.125;
double r12081 = r12080 / r12079;
double r12082 = r12079 - r12081;
double r12083 = 0.5;
double r12084 = r12082 - r12083;
return r12084;
}



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