\frac{\log \left(1 - x\right)}{\log \left(1 + x\right)}-1 - \left(x + \left(0.5 \cdot {x}^{2} + \left(0.4166666666666667 \cdot {x}^{3} + 0.2916666666666667 \cdot {x}^{4}\right)\right)\right)(FPCore (x) :precision binary64 (/ (log (- 1.0 x)) (log (+ 1.0 x))))
(FPCore (x)
:precision binary64
(-
-1.0
(+
x
(+
(* 0.5 (pow x 2.0))
(+
(* 0.4166666666666667 (pow x 3.0))
(* 0.2916666666666667 (pow x 4.0)))))))double code(double x) {
return log(1.0 - x) / log(1.0 + x);
}
double code(double x) {
return -1.0 - (x + ((0.5 * pow(x, 2.0)) + ((0.4166666666666667 * pow(x, 3.0)) + (0.2916666666666667 * pow(x, 4.0)))));
}




Bits error versus x
Results
| Original | 61.4 |
|---|---|
| Target | 0.3 |
| Herbie | 0.3 |
Initial program 61.4
rmApplied flip--_binary6461.0
Applied log-div_binary6461.2
Applied div-sub_binary6461.2
Simplified61.2
Simplified61.2
Taylor expanded around 0 0.3
Final simplification0.3
herbie shell --seed 2021173
(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))))