Average Error: 29.6 → 0.6
Time: 31.0s
Precision: 64
Internal Precision: 1408
\[\log \left(N + 1\right) - \log N\]
\[\begin{array}{l} \mathbf{if}\;\log \left(N + 1\right) - \log N \le 9.045593907618368 \cdot 10^{-32}:\\ \;\;\;\;(\left(\frac{1}{N \cdot N}\right) \cdot \left(\frac{\frac{1}{3}}{N} - \frac{1}{2}\right) + \left(\frac{1}{N}\right))_*\\ \mathbf{else}:\\ \;\;\;\;\log \left(\frac{e^{\log_* (1 + N)}}{N}\right)\\ \end{array}\]

Error

Bits error versus N

Derivation

  1. Split input into 2 regimes
  2. if (- (log (+ N 1)) (log N)) < 9.045593907618368e-32

    1. Initial program 60.5

      \[\log \left(N + 1\right) - \log N\]
    2. Applied simplify60.5

      \[\leadsto \color{blue}{\log_* (1 + N) - \log N}\]
    3. Taylor expanded around inf 0.0

      \[\leadsto \color{blue}{\left(\frac{1}{3} \cdot \frac{1}{{N}^{3}} + \frac{1}{N}\right) - \frac{1}{2} \cdot \frac{1}{{N}^{2}}}\]
    4. Applied simplify0.0

      \[\leadsto \color{blue}{(\left(\frac{1}{N \cdot N}\right) \cdot \left(\frac{\frac{1}{3}}{N} - \frac{1}{2}\right) + \left(\frac{1}{N}\right))_*}\]

    if 9.045593907618368e-32 < (- (log (+ N 1)) (log N))

    1. Initial program 1.2

      \[\log \left(N + 1\right) - \log N\]
    2. Applied simplify1.2

      \[\leadsto \color{blue}{\log_* (1 + N) - \log N}\]
    3. Using strategy rm
    4. Applied add-log-exp1.2

      \[\leadsto \color{blue}{\log \left(e^{\log_* (1 + N)}\right)} - \log N\]
    5. Applied diff-log1.2

      \[\leadsto \color{blue}{\log \left(\frac{e^{\log_* (1 + N)}}{N}\right)}\]
  3. Recombined 2 regimes into one program.

Runtime

Time bar (total: 31.0s)Debug logProfile

herbie shell --seed '#(1070833653 108281690 3330367898 3632331308 3494323072 43156186)' +o rules:numerics
(FPCore (N)
  :name "2log (problem 3.3.6)"
  (- (log (+ N 1)) (log N)))