\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 r92543 = 1.0;
double r92544 = eps;
double r92545 = r92543 - r92544;
double r92546 = r92543 + r92544;
double r92547 = r92545 / r92546;
double r92548 = log(r92547);
return r92548;
}
double f(double eps) {
double r92549 = 2.0;
double r92550 = eps;
double r92551 = 2.0;
double r92552 = pow(r92550, r92551);
double r92553 = 1.0;
double r92554 = pow(r92553, r92551);
double r92555 = r92552 / r92554;
double r92556 = r92555 + r92550;
double r92557 = r92552 - r92556;
double r92558 = r92549 * r92557;
double r92559 = log(r92553);
double r92560 = r92558 + r92559;
return r92560;
}




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