\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 r90041 = 1.0;
double r90042 = eps;
double r90043 = r90041 - r90042;
double r90044 = r90041 + r90042;
double r90045 = r90043 / r90044;
double r90046 = log(r90045);
return r90046;
}
double f(double eps) {
double r90047 = 2.0;
double r90048 = eps;
double r90049 = 2.0;
double r90050 = pow(r90048, r90049);
double r90051 = 1.0;
double r90052 = pow(r90051, r90049);
double r90053 = r90050 / r90052;
double r90054 = r90053 + r90048;
double r90055 = r90050 - r90054;
double r90056 = r90047 * r90055;
double r90057 = log(r90051);
double r90058 = r90056 + r90057;
return r90058;
}




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