\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 r86614 = 1.0;
double r86615 = eps;
double r86616 = r86614 - r86615;
double r86617 = r86614 + r86615;
double r86618 = r86616 / r86617;
double r86619 = log(r86618);
return r86619;
}
double f(double eps) {
double r86620 = 2.0;
double r86621 = eps;
double r86622 = 2.0;
double r86623 = pow(r86621, r86622);
double r86624 = 1.0;
double r86625 = pow(r86624, r86622);
double r86626 = r86623 / r86625;
double r86627 = r86626 + r86621;
double r86628 = r86623 - r86627;
double r86629 = r86620 * r86628;
double r86630 = log(r86624);
double r86631 = r86629 + r86630;
return r86631;
}




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