\sqrt{x + 1} - \sqrt{x}\frac{1}{\sqrt{\sqrt{x + 1}} \cdot \sqrt{\sqrt{x + 1}} + \sqrt{x}}double f(double x) {
double r406004 = x;
double r406005 = 1.0;
double r406006 = r406004 + r406005;
double r406007 = sqrt(r406006);
double r406008 = sqrt(r406004);
double r406009 = r406007 - r406008;
return r406009;
}
double f(double x) {
double r406010 = 1.0;
double r406011 = x;
double r406012 = r406011 + r406010;
double r406013 = sqrt(r406012);
double r406014 = sqrt(r406013);
double r406015 = r406014 * r406014;
double r406016 = sqrt(r406011);
double r406017 = r406015 + r406016;
double r406018 = r406010 / r406017;
return r406018;
}




Bits error versus x
Results
| Original | 30.0 |
|---|---|
| Target | 0.2 |
| Herbie | 0.3 |
Initial program 30.0
rmApplied flip--29.9
Simplified0.2
Simplified0.2
rmApplied add-sqr-sqrt0.2
Applied sqrt-prod0.3
Final simplification0.3
herbie shell --seed 2019194 +o rules:numerics
(FPCore (x)
:name "Main:bigenough3 from C"
:herbie-target
(/ 1.0 (+ (sqrt (+ x 1.0)) (sqrt x)))
(- (sqrt (+ x 1.0)) (sqrt x)))