\sqrt{x + 1} - \sqrt{x}\mathsf{expm1}\left(\mathsf{log1p}\left(\frac{1}{\sqrt{x + 1} + \sqrt{x}}\right)\right)double f(double x) {
double r462556 = x;
double r462557 = 1.0;
double r462558 = r462556 + r462557;
double r462559 = sqrt(r462558);
double r462560 = sqrt(r462556);
double r462561 = r462559 - r462560;
return r462561;
}
double f(double x) {
double r462562 = 1.0;
double r462563 = x;
double r462564 = r462563 + r462562;
double r462565 = sqrt(r462564);
double r462566 = sqrt(r462563);
double r462567 = r462565 + r462566;
double r462568 = r462562 / r462567;
double r462569 = log1p(r462568);
double r462570 = expm1(r462569);
return r462570;
}




Bits error versus x
Results
| Original | 30.3 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
Initial program 30.3
Simplified30.3
rmApplied flip--30.1
Simplified0.2
rmApplied expm1-log1p-u0.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2019179 +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)))