\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 r92153 = 1.0;
double r92154 = eps;
double r92155 = r92153 - r92154;
double r92156 = r92153 + r92154;
double r92157 = r92155 / r92156;
double r92158 = log(r92157);
return r92158;
}
double f(double eps) {
double r92159 = 2.0;
double r92160 = eps;
double r92161 = 2.0;
double r92162 = pow(r92160, r92161);
double r92163 = 1.0;
double r92164 = pow(r92163, r92161);
double r92165 = r92162 / r92164;
double r92166 = r92165 + r92160;
double r92167 = r92162 - r92166;
double r92168 = r92159 * r92167;
double r92169 = log(r92163);
double r92170 = r92168 + r92169;
return r92170;
}




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