\log \left(\frac{1 - \varepsilon}{1 + \varepsilon}\right){\varepsilon}^{5} \cdot \frac{-2}{5} - \frac{\left(\left(\left(\varepsilon \cdot \varepsilon\right) \cdot \frac{8}{27}\right) \cdot \left(\left(\varepsilon \cdot \varepsilon\right) \cdot \left(\varepsilon \cdot \varepsilon\right)\right) + 8\right) \cdot \varepsilon}{\left(4 - \left(\varepsilon \cdot \left(\varepsilon \cdot \frac{2}{3}\right)\right) \cdot 2\right) + \left(\varepsilon \cdot \left(\varepsilon \cdot \frac{2}{3}\right)\right) \cdot \left(\varepsilon \cdot \left(\varepsilon \cdot \frac{2}{3}\right)\right)}double f(double eps) {
double r5290553 = 1.0;
double r5290554 = eps;
double r5290555 = r5290553 - r5290554;
double r5290556 = r5290553 + r5290554;
double r5290557 = r5290555 / r5290556;
double r5290558 = log(r5290557);
return r5290558;
}
double f(double eps) {
double r5290559 = eps;
double r5290560 = 5.0;
double r5290561 = pow(r5290559, r5290560);
double r5290562 = -0.4;
double r5290563 = r5290561 * r5290562;
double r5290564 = r5290559 * r5290559;
double r5290565 = 0.2962962962962963;
double r5290566 = r5290564 * r5290565;
double r5290567 = r5290564 * r5290564;
double r5290568 = r5290566 * r5290567;
double r5290569 = 8.0;
double r5290570 = r5290568 + r5290569;
double r5290571 = r5290570 * r5290559;
double r5290572 = 4.0;
double r5290573 = 0.6666666666666666;
double r5290574 = r5290559 * r5290573;
double r5290575 = r5290559 * r5290574;
double r5290576 = 2.0;
double r5290577 = r5290575 * r5290576;
double r5290578 = r5290572 - r5290577;
double r5290579 = r5290575 * r5290575;
double r5290580 = r5290578 + r5290579;
double r5290581 = r5290571 / r5290580;
double r5290582 = r5290563 - r5290581;
return r5290582;
}




Bits error versus eps
Results
| Original | 58.6 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
Initial program 58.6
Taylor expanded around 0 0.2
Simplified0.2
rmApplied flip3-+0.2
Applied associate-*l/0.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2019104
(FPCore (eps)
:name "logq (problem 3.4.3)"
:herbie-target
(* -2 (+ (+ eps (/ (pow eps 3) 3)) (/ (pow eps 5) 5)))
(log (/ (- 1 eps) (+ 1 eps))))