\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 r151044 = 1.0;
double r151045 = eps;
double r151046 = r151044 - r151045;
double r151047 = r151044 + r151045;
double r151048 = r151046 / r151047;
double r151049 = log(r151048);
return r151049;
}
double f(double eps) {
double r151050 = 2.0;
double r151051 = eps;
double r151052 = 2.0;
double r151053 = pow(r151051, r151052);
double r151054 = 1.0;
double r151055 = pow(r151054, r151052);
double r151056 = r151053 / r151055;
double r151057 = r151056 + r151051;
double r151058 = r151053 - r151057;
double r151059 = r151050 * r151058;
double r151060 = log(r151054);
double r151061 = r151059 + r151060;
return r151061;
}




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