\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 r97171 = 1.0;
double r97172 = eps;
double r97173 = r97171 - r97172;
double r97174 = r97171 + r97172;
double r97175 = r97173 / r97174;
double r97176 = log(r97175);
return r97176;
}
double f(double eps) {
double r97177 = 2.0;
double r97178 = eps;
double r97179 = 2.0;
double r97180 = pow(r97178, r97179);
double r97181 = 1.0;
double r97182 = pow(r97181, r97179);
double r97183 = r97180 / r97182;
double r97184 = r97183 + r97178;
double r97185 = r97180 - r97184;
double r97186 = r97177 * r97185;
double r97187 = log(r97181);
double r97188 = r97186 + r97187;
return r97188;
}




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