\left(\left(n + 1\right) \cdot \log \left(n + 1\right) - n \cdot \log n\right) - 1
\left(\frac{0.5}{n} - \frac{0.1666666666666666851703837437526090070605}{n \cdot n}\right) + \log n \cdot 1double f(double n) {
double r88397 = n;
double r88398 = 1.0;
double r88399 = r88397 + r88398;
double r88400 = log(r88399);
double r88401 = r88399 * r88400;
double r88402 = log(r88397);
double r88403 = r88397 * r88402;
double r88404 = r88401 - r88403;
double r88405 = r88404 - r88398;
return r88405;
}
double f(double n) {
double r88406 = 0.5;
double r88407 = n;
double r88408 = r88406 / r88407;
double r88409 = 0.16666666666666669;
double r88410 = r88407 * r88407;
double r88411 = r88409 / r88410;
double r88412 = r88408 - r88411;
double r88413 = log(r88407);
double r88414 = 1.0;
double r88415 = r88413 * r88414;
double r88416 = r88412 + r88415;
return r88416;
}




Bits error versus n
Results
| 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 2019194
(FPCore (n)
:name "logs (example 3.8)"
:pre (> n 6.8e+15)
:herbie-target
(- (log (+ n 1.0)) (- (/ 1.0 (* 2.0 n)) (- (/ 1.0 (* 3.0 (* n n))) (/ 4.0 (pow n 3.0)))))
(- (- (* (+ n 1.0) (log (+ n 1.0))) (* n (log n))) 1.0))