\left(\left(n + 1\right) \cdot \log \left(n + 1\right) - n \cdot \log n\right) - 1
1 \cdot \log n + \left(0.5 \cdot \frac{1}{n} - \frac{0.16666666666666669}{{n}^{2}}\right)double f(double n) {
double r72131 = n;
double r72132 = 1.0;
double r72133 = r72131 + r72132;
double r72134 = log(r72133);
double r72135 = r72133 * r72134;
double r72136 = log(r72131);
double r72137 = r72131 * r72136;
double r72138 = r72135 - r72137;
double r72139 = r72138 - r72132;
return r72139;
}
double f(double n) {
double r72140 = 1.0;
double r72141 = n;
double r72142 = log(r72141);
double r72143 = r72140 * r72142;
double r72144 = 0.5;
double r72145 = 1.0;
double r72146 = r72145 / r72141;
double r72147 = r72144 * r72146;
double r72148 = 0.16666666666666669;
double r72149 = 2.0;
double r72150 = pow(r72141, r72149);
double r72151 = r72148 / r72150;
double r72152 = r72147 - r72151;
double r72153 = r72143 + r72152;
return r72153;
}




Bits error versus n
Results
| Original | 63.0 |
|---|---|
| Target | 0 |
| Herbie | 0 |
Initial program 63.0
Taylor expanded around inf 0.0
Simplified0.0
Taylor expanded around 0 0
Simplified0
Final simplification0
herbie shell --seed 2020035
(FPCore (n)
:name "logs (example 3.8)"
:precision binary64
:pre (> n 6.8e+15)
:herbie-target
(- (log (+ n 1)) (- (/ 1 (* 2 n)) (- (/ 1 (* 3 (* n n))) (/ 4 (pow n 3)))))
(- (- (* (+ n 1) (log (+ n 1))) (* n (log n))) 1))