\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 r90443 = 1.0;
double r90444 = eps;
double r90445 = r90443 - r90444;
double r90446 = r90443 + r90444;
double r90447 = r90445 / r90446;
double r90448 = log(r90447);
return r90448;
}
double f(double eps) {
double r90449 = 2.0;
double r90450 = eps;
double r90451 = 2.0;
double r90452 = pow(r90450, r90451);
double r90453 = 1.0;
double r90454 = pow(r90453, r90451);
double r90455 = r90452 / r90454;
double r90456 = r90455 + r90450;
double r90457 = r90452 - r90456;
double r90458 = r90449 * r90457;
double r90459 = log(r90453);
double r90460 = r90458 + r90459;
return r90460;
}




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 2020060
(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))))