Average Error: 0.4 → 0.4
Time: 1.0m
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)}\]
\[{\left(2 \cdot n\right)}^{\left(\frac{1}{2} - \frac{k}{2}\right)} \cdot \frac{\frac{1}{\frac{1}{{\pi}^{\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 clear-num0.4

    \[\leadsto \color{blue}{\frac{1}{\frac{\sqrt{k}}{{\left(\pi \cdot \left(n \cdot 2\right)\right)}^{\left(\frac{1}{2} - \frac{k}{2}\right)}}}}\]
  5. Using strategy rm
  6. Applied unpow-prod-down0.5

    \[\leadsto \frac{1}{\frac{\sqrt{k}}{\color{blue}{{\pi}^{\left(\frac{1}{2} - \frac{k}{2}\right)} \cdot {\left(n \cdot 2\right)}^{\left(\frac{1}{2} - \frac{k}{2}\right)}}}}\]
  7. Applied *-un-lft-identity0.5

    \[\leadsto \frac{1}{\frac{\color{blue}{1 \cdot \sqrt{k}}}{{\pi}^{\left(\frac{1}{2} - \frac{k}{2}\right)} \cdot {\left(n \cdot 2\right)}^{\left(\frac{1}{2} - \frac{k}{2}\right)}}}\]
  8. Applied times-frac0.4

    \[\leadsto \frac{1}{\color{blue}{\frac{1}{{\pi}^{\left(\frac{1}{2} - \frac{k}{2}\right)}} \cdot \frac{\sqrt{k}}{{\left(n \cdot 2\right)}^{\left(\frac{1}{2} - \frac{k}{2}\right)}}}}\]
  9. Applied associate-/r*0.4

    \[\leadsto \color{blue}{\frac{\frac{1}{\frac{1}{{\pi}^{\left(\frac{1}{2} - \frac{k}{2}\right)}}}}{\frac{\sqrt{k}}{{\left(n \cdot 2\right)}^{\left(\frac{1}{2} - \frac{k}{2}\right)}}}}\]
  10. Using strategy rm
  11. Applied associate-/r/0.4

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

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

Reproduce

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

Details

Time bar (total: 1.0m)Debug log

sample296.0ms

Algorithm
intervals

simplify48.0ms

Counts
1 → 1
Calls

1 calls. Slowest were:

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

prune13.0ms

Pruning

2 alts after pruning (2 fresh and 0 done)

Merged error: 0.3b

localize41.0ms

Local error

Found 3 expressions with local error:

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

rewrite16.0ms

Algorithm
rewrite-expression-head
Counts
3 → 68
Calls

3 calls. Slowest were:

9.0ms
(/ (pow (* PI (* n 2)) (- 1/2 (/ k 2))) (sqrt k))
4.0ms
(pow (* PI (* n 2)) (- 1/2 (/ k 2)))
2.0ms
(* PI (* n 2))

series428.0ms

Counts
3 → 9
Calls

3 calls. Slowest were:

223.0ms
(/ (pow (* PI (* n 2)) (- 1/2 (/ k 2))) (sqrt k))
122.0ms
(pow (* PI (* n 2)) (- 1/2 (/ k 2)))
83.0ms
(* PI (* n 2))

simplify3.6s

Counts
67 → 77
Calls

67 calls. Slowest were:

874.0ms
(- (* (log (* PI (* n 2))) (- 1/2 (/ k 2))) (log (sqrt k)))
588.0ms
(- (* (log (* PI (* n 2))) (- 1/2 (/ k 2))) (log (sqrt k)))
287.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)))

prune998.0ms

Pruning

10 alts after pruning (10 fresh and 0 done)

Merged error: 0.0b

localize12.0ms

Local error

Found 4 expressions with local error:

1.4b
(pow (* PI (* n 2)) (- 1/2 (/ k 2)))
0.3b
(* PI (* n 2))
0.1b
(/ (sqrt k) (pow (* PI (* n 2)) (- 1/2 (/ k 2))))
0.1b
(/ 1 (/ (sqrt k) (pow (* PI (* n 2)) (- 1/2 (/ k 2)))))

rewrite22.0ms

Algorithm
rewrite-expression-head
Counts
4 → 108
Calls

4 calls. Slowest were:

9.0ms
(/ (sqrt k) (pow (* PI (* n 2)) (- 1/2 (/ k 2))))
5.0ms
(/ 1 (/ (sqrt k) (pow (* PI (* n 2)) (- 1/2 (/ k 2)))))
4.0ms
(pow (* PI (* n 2)) (- 1/2 (/ k 2)))

series549.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

171.0ms
(/ 1 (/ (sqrt k) (pow (* PI (* n 2)) (- 1/2 (/ k 2)))))
155.0ms
(/ (sqrt k) (pow (* PI (* n 2)) (- 1/2 (/ k 2))))
147.0ms
(pow (* PI (* n 2)) (- 1/2 (/ k 2)))
76.0ms
(* PI (* n 2))

simplify3.5s

Counts
95 → 120
Calls

95 calls. Slowest were:

405.0ms
(- (log (sqrt k)) (* (log (* PI (* n 2))) (- 1/2 (/ k 2))))
381.0ms
(- (log (sqrt k)) (* (log (* PI (* n 2))) (- 1/2 (/ k 2))))
346.0ms
(* (log (* PI (* n 2))) (- 1/2 (/ k 2)))

prune1.6s

Pruning

8 alts after pruning (8 fresh and 0 done)

Merged error: 0.0b

localize9.0ms

Local error

Found 4 expressions with local error:

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

rewrite11.0ms

Algorithm
rewrite-expression-head
Counts
4 → 93
Calls

4 calls. Slowest were:

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

series429.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

172.0ms
(/ (/ 1 (sqrt k)) (/ 1 (pow (* PI (* n 2)) (- 1/2 (/ k 2)))))
167.0ms
(pow (* PI (* n 2)) (- 1/2 (/ k 2)))
71.0ms
(* PI (* n 2))
19.0ms
(/ 1 (sqrt k))

simplify5.0s

Counts
83 → 105
Calls

83 calls. Slowest were:

1.4s
(/ (* (* (/ 1 (sqrt k)) (/ 1 (sqrt k))) (/ 1 (sqrt k))) (* (* (/ 1 (pow (* PI (* n 2)) (- 1/2 (/ k 2)))) (/ 1 (pow (* PI (* n 2)) (- 1/2 (/ k 2))))) (/ 1 (pow (* PI (* n 2)) (- 1/2 (/ k 2))))))
338.0ms
(- (- (log (sqrt k))) (- (log (pow (* PI (* n 2)) (- 1/2 (/ k 2))))))
276.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 (7 fresh and 0 done)

Merged error: 0.0b

localize44.0ms

Local error

Found 4 expressions with local error:

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

rewrite35.0ms

Algorithm
rewrite-expression-head
Counts
4 → 246
Calls

4 calls. Slowest were:

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

series562.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

298.0ms
(/ (/ 1 (/ 1 (pow PI (- 1/2 (/ k 2))))) (/ (sqrt k) (pow (* n 2) (- 1/2 (/ k 2)))))
132.0ms
(/ (sqrt k) (pow (* n 2) (- 1/2 (/ k 2))))
107.0ms
(pow (* n 2) (- 1/2 (/ k 2)))
24.0ms
(pow PI (- 1/2 (/ k 2)))

simplify30.8s

Counts
351 → 258
Calls

351 calls. Slowest were:

773.0ms
(/ (/ 1 (/ 1 (pow PI (- 1/2 (/ k 2))))) (/ (sqrt k) (pow (* n 2) (- 1/2 (/ k 2)))))
524.0ms
(- (- (log (/ 1 (pow PI (- 1/2 (/ k 2)))))) (- (log (sqrt k)) (log (pow (* n 2) (- 1/2 (/ k 2))))))
486.0ms
(/ (/ 1 (/ 1 (pow PI (- 1/2 (/ k 2))))) (/ (sqrt (sqrt k)) (pow (* n 2) (- (/ k 2)))))

prune3.7s

Pruning

6 alts after pruning (6 fresh and 0 done)

Merged error: 0b

regimes138.0ms

Accuracy

0% (0.4b remaining)

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

bsearch4.0ms

end0.0ms

sample8.0s

Algorithm
intervals