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




Bits error versus eps
Results
| Original | 58.7 |
|---|---|
| Target | 0.2 |
| Herbie | 0.1 |
Initial program 58.7
Taylor expanded around 0 0.1
Simplified0.1
rmApplied associate--l-_binary640.1
Simplified0.1
Final simplification0.1
herbie shell --seed 2021205
(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))))