double code(double x) {
return ((double) log(((double) (1.0 + x))));
}
double code(double x) {
double VAR;
if ((((double) (1.0 + x)) <= 1.0)) {
VAR = ((double) (((double) log(1.0)) + ((double) (x * ((double) (1.0 - ((double) (((double) (x * 0.5)) / ((double) (1.0 * 1.0))))))))));
} else {
VAR = ((double) (((double) log(((double) sqrt(((double) (1.0 + x)))))) + ((double) log(((double) sqrt(((double) (1.0 + x))))))));
}
return VAR;
}




Bits error versus x
Results
| Original | 39.5 |
|---|---|
| Target | 0.2 |
| Herbie | 0.6 |
if (+ 1.0 x) < 1Initial program 59.6
Taylor expanded around 0 0.3
Simplified0.3
if 1 < (+ 1.0 x) Initial program 1.1
rmApplied add-sqr-sqrt1.2
Applied log-prod1.2
Final simplification0.6
herbie shell --seed 2020181
(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)))