\frac{\log \left(1 - x\right)}{\log \left(1 + x\right)}\frac{\mathsf{log1p}\left(-x\right)}{\mathsf{log1p}\left(x\right)}double f(double x) {
double r5914010 = 1.0;
double r5914011 = x;
double r5914012 = r5914010 - r5914011;
double r5914013 = log(r5914012);
double r5914014 = r5914010 + r5914011;
double r5914015 = log(r5914014);
double r5914016 = r5914013 / r5914015;
return r5914016;
}
double f(double x) {
double r5914017 = x;
double r5914018 = -r5914017;
double r5914019 = log1p(r5914018);
double r5914020 = log1p(r5914017);
double r5914021 = r5914019 / r5914020;
return r5914021;
}




Bits error versus x
Results
| Original | 60.9 |
|---|---|
| Target | 0.3 |
| Herbie | 0.0 |
Initial program 60.9
Simplified60.0
rmApplied log1p-expm1-u60.0
Simplified0.0
Final simplification0.0
herbie shell --seed 2019146 +o rules:numerics
(FPCore (x)
:name "qlog (example 3.10)"
:pre (and (< -1 x) (< x 1))
:herbie-target
(- (+ (+ (+ 1 x) (/ (* x x) 2)) (* 5/12 (pow x 3))))
(/ (log (- 1 x)) (log (+ 1 x))))