\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 r93624 = n;
double r93625 = 1.0;
double r93626 = r93624 + r93625;
double r93627 = log(r93626);
double r93628 = r93626 * r93627;
double r93629 = log(r93624);
double r93630 = r93624 * r93629;
double r93631 = r93628 - r93630;
double r93632 = r93631 - r93625;
return r93632;
}
double f(double n) {
double r93633 = 1.0;
double r93634 = n;
double r93635 = log(r93634);
double r93636 = r93633 * r93635;
double r93637 = 0.5;
double r93638 = 1.0;
double r93639 = r93638 / r93634;
double r93640 = r93637 * r93639;
double r93641 = 0.16666666666666669;
double r93642 = 2.0;
double r93643 = pow(r93634, r93642);
double r93644 = r93641 / r93643;
double r93645 = r93640 - r93644;
double r93646 = r93636 + r93645;
return r93646;
}




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 2020064
(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))