\log \left(\frac{1 - \varepsilon}{1 + \varepsilon}\right)2 \cdot \left(\varepsilon \cdot \varepsilon - \left(\frac{{\varepsilon}^{2}}{{1}^{2}} + \varepsilon\right)\right) + \log 1double f(double eps) {
double r114408 = 1.0;
double r114409 = eps;
double r114410 = r114408 - r114409;
double r114411 = r114408 + r114409;
double r114412 = r114410 / r114411;
double r114413 = log(r114412);
return r114413;
}
double f(double eps) {
double r114414 = 2.0;
double r114415 = eps;
double r114416 = r114415 * r114415;
double r114417 = 2.0;
double r114418 = pow(r114415, r114417);
double r114419 = 1.0;
double r114420 = pow(r114419, r114417);
double r114421 = r114418 / r114420;
double r114422 = r114421 + r114415;
double r114423 = r114416 - r114422;
double r114424 = r114414 * r114423;
double r114425 = log(r114419);
double r114426 = r114424 + r114425;
return r114426;
}




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