\log \left(\frac{1 - \varepsilon}{1 + \varepsilon}\right)\left({\left(\frac{\varepsilon}{1}\right)}^{3} \cdot \frac{-2}{3} - \frac{2}{5} \cdot \frac{{\varepsilon}^{5}}{{1}^{5}}\right) - 2 \cdot \varepsilondouble f(double eps) {
double r34107 = 1.0;
double r34108 = eps;
double r34109 = r34107 - r34108;
double r34110 = r34107 + r34108;
double r34111 = r34109 / r34110;
double r34112 = log(r34111);
return r34112;
}
double f(double eps) {
double r34113 = eps;
double r34114 = 1.0;
double r34115 = r34113 / r34114;
double r34116 = 3.0;
double r34117 = pow(r34115, r34116);
double r34118 = -0.6666666666666666;
double r34119 = r34117 * r34118;
double r34120 = 0.4;
double r34121 = 5.0;
double r34122 = pow(r34113, r34121);
double r34123 = pow(r34114, r34121);
double r34124 = r34122 / r34123;
double r34125 = r34120 * r34124;
double r34126 = r34119 - r34125;
double r34127 = 2.0;
double r34128 = r34127 * r34113;
double r34129 = r34126 - r34128;
return r34129;
}




Bits error versus eps
Results
| Original | 58.7 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
Initial program 58.7
rmApplied log-div58.6
Simplified58.6
Taylor expanded around 0 0.2
Simplified0.2
rmApplied associate--r+0.2
Final simplification0.2
herbie shell --seed 2019323
(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))))