\log \left(\frac{1 - \varepsilon}{1 + \varepsilon}\right)2 \cdot \left({\varepsilon}^{2} - \left(\frac{{\varepsilon}^{2}}{{1}^{2}} + \varepsilon\right)\right) + \log 1double f(double eps) {
double r81982 = 1.0;
double r81983 = eps;
double r81984 = r81982 - r81983;
double r81985 = r81982 + r81983;
double r81986 = r81984 / r81985;
double r81987 = log(r81986);
return r81987;
}
double f(double eps) {
double r81988 = 2.0;
double r81989 = eps;
double r81990 = 2.0;
double r81991 = pow(r81989, r81990);
double r81992 = 1.0;
double r81993 = pow(r81992, r81990);
double r81994 = r81991 / r81993;
double r81995 = r81994 + r81989;
double r81996 = r81991 - r81995;
double r81997 = r81988 * r81996;
double r81998 = log(r81992);
double r81999 = r81997 + r81998;
return r81999;
}




Bits error versus eps
Results
| Original | 58.7 |
|---|---|
| Target | 0.2 |
| Herbie | 0.6 |
Initial program 58.7
Taylor expanded around 0 0.6
Simplified0.6
Final simplification0.6
herbie shell --seed 2020035
(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))))