\frac{\log \left(1 - x\right)}{\log \left(1 + x\right)}
\frac{1}{\frac{\mathsf{log1p}\left(x\right)}{\mathsf{log1p}\left(-x\right)}}
(FPCore (x) :precision binary64 (/ (log (- 1.0 x)) (log (+ 1.0 x))))
(FPCore (x) :precision binary64 (/ 1.0 (/ (log1p x) (log1p (- x)))))
double code(double x) {
return log(1.0 - x) / log(1.0 + x);
}
double code(double x) {
return 1.0 / (log1p(x) / log1p(-x));
}




Bits error versus x
Results
| Original | 61.4 |
|---|---|
| Target | 0.4 |
| Herbie | 0.0 |
Initial program 61.4
Simplified0.0
Applied clear-num_binary640.0
Final simplification0.0
herbie shell --seed 2022077
(FPCore (x)
:name "qlog (example 3.10)"
:precision binary64
:pre (and (< -1.0 x) (< x 1.0))
:herbie-target
(- (+ (+ (+ 1.0 x) (/ (* x x) 2.0)) (* 0.4166666666666667 (pow x 3.0))))
(/ (log (- 1.0 x)) (log (+ 1.0 x))))