\frac{\log \left(1 - x\right)}{\log \left(1 + x\right)}\frac{\mathsf{log1p}\left(\left(-x\right)\right)}{\mathsf{log1p}\left(x\right)}double f(double x) {
double r2302449 = 1.0;
double r2302450 = x;
double r2302451 = r2302449 - r2302450;
double r2302452 = log(r2302451);
double r2302453 = r2302449 + r2302450;
double r2302454 = log(r2302453);
double r2302455 = r2302452 / r2302454;
return r2302455;
}
double f(double x) {
double r2302456 = x;
double r2302457 = -r2302456;
double r2302458 = log1p(r2302457);
double r2302459 = log1p(r2302456);
double r2302460 = r2302458 / r2302459;
return r2302460;
}




Bits error versus x
Results
| Original | 60.9 |
|---|---|
| Target | 0.3 |
| Herbie | 0.0 |
Initial program 60.9
Simplified60.0
rmApplied sub-neg60.0
Applied log1p-def0.0
rmApplied expm1-log1p-u0.0
Applied log1p-expm10.0
Final simplification0.0
herbie shell --seed 2019130 +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))))