\log \left(\frac{1 - \varepsilon}{1 + \varepsilon}\right)2 \cdot \left(\varepsilon \cdot \varepsilon - \left(\varepsilon + \frac{\varepsilon}{1} \cdot \frac{\varepsilon}{1}\right)\right) + \log 1double f(double eps) {
double r3563622 = 1.0;
double r3563623 = eps;
double r3563624 = r3563622 - r3563623;
double r3563625 = r3563622 + r3563623;
double r3563626 = r3563624 / r3563625;
double r3563627 = log(r3563626);
return r3563627;
}
double f(double eps) {
double r3563628 = 2.0;
double r3563629 = eps;
double r3563630 = r3563629 * r3563629;
double r3563631 = 1.0;
double r3563632 = r3563629 / r3563631;
double r3563633 = r3563632 * r3563632;
double r3563634 = r3563629 + r3563633;
double r3563635 = r3563630 - r3563634;
double r3563636 = r3563628 * r3563635;
double r3563637 = log(r3563631);
double r3563638 = r3563636 + r3563637;
return r3563638;
}




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 2019172
(FPCore (eps)
:name "logq (problem 3.4.3)"
:herbie-target
(* -2.0 (+ (+ eps (/ (pow eps 3.0) 3.0)) (/ (pow eps 5.0) 5.0)))
(log (/ (- 1.0 eps) (+ 1.0 eps))))