\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 r59916 = 1.0;
double r59917 = eps;
double r59918 = r59916 - r59917;
double r59919 = r59916 + r59917;
double r59920 = r59918 / r59919;
double r59921 = log(r59920);
return r59921;
}
double f(double eps) {
double r59922 = 2.0;
double r59923 = eps;
double r59924 = 2.0;
double r59925 = pow(r59923, r59924);
double r59926 = 1.0;
double r59927 = pow(r59926, r59924);
double r59928 = r59925 / r59927;
double r59929 = r59928 + r59923;
double r59930 = r59925 - r59929;
double r59931 = r59922 * r59930;
double r59932 = log(r59926);
double r59933 = r59931 + r59932;
return r59933;
}




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