\log \left(\frac{1 - \varepsilon}{1 + \varepsilon}\right)2 \cdot \left({\varepsilon}^{2} - \left(\frac{{\varepsilon}^{2}}{{1}^{2}} + \varepsilon\right)\right) + \log 1double f(double eps) {
double r84285 = 1.0;
double r84286 = eps;
double r84287 = r84285 - r84286;
double r84288 = r84285 + r84286;
double r84289 = r84287 / r84288;
double r84290 = log(r84289);
return r84290;
}
double f(double eps) {
double r84291 = 2.0;
double r84292 = eps;
double r84293 = 2.0;
double r84294 = pow(r84292, r84293);
double r84295 = 1.0;
double r84296 = pow(r84295, r84293);
double r84297 = r84294 / r84296;
double r84298 = r84297 + r84292;
double r84299 = r84294 - r84298;
double r84300 = r84291 * r84299;
double r84301 = log(r84295);
double r84302 = r84300 + r84301;
return r84302;
}




Bits error versus eps
Results
| Original | 58.5 |
|---|---|
| Target | 0.2 |
| Herbie | 0.7 |
Initial program 58.5
Taylor expanded around 0 0.7
Simplified0.7
Final simplification0.7
herbie shell --seed 2020021
(FPCore (eps)
:name "logq (problem 3.4.3)"
:precision binary64
:herbie-target
(* -2 (+ (+ eps (/ (pow eps 3) 3)) (/ (pow eps 5) 5)))
(log (/ (- 1 eps) (+ 1 eps))))