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




Bits error versus eps
Results
| Original | 58.5 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
Initial program 58.5
rmApplied log-div58.5
Taylor expanded around 0 0.2
Simplified0.2
rmApplied add-log-exp0.3
Simplified0.2
Final simplification0.2
herbie shell --seed 2020196
(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))))