\sqrt{x + 1} - \sqrt{x}\sqrt{1} \cdot \frac{\sqrt{1}}{\sqrt{x + 1} + \sqrt{x}}double f(double x) {
double r111030 = x;
double r111031 = 1.0;
double r111032 = r111030 + r111031;
double r111033 = sqrt(r111032);
double r111034 = sqrt(r111030);
double r111035 = r111033 - r111034;
return r111035;
}
double f(double x) {
double r111036 = 1.0;
double r111037 = sqrt(r111036);
double r111038 = x;
double r111039 = r111038 + r111036;
double r111040 = sqrt(r111039);
double r111041 = sqrt(r111038);
double r111042 = r111040 + r111041;
double r111043 = r111037 / r111042;
double r111044 = r111037 * r111043;
return r111044;
}




Bits error versus x
Results
| Original | 29.7 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
Initial program 29.7
rmApplied flip--29.5
Simplified0.2
rmApplied *-un-lft-identity0.2
Applied add-sqr-sqrt0.2
Applied times-frac0.2
Simplified0.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2020049 +o rules:numerics
(FPCore (x)
:name "2sqrt (example 3.1)"
:precision binary64
:herbie-target
(/ 1 (+ (sqrt (+ x 1)) (sqrt x)))
(- (sqrt (+ x 1)) (sqrt x)))