\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 r3660092 = 1.0;
double r3660093 = eps;
double r3660094 = r3660092 - r3660093;
double r3660095 = r3660092 + r3660093;
double r3660096 = r3660094 / r3660095;
double r3660097 = log(r3660096);
return r3660097;
}
double f(double eps) {
double r3660098 = 2.0;
double r3660099 = eps;
double r3660100 = r3660099 * r3660099;
double r3660101 = 1.0;
double r3660102 = r3660099 / r3660101;
double r3660103 = r3660102 * r3660102;
double r3660104 = r3660099 + r3660103;
double r3660105 = r3660100 - r3660104;
double r3660106 = r3660098 * r3660105;
double r3660107 = log(r3660101);
double r3660108 = r3660106 + r3660107;
return r3660108;
}




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