\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 r101145 = 1.0;
double r101146 = eps;
double r101147 = r101145 - r101146;
double r101148 = r101145 + r101146;
double r101149 = r101147 / r101148;
double r101150 = log(r101149);
return r101150;
}
double f(double eps) {
double r101151 = 2.0;
double r101152 = eps;
double r101153 = 2.0;
double r101154 = pow(r101152, r101153);
double r101155 = 1.0;
double r101156 = pow(r101155, r101153);
double r101157 = r101154 / r101156;
double r101158 = r101157 + r101152;
double r101159 = r101154 - r101158;
double r101160 = r101151 * r101159;
double r101161 = log(r101155);
double r101162 = r101160 + r101161;
return r101162;
}




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