\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 r85837 = 1.0;
double r85838 = eps;
double r85839 = r85837 - r85838;
double r85840 = r85837 + r85838;
double r85841 = r85839 / r85840;
double r85842 = log(r85841);
return r85842;
}
double f(double eps) {
double r85843 = 2.0;
double r85844 = eps;
double r85845 = r85844 * r85844;
double r85846 = 1.0;
double r85847 = r85844 / r85846;
double r85848 = r85847 * r85847;
double r85849 = r85844 + r85848;
double r85850 = r85845 - r85849;
double r85851 = r85843 * r85850;
double r85852 = log(r85846);
double r85853 = r85851 + r85852;
return r85853;
}




Bits error versus eps
Results
| Original | 58.4 |
|---|---|
| Target | 0.2 |
| Herbie | 0.7 |
Initial program 58.4
Simplified58.4
Taylor expanded around 0 0.7
Simplified0.7
Final simplification0.7
herbie shell --seed 2019174
(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))))