\frac{\log \left(1 - x\right)}{\log \left(1 + x\right)}\frac{\mathsf{log1p}\left(x \cdot \left(-x\right)\right)}{\mathsf{log1p}\left(x\right)} - \frac{\mathsf{log1p}\left(x\right)}{\mathsf{log1p}\left(x\right)}double f(double x) {
double r1511036 = 1.0;
double r1511037 = x;
double r1511038 = r1511036 - r1511037;
double r1511039 = log(r1511038);
double r1511040 = r1511036 + r1511037;
double r1511041 = log(r1511040);
double r1511042 = r1511039 / r1511041;
return r1511042;
}
double f(double x) {
double r1511043 = x;
double r1511044 = -r1511043;
double r1511045 = r1511043 * r1511044;
double r1511046 = log1p(r1511045);
double r1511047 = log1p(r1511043);
double r1511048 = r1511046 / r1511047;
double r1511049 = r1511047 / r1511047;
double r1511050 = r1511048 - r1511049;
return r1511050;
}




Bits error versus x
Results
| Original | 60.9 |
|---|---|
| Target | 0.4 |
| Herbie | 0.0 |
Initial program 60.9
Simplified59.9
rmApplied flip--59.9
Applied log-div59.9
Applied div-sub59.9
Simplified59.9
Simplified0.6
rmApplied log1p-expm1-u0.6
Simplified0.0
Final simplification0.0
herbie shell --seed 2019153 +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))))