\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 r92966 = 1.0;
double r92967 = eps;
double r92968 = r92966 - r92967;
double r92969 = r92966 + r92967;
double r92970 = r92968 / r92969;
double r92971 = log(r92970);
return r92971;
}
double f(double eps) {
double r92972 = 2.0;
double r92973 = eps;
double r92974 = 2.0;
double r92975 = pow(r92973, r92974);
double r92976 = 1.0;
double r92977 = pow(r92976, r92974);
double r92978 = r92975 / r92977;
double r92979 = r92978 + r92973;
double r92980 = r92975 - r92979;
double r92981 = r92972 * r92980;
double r92982 = log(r92976);
double r92983 = r92981 + r92982;
return r92983;
}




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