double f(double n) {
double r8618005 = n;
double r8618006 = 1.0;
double r8618007 = r8618005 + r8618006;
double r8618008 = log(r8618007);
double r8618009 = r8618007 * r8618008;
double r8618010 = log(r8618005);
double r8618011 = r8618005 * r8618010;
double r8618012 = r8618009 - r8618011;
double r8618013 = r8618012 - r8618006;
return r8618013;
}
double f(double n) {
double r8618014 = 0.5;
double r8618015 = n;
double r8618016 = r8618014 / r8618015;
double r8618017 = 0.16666666666666666;
double r8618018 = r8618015 * r8618015;
double r8618019 = r8618017 / r8618018;
double r8618020 = log(r8618015);
double r8618021 = r8618019 - r8618020;
double r8618022 = r8618016 - r8618021;
return r8618022;
}
\left(\left(n + 1\right) \cdot \log \left(n + 1\right) - n \cdot \log n\right) - 1
\frac{\frac{1}{2}}{n} - \left(\frac{\frac{1}{6}}{n \cdot n} - \log n\right)



Bits error versus n
| Original | 63.0 |
|---|---|
| Target | 0 |
| Herbie | 0 |
Initial program 63.0
Simplified62.0
Taylor expanded around inf 0.0
Simplified0
Final simplification0
herbie shell --seed 2019102 +o rules:numerics
(FPCore (n)
:name "logs (example 3.8)"
: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))