\log \left(\frac{1 - \varepsilon}{1 + \varepsilon}\right)\left(\varepsilon \cdot -2 - \left(0.6666666666666666 \cdot {\varepsilon}^{3} + 0.4 \cdot {\varepsilon}^{5}\right)\right) - 0.2857142857142857 \cdot {\varepsilon}^{7}(FPCore (eps) :precision binary64 (log (/ (- 1.0 eps) (+ 1.0 eps))))
(FPCore (eps) :precision binary64 (- (- (* eps -2.0) (+ (* 0.6666666666666666 (pow eps 3.0)) (* 0.4 (pow eps 5.0)))) (* 0.2857142857142857 (pow eps 7.0))))
double code(double eps) {
return log((1.0 - eps) / (1.0 + eps));
}
double code(double eps) {
return ((eps * -2.0) - ((0.6666666666666666 * pow(eps, 3.0)) + (0.4 * pow(eps, 5.0)))) - (0.2857142857142857 * pow(eps, 7.0));
}




Bits error versus eps
Results
| Original | 58.6 |
|---|---|
| Target | 0.3 |
| Herbie | 0.2 |
Initial program 58.6
Taylor expanded around 0 0.2
Simplified0.2
rmApplied associate--l-_binary64_10390.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2021007
(FPCore (eps)
:name "logq (problem 3.4.3)"
:precision binary64
:herbie-target
(* -2.0 (+ (+ eps (/ (pow eps 3.0) 3.0)) (/ (pow eps 5.0) 5.0)))
(log (/ (- 1.0 eps) (+ 1.0 eps))))