\sqrt{x - 1} \cdot \sqrt{x}x + \left(\frac{-1}{2} - \frac{\frac{1}{8}}{x}\right)double f(double x) {
double r197640 = x;
double r197641 = 1.0;
double r197642 = r197640 - r197641;
double r197643 = sqrt(r197642);
double r197644 = sqrt(r197640);
double r197645 = r197643 * r197644;
return r197645;
}
double f(double x) {
double r197646 = x;
double r197647 = -0.5;
double r197648 = 0.125;
double r197649 = r197648 / r197646;
double r197650 = r197647 - r197649;
double r197651 = r197646 + r197650;
return r197651;
}



Bits error versus x
Results
Initial program 0.5
Taylor expanded around inf 0.3
Simplified0.3
Final simplification0.3
herbie shell --seed 2019156
(FPCore (x)
:name "sqrt times"
(* (sqrt (- x 1)) (sqrt x)))