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








Bits error versus eps
Results
| Original | 58.6 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
| Alternative 1 | |
|---|---|
| Error | 0.3 |
| Cost | 6912 |
| Alternative 2 | |
|---|---|
| Error | 0.7 |
| Cost | 192 |
| Alternative 3 | |
|---|---|
| Error | 60.6 |
| Cost | 64 |
| Alternative 4 | |
|---|---|
| Error | 61.6 |
| Cost | 64 |

Initial program 58.6
Taylor expanded around 0 0.2
Simplified0.2
rmApplied unpow3_binary64_15080.2
Applied associate-*l*_binary64_13830.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2021044
(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))))