double code(double eps) {
return ((double) log(((double) (((double) (1.0 - eps)) / ((double) (1.0 + eps))))));
}
double code(double eps) {
return ((double) (((double) log(1.0)) + ((double) (2.0 * ((double) (((double) (eps * ((double) (eps - ((double) (eps / ((double) (1.0 * 1.0)))))))) - eps))))));
}




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