Average Error: 29.8 → 0.1
Time: 15.1s
Precision: 64
Internal Precision: 128
\[\log \left(N + 1\right) - \log N\]
\[\begin{array}{l} \mathbf{if}\;N \le 10502.73743544022:\\ \;\;\;\;\log \left(\sqrt{\frac{1 + N}{N}}\right) + \log \left(\frac{1 + N}{N}\right) \cdot \frac{1}{2}\\ \mathbf{else}:\\ \;\;\;\;\frac{1}{N} - \frac{1}{N \cdot N} \cdot \left(\frac{1}{2} - \frac{\frac{1}{3}}{N}\right)\\ \end{array}\]

Error

Bits error versus N

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Split input into 2 regimes
  2. if N < 10502.73743544022

    1. Initial program 0.1

      \[\log \left(N + 1\right) - \log N\]
    2. Using strategy rm
    3. Applied diff-log0.1

      \[\leadsto \color{blue}{\log \left(\frac{N + 1}{N}\right)}\]
    4. Using strategy rm
    5. Applied add-sqr-sqrt0.1

      \[\leadsto \log \color{blue}{\left(\sqrt{\frac{N + 1}{N}} \cdot \sqrt{\frac{N + 1}{N}}\right)}\]
    6. Applied log-prod0.1

      \[\leadsto \color{blue}{\log \left(\sqrt{\frac{N + 1}{N}}\right) + \log \left(\sqrt{\frac{N + 1}{N}}\right)}\]
    7. Using strategy rm
    8. Applied pow1/20.1

      \[\leadsto \log \left(\sqrt{\frac{N + 1}{N}}\right) + \log \color{blue}{\left({\left(\frac{N + 1}{N}\right)}^{\frac{1}{2}}\right)}\]
    9. Applied log-pow0.1

      \[\leadsto \log \left(\sqrt{\frac{N + 1}{N}}\right) + \color{blue}{\frac{1}{2} \cdot \log \left(\frac{N + 1}{N}\right)}\]

    if 10502.73743544022 < N

    1. Initial program 59.4

      \[\log \left(N + 1\right) - \log N\]
    2. Using strategy rm
    3. Applied diff-log59.2

      \[\leadsto \color{blue}{\log \left(\frac{N + 1}{N}\right)}\]
    4. 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}}}\]
    5. Simplified0.0

      \[\leadsto \color{blue}{\frac{1}{N} - \frac{1}{N \cdot N} \cdot \left(\frac{1}{2} - \frac{\frac{1}{3}}{N}\right)}\]
  3. Recombined 2 regimes into one program.
  4. Final simplification0.1

    \[\leadsto \begin{array}{l} \mathbf{if}\;N \le 10502.73743544022:\\ \;\;\;\;\log \left(\sqrt{\frac{1 + N}{N}}\right) + \log \left(\frac{1 + N}{N}\right) \cdot \frac{1}{2}\\ \mathbf{else}:\\ \;\;\;\;\frac{1}{N} - \frac{1}{N \cdot N} \cdot \left(\frac{1}{2} - \frac{\frac{1}{3}}{N}\right)\\ \end{array}\]

Reproduce

herbie shell --seed 1042149663 
(FPCore (N)
  :name "2log (problem 3.3.6)"
  (- (log (+ N 1)) (log N)))

Details

Time bar (total: 14.3s)Debug log

sample160.0ms

Algorithm
intervals
Results
472×(pre true 80)
215×(body nan 80)
140×(body real 80)
52×(body real 1280)
33×(body real 640)
26×(body real 320)
(body real 160)

simplify10.0ms

Counts
1 → 1
Calls
1 calls:
Slowest
9.0ms
(- (log (+ N 1)) (log N))

prune11.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 27.2b

localize31.0ms

Local error

Found 2 expressions with local error:

3.1b
(log (+ N 1))
2.2b
(- (log (+ N 1)) (log N))

rewrite14.0ms

Algorithm
rewrite-expression-head
Rules
*-un-lft-identity
log-prod
add-sqr-sqrt
add-cube-cbrt
log-div
pow1
add-exp-log
associate--l+
add-log-exp
distribute-lft-out--
flip-+
log-pow
add-cbrt-cube
flip3-+
associate--l-
difference-of-squares
flip--
rem-log-exp
diff-log
flip3--
sub-neg
Counts
2 → 34
Calls
2 calls:
Slowest
11.0ms
(- (log (+ N 1)) (log N))
2.0ms
(log (+ N 1))

series168.0ms

Counts
2 → 6
Calls
2 calls:
Slowest
106.0ms
(- (log (+ N 1)) (log N))
62.0ms
(log (+ N 1))

simplify1.4s

Counts
22 → 40
Calls
22 calls:
Slowest
380.0ms
(- (+ (/ 1 N) (log -1)) (+ (* 1/2 (/ 1 (pow N 2))) (log (/ -1 N))))
230.0ms
(- (/ 1 N) (+ (* 1/2 (/ 1 (pow N 2))) (log (/ 1 N))))
187.0ms
(- (+ (* 1/3 (/ 1 (pow N 3))) (/ 1 N)) (* 1/2 (/ 1 (pow N 2))))
184.0ms
(- (+ (* 1/3 (/ 1 (pow N 3))) (/ 1 N)) (* 1/2 (/ 1 (pow N 2))))
168.0ms
(- (+ N (* 1/3 (pow N 3))) (* 1/2 (pow N 2)))

prune555.0ms

Pruning

4 alts after pruning (4 fresh and 0 done)

Merged error: 0.1b

localize15.0ms

Local error

Found 2 expressions with local error:

2.1b
(log (/ (+ N 1) N))
0.0b
(/ (+ N 1) N)

rewrite8.0ms

Algorithm
rewrite-expression-head
Rules
add-cube-cbrt
*-un-lft-identity
log-prod
add-sqr-sqrt
associate-/l*
add-exp-log
pow1
associate-/l/
add-log-exp
div-inv
add-cbrt-cube
flip-+
log-pow
rem-log-exp
flip3-+
frac-2neg
clear-num
log-div
Counts
2 → 29
Calls
2 calls:
Slowest
5.0ms
(/ (+ N 1) N)
2.0ms
(log (/ (+ N 1) N))

series106.0ms

Counts
2 → 6
Calls
2 calls:
Slowest
66.0ms
(log (/ (+ N 1) N))
40.0ms
(/ (+ N 1) N)

simplify564.0ms

Counts
16 → 35
Calls
16 calls:
Slowest
199.0ms
(- (+ (* 1/3 (/ 1 (pow N 3))) (/ 1 N)) (* 1/2 (/ 1 (pow N 2))))
196.0ms
(- (+ (* 1/3 (/ 1 (pow N 3))) (/ 1 N)) (* 1/2 (/ 1 (pow N 2))))
95.0ms
(- N (+ (log N) (* 1/2 (pow N 2))))
17.0ms
(- (* N N) (* 1 1))
7.0ms
(+ (pow N 3) (pow 1 3))

prune446.0ms

Pruning

5 alts after pruning (4 fresh and 1 done)

Merged error: 0.1b

localize27.0ms

Local error

Found 4 expressions with local error:

2.1b
(log (sqrt (/ (+ N 1) N)))
2.1b
(log (sqrt (/ (+ N 1) N)))
0.1b
(sqrt (/ (+ N 1) N))
0.1b
(sqrt (/ (+ N 1) N))

rewrite10.0ms

Algorithm
rewrite-expression-head
Rules
16×sqrt-prod
14×log-prod
12×add-sqr-sqrt
10×add-cube-cbrt
10×*-un-lft-identity
add-exp-log
pow1
add-log-exp
div-inv
log-pow
sqrt-div
pow1/2
add-cbrt-cube
rem-log-exp
log-div
rem-sqrt-square
Counts
4 → 64
Calls
4 calls:
Slowest
3.0ms
(log (sqrt (/ (+ N 1) N)))
3.0ms
(log (sqrt (/ (+ N 1) N)))
2.0ms
(sqrt (/ (+ N 1) N))
2.0ms
(sqrt (/ (+ N 1) N))

series250.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
87.0ms
(log (sqrt (/ (+ N 1) N)))
86.0ms
(log (sqrt (/ (+ N 1) N)))
39.0ms
(sqrt (/ (+ N 1) N))
37.0ms
(sqrt (/ (+ N 1) N))

simplify1.5s

Counts
42 → 76
Calls
42 calls:
Slowest
217.0ms
(- (+ (* 1/6 (/ 1 (pow N 3))) (* 1/2 (/ 1 N))) (* 1/4 (/ 1 (pow N 2))))
215.0ms
(- (+ (* 1/6 (/ 1 (pow N 3))) (* 1/2 (/ 1 N))) (* 1/4 (/ 1 (pow N 2))))
213.0ms
(- (+ (* 1/6 (/ 1 (pow N 3))) (* 1/2 (/ 1 N))) (* 1/4 (/ 1 (pow N 2))))
207.0ms
(- (+ (* 1/6 (/ 1 (pow N 3))) (* 1/2 (/ 1 N))) (* 1/4 (/ 1 (pow N 2))))
139.0ms
(- (+ (* 1/2 (/ 1 N)) 1) (* 1/8 (/ 1 (pow N 2))))

prune1.0s

Pruning

4 alts after pruning (3 fresh and 1 done)

Merged error: 0.1b

localize16.0ms

Local error

Found 4 expressions with local error:

2.1b
(log (sqrt (/ (+ N 1) N)))
2.1b
(log (/ (+ N 1) N))
0.1b
(sqrt (/ (+ N 1) N))
0.0b
(/ (+ N 1) N)

rewrite12.0ms

Algorithm
rewrite-expression-head
Rules
11×log-prod
10×add-sqr-sqrt
add-cube-cbrt
*-un-lft-identity
sqrt-prod
add-exp-log
pow1
add-log-exp
div-inv
add-cbrt-cube
associate-/l*
log-pow
associate-/l/
rem-log-exp
sqrt-div
pow1/2
log-div
flip-+
flip3-+
frac-2neg
clear-num
rem-sqrt-square
Counts
4 → 61
Calls
4 calls:
Slowest
5.0ms
(/ (+ N 1) N)
2.0ms
(log (sqrt (/ (+ N 1) N)))
2.0ms
(log (/ (+ N 1) N))
2.0ms
(sqrt (/ (+ N 1) N))

series219.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
80.0ms
(log (sqrt (/ (+ N 1) N)))
69.0ms
(log (/ (+ N 1) N))
39.0ms
(/ (+ N 1) N)
31.0ms
(sqrt (/ (+ N 1) N))

simplify1.3s

Counts
37 → 73
Calls
37 calls:
Slowest
214.0ms
(- (+ (* 1/6 (/ 1 (pow N 3))) (* 1/2 (/ 1 N))) (* 1/4 (/ 1 (pow N 2))))
208.0ms
(- (+ (* 1/6 (/ 1 (pow N 3))) (* 1/2 (/ 1 N))) (* 1/4 (/ 1 (pow N 2))))
189.0ms
(- (+ (* 1/3 (/ 1 (pow N 3))) (/ 1 N)) (* 1/2 (/ 1 (pow N 2))))
183.0ms
(- (+ (* 1/3 (/ 1 (pow N 3))) (/ 1 N)) (* 1/2 (/ 1 (pow N 2))))
141.0ms
(- (+ (* 1/2 (/ 1 N)) 1) (* 1/8 (/ 1 (pow N 2))))

prune1.0s

Pruning

4 alts after pruning (2 fresh and 2 done)

Merged error: 0.1b

regimes60.0ms

Accuracy

100% (0.0b remaining)

Error of 0.1b against oracle of 0.1b and baseline of 29.6b

bsearch78.0ms

end0.0ms

sample5.2s

Algorithm
intervals
Results
16095×(pre true 80)
8094×(body nan 80)
4034×(body real 80)
1748×(body real 1280)
1291×(body real 640)
638×(body real 320)
290×(body real 160)