\sqrt{x - 1} \cdot \sqrt{x}\left(x + \frac{-1}{2}\right) + \frac{\frac{-1}{8}}{x}double f(double x) {
double r62034 = x;
double r62035 = 1.0;
double r62036 = r62034 - r62035;
double r62037 = sqrt(r62036);
double r62038 = sqrt(r62034);
double r62039 = r62037 * r62038;
return r62039;
}
double f(double x) {
double r62040 = x;
double r62041 = -0.5;
double r62042 = r62040 + r62041;
double r62043 = -0.125;
double r62044 = r62043 / r62040;
double r62045 = r62042 + r62044;
return r62045;
}



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