Average Error: 0.4 → 0.5
Time: 30.0s
Precision: 64
Internal Precision: 128
\[\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}\]
↓
\[\frac{\sqrt{{\left(\left(2 \cdot n\right) \cdot \pi\right)}^{\left(\frac{1}{2} - \frac{k}{2}\right)}} \cdot \sqrt{{\left(\left(2 \cdot n\right) \cdot \pi\right)}^{\left(\frac{1}{2} - \frac{k}{2}\right)}}}{\sqrt{k}}\]

Error

Bits error versus k

Bits error versus n

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Initial program 0.4

    \[\frac{1}{\sqrt{k}} \cdot {\left(\left(2 \cdot \pi\right) \cdot n\right)}^{\left(\frac{1 - k}{2}\right)}\]

  2. Simplified0.4

    \[\leadsto \color{blue}{\frac{{\left(\pi \cdot \left(n \cdot 2\right)\right)}^{\left(\frac{1}{2} - \frac{k}{2}\right)}}{\sqrt{k}}}\]
  3. Using strategy rm

  4. Applied add-sqr-sqrt0.5

    \[\leadsto \frac{\color{blue}{\sqrt{{\left(\pi \cdot \left(n \cdot 2\right)\right)}^{\left(\frac{1}{2} - \frac{k}{2}\right)}} \cdot \sqrt{{\left(\pi \cdot \left(n \cdot 2\right)\right)}^{\left(\frac{1}{2} - \frac{k}{2}\right)}}}}{\sqrt{k}}\]

  5. Final simplification0.5

    \[\leadsto \frac{\sqrt{{\left(\left(2 \cdot n\right) \cdot \pi\right)}^{\left(\frac{1}{2} - \frac{k}{2}\right)}} \cdot \sqrt{{\left(\left(2 \cdot n\right) \cdot \pi\right)}^{\left(\frac{1}{2} - \frac{k}{2}\right)}}}{\sqrt{k}}\]

Reproduce

herbie shell --seed 2019008 
(FPCore (k n)
  :name "Migdal et al, Equation (51)"
  (* (/ 1 (sqrt k)) (pow (* (* 2 PI) n) (/ (- 1 k) 2))))

Details

Time bar (total: 29.2s)Debug log

sample326.0ms

Algorithm
intervals

simplify73.0ms

Counts
1 → 1
Calls

1 calls. Slowest were:

72.0ms
(* (/ 1 (sqrt k)) (pow (* (* 2 PI) n) (/ (- 1 k) 2)))

prune15.0ms

Pruning

2 alts after pruning (2 fresh and 0 done)

Merged error: 0.3b

localize22.0ms

Local error

Found 3 expressions with local error:

0.5b
(pow (* PI (* n 2)) (- 1/2 (/ k 2)))
0.2b
(* PI (* n 2))
0.1b
(/ (pow (* PI (* n 2)) (- 1/2 (/ k 2))) (sqrt k))

rewrite9.0ms

Algorithm
rewrite-expression-head
Counts
3 → 68
Calls

3 calls. Slowest were:

5.0ms
(/ (pow (* PI (* n 2)) (- 1/2 (/ k 2))) (sqrt k))
2.0ms
(pow (* PI (* n 2)) (- 1/2 (/ k 2)))
1.0ms
(* PI (* n 2))

series384.0ms

Counts
3 → 9
Calls

3 calls. Slowest were:

202.0ms
(/ (pow (* PI (* n 2)) (- 1/2 (/ k 2))) (sqrt k))
124.0ms
(pow (* PI (* n 2)) (- 1/2 (/ k 2)))
58.0ms
(* PI (* n 2))

simplify3.2s

Counts
67 → 77
Calls

67 calls. Slowest were:

607.0ms
(- (* (log (* PI (* n 2))) (- 1/2 (/ k 2))) (log (sqrt k)))
566.0ms
(- (* (log (* PI (* n 2))) (- 1/2 (/ k 2))) (log (sqrt k)))
260.0ms
(/ (* (* (pow (* PI (* n 2)) (- 1/2 (/ k 2))) (pow (* PI (* n 2)) (- 1/2 (/ k 2)))) (pow (* PI (* n 2)) (- 1/2 (/ k 2)))) (* (* (sqrt k) (sqrt k)) (sqrt k)))

prune787.0ms

Pruning

7 alts after pruning (7 fresh and 0 done)

Merged error: 0.0b

localize16.0ms

Local error

Found 4 expressions with local error:

0.5b
(pow (* PI (* n 2)) (- 1/2 (/ k 2)))
0.5b
(pow (* PI (* n 2)) (- 1/2 (/ k 2)))
0.2b
(* (sqrt (pow (* PI (* n 2)) (- 1/2 (/ k 2)))) (sqrt (pow (* PI (* n 2)) (- 1/2 (/ k 2)))))
0.2b
(* PI (* n 2))

rewrite21.0ms

Algorithm
rewrite-expression-head
Counts
4 → 82
Calls

4 calls. Slowest were:

13.0ms
(* (sqrt (pow (* PI (* n 2)) (- 1/2 (/ k 2)))) (sqrt (pow (* PI (* n 2)) (- 1/2 (/ k 2)))))
3.0ms
(pow (* PI (* n 2)) (- 1/2 (/ k 2)))
3.0ms
(pow (* PI (* n 2)) (- 1/2 (/ k 2)))

series451.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

132.0ms
(* (sqrt (pow (* PI (* n 2)) (- 1/2 (/ k 2)))) (sqrt (pow (* PI (* n 2)) (- 1/2 (/ k 2)))))
125.0ms
(pow (* PI (* n 2)) (- 1/2 (/ k 2)))
112.0ms
(pow (* PI (* n 2)) (- 1/2 (/ k 2)))
82.0ms
(* PI (* n 2))

simplify2.7s

Counts
63 → 94
Calls

63 calls. Slowest were:

294.0ms
(- (+ (* 1/4 (* (log (* 2 PI)) (* (exp (* 1/2 (+ (log n) (log (* 2 PI))))) (* (log n) (pow k 2))))) (+ (* 1/8 (* (exp (* 1/2 (+ (log n) (log (* 2 PI))))) (* (pow (log n) 2) (pow k 2)))) (+ (exp (* 1/2 (+ (log n) (log (* 2 PI))))) (* 1/8 (* (pow (log (* 2 PI)) 2) (* (exp (* 1/2 (+ (log n) (log (* 2 PI))))) (pow k 2))))))) (+ (* 1/2 (* (exp (* 1/2 (+ (log n) (log (* 2 PI))))) (* (log n) k))) (* 1/2 (* (log (* 2 PI)) (* (exp (* 1/2 (+ (log n) (log (* 2 PI))))) k)))))
276.0ms
(* (log (* PI (* n 2))) (- 1/2 (/ k 2)))
274.0ms
(- (+ (* 1/4 (* (log (* 2 PI)) (* (exp (* 1/2 (+ (log n) (log (* 2 PI))))) (* (log n) (pow k 2))))) (+ (* 1/8 (* (exp (* 1/2 (+ (log n) (log (* 2 PI))))) (* (pow (log n) 2) (pow k 2)))) (+ (exp (* 1/2 (+ (log n) (log (* 2 PI))))) (* 1/8 (* (pow (log (* 2 PI)) 2) (* (exp (* 1/2 (+ (log n) (log (* 2 PI))))) (pow k 2))))))) (+ (* 1/2 (* (exp (* 1/2 (+ (log n) (log (* 2 PI))))) (* (log n) k))) (* 1/2 (* (log (* 2 PI)) (* (exp (* 1/2 (+ (log n) (log (* 2 PI))))) k)))))

prune1.4s

Pruning

7 alts after pruning (6 fresh and 1 done)

Merged error: 0.0b

localize13.0ms

Local error

Found 4 expressions with local error:

0.5b
(pow (* PI (* n 2)) (- 1/2 (/ k 2)))
0.5b
(pow (* PI (* n 2)) (- 1/2 (/ k 2)))
0.3b
(* (sqrt (/ (pow (* PI (* n 2)) (- 1/2 (/ k 2))) (sqrt k))) (sqrt (/ (pow (* PI (* n 2)) (- 1/2 (/ k 2))) (sqrt k))))
0.2b
(* PI (* n 2))

rewrite26.0ms

Algorithm
rewrite-expression-head
Counts
4 → 114
Calls

4 calls. Slowest were:

14.0ms
(* (sqrt (/ (pow (* PI (* n 2)) (- 1/2 (/ k 2))) (sqrt k))) (sqrt (/ (pow (* PI (* n 2)) (- 1/2 (/ k 2))) (sqrt k))))
3.0ms
(pow (* PI (* n 2)) (- 1/2 (/ k 2)))
3.0ms
(pow (* PI (* n 2)) (- 1/2 (/ k 2)))

series577.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

216.0ms
(* (sqrt (/ (pow (* PI (* n 2)) (- 1/2 (/ k 2))) (sqrt k))) (sqrt (/ (pow (* PI (* n 2)) (- 1/2 (/ k 2))) (sqrt k))))
166.0ms
(pow (* PI (* n 2)) (- 1/2 (/ k 2)))
130.0ms
(pow (* PI (* n 2)) (- 1/2 (/ k 2)))
66.0ms
(* PI (* n 2))

simplify2.7s

Counts
95 → 126
Calls

95 calls. Slowest were:

242.0ms
(- (+ (* 1/4 (* (log (* 2 PI)) (* (exp (* 1/2 (+ (log n) (log (* 2 PI))))) (* (log n) (pow k 2))))) (+ (* 1/8 (* (exp (* 1/2 (+ (log n) (log (* 2 PI))))) (* (pow (log n) 2) (pow k 2)))) (+ (exp (* 1/2 (+ (log n) (log (* 2 PI))))) (* 1/8 (* (pow (log (* 2 PI)) 2) (* (exp (* 1/2 (+ (log n) (log (* 2 PI))))) (pow k 2))))))) (+ (* 1/2 (* (exp (* 1/2 (+ (log n) (log (* 2 PI))))) (* (log n) k))) (* 1/2 (* (log (* 2 PI)) (* (exp (* 1/2 (+ (log n) (log (* 2 PI))))) k)))))
241.0ms
(* (log (* PI (* n 2))) (- 1/2 (/ k 2)))
204.0ms
(* (log (* PI (* n 2))) (- 1/2 (/ k 2)))

prune1.9s

Pruning

7 alts after pruning (5 fresh and 2 done)

Merged error: 0.0b

localize26.0ms

Local error

Found 4 expressions with local error:

0.4b
(pow PI (- 1/2 (/ k 2)))
0.4b
(pow (* n 2) (- 1/2 (/ k 2)))
0.2b
(/ (pow (* n 2) (- 1/2 (/ k 2))) (sqrt k))
0.2b
(* (pow PI (- 1/2 (/ k 2))) (/ (pow (* n 2) (- 1/2 (/ k 2))) (sqrt k)))

rewrite30.0ms

Algorithm
rewrite-expression-head
Counts
4 → 121
Calls

4 calls. Slowest were:

13.0ms
(* (pow PI (- 1/2 (/ k 2))) (/ (pow (* n 2) (- 1/2 (/ k 2))) (sqrt k)))
9.0ms
(/ (pow (* n 2) (- 1/2 (/ k 2))) (sqrt k))
3.0ms
(pow (* n 2) (- 1/2 (/ k 2)))

series540.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

224.0ms
(* (pow PI (- 1/2 (/ k 2))) (/ (pow (* n 2) (- 1/2 (/ k 2))) (sqrt k)))
153.0ms
(/ (pow (* n 2) (- 1/2 (/ k 2))) (sqrt k))
118.0ms
(pow (* n 2) (- 1/2 (/ k 2)))
45.0ms
(pow PI (- 1/2 (/ k 2)))

simplify4.8s

Counts
116 → 133
Calls

116 calls. Slowest were:

455.0ms
(* (* (* (pow PI (- 1/2 (/ k 2))) (pow PI (- 1/2 (/ k 2)))) (pow PI (- 1/2 (/ k 2)))) (/ (* (* (pow (* n 2) (- 1/2 (/ k 2))) (pow (* n 2) (- 1/2 (/ k 2)))) (pow (* n 2) (- 1/2 (/ k 2)))) (* (* (sqrt k) (sqrt k)) (sqrt k))))
424.0ms
(+ (* (log PI) (- 1/2 (/ k 2))) (- (log (pow (* n 2) (- 1/2 (/ k 2)))) (log (sqrt k))))
382.0ms
(- (log (pow (* n 2) (- 1/2 (/ k 2)))) (log (sqrt k)))

prune1.7s

Pruning

6 alts after pruning (5 fresh and 1 done)

Merged error: 0b

regimes151.0ms

Accuracy

0% (0.5b remaining)

Error of 0.5b against oracle of 0.0b and baseline of 0.5b

bsearch2.0ms

end0.0ms

sample7.3s

Algorithm
intervals