double code(double x) {
return ((double) log(((double) (1.0 + x))));
}
double code(double x) {
double VAR;
if ((((double) (1.0 + x)) <= 1.0000054502705635)) {
VAR = ((double) (((double) (1.0 * x)) + ((double) (((double) log(1.0)) + ((double) (((double) (((double) (-0.5 / ((double) (1.0 * 1.0)))) * ((double) (2.0 * ((double) (x * ((double) log(((double) cbrt(((double) exp(x)))))))))))) + ((double) (x * ((double) (((double) (-0.5 / ((double) (1.0 * 1.0)))) * ((double) (x * 0.3333333333333333))))))))))));
} else {
VAR = ((double) log(((double) (1.0 + x))));
}
return VAR;
}




Bits error versus x
Results
| Original | 39.6 |
|---|---|
| Target | 0.2 |
| Herbie | 0.3 |
if (+ 1.0 x) < 1.00000545027056353Initial program 59.1
Taylor expanded around 0 0.3
Simplified0.3
rmApplied add-log-exp0.7
Simplified0.7
rmApplied add-cube-cbrt0.7
Applied unpow-prod-down0.7
Applied unpow-prod-down0.7
Applied log-prod0.7
Simplified0.7
Simplified0.3
rmApplied pow1/30.3
Applied log-pow0.3
Simplified0.3
if 1.00000545027056353 < (+ 1.0 x) Initial program 0.1
Final simplification0.3
herbie shell --seed 2020190
(FPCore (x)
:name "ln(1 + x)"
:precision binary64
:herbie-target
(if (== (+ 1.0 x) 1.0) x (/ (* x (log (+ 1.0 x))) (- (+ 1.0 x) 1.0)))
(log (+ 1.0 x)))