\frac{\log \left(1 - x\right)}{\log \left(1 + x\right)}\mathsf{log1p}\left(\mathsf{expm1}\left(\frac{\mathsf{log1p}\left(-x\right)}{\mathsf{log1p}\left(x\right)}\right)\right)double f(double x) {
double r2112985 = 1.0;
double r2112986 = x;
double r2112987 = r2112985 - r2112986;
double r2112988 = log(r2112987);
double r2112989 = r2112985 + r2112986;
double r2112990 = log(r2112989);
double r2112991 = r2112988 / r2112990;
return r2112991;
}
double f(double x) {
double r2112992 = x;
double r2112993 = -r2112992;
double r2112994 = log1p(r2112993);
double r2112995 = log1p(r2112992);
double r2112996 = r2112994 / r2112995;
double r2112997 = expm1(r2112996);
double r2112998 = log1p(r2112997);
return r2112998;
}




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