\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 r78780 = 1.0;
double r78781 = eps;
double r78782 = r78780 - r78781;
double r78783 = r78780 + r78781;
double r78784 = r78782 / r78783;
double r78785 = log(r78784);
return r78785;
}
double f(double eps) {
double r78786 = 2.0;
double r78787 = eps;
double r78788 = 2.0;
double r78789 = pow(r78787, r78788);
double r78790 = 1.0;
double r78791 = pow(r78790, r78788);
double r78792 = r78789 / r78791;
double r78793 = r78792 + r78787;
double r78794 = r78789 - r78793;
double r78795 = r78786 * r78794;
double r78796 = log(r78790);
double r78797 = r78795 + r78796;
return r78797;
}




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