\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 r132534 = 1.0;
double r132535 = eps;
double r132536 = r132534 - r132535;
double r132537 = r132534 + r132535;
double r132538 = r132536 / r132537;
double r132539 = log(r132538);
return r132539;
}
double f(double eps) {
double r132540 = 2.0;
double r132541 = eps;
double r132542 = 2.0;
double r132543 = pow(r132541, r132542);
double r132544 = 1.0;
double r132545 = pow(r132544, r132542);
double r132546 = r132543 / r132545;
double r132547 = r132546 + r132541;
double r132548 = r132543 - r132547;
double r132549 = r132540 * r132548;
double r132550 = log(r132544);
double r132551 = r132549 + r132550;
return r132551;
}




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