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




Bits error versus eps
Results
| Original | 58.6 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
Initial program 58.6
Taylor expanded around 0 0.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2020219
(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))))