\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 r129572 = 1.0;
double r129573 = eps;
double r129574 = r129572 - r129573;
double r129575 = r129572 + r129573;
double r129576 = r129574 / r129575;
double r129577 = log(r129576);
return r129577;
}
double f(double eps) {
double r129578 = 2.0;
double r129579 = eps;
double r129580 = 2.0;
double r129581 = pow(r129579, r129580);
double r129582 = 1.0;
double r129583 = pow(r129582, r129580);
double r129584 = r129581 / r129583;
double r129585 = r129584 + r129579;
double r129586 = r129581 - r129585;
double r129587 = r129578 * r129586;
double r129588 = log(r129582);
double r129589 = r129587 + r129588;
return r129589;
}




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