Average Error: 0.3 → 0.4
Time: 1.1m
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{e^{\left(\frac{1}{2} - \frac{k}{2}\right) \cdot \log \pi}}{\frac{\sqrt{k}}{{\left(2 \cdot n\right)}^{\left(\frac{1}{2} - \frac{k}{2}\right)}}}\]

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.3

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

  2. Simplified0.3

    \[\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 unpow-prod-down0.5

    \[\leadsto \frac{\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)}}}{\sqrt{k}}\]

  5. Applied associate-/l*0.4

    \[\leadsto \color{blue}{\frac{{\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)}}}}\]
  6. Using strategy rm

  7. Applied pow-to-exp0.4

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

  8. Final simplification0.4

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

Reproduce

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

Details

Time bar (total: 1.1m)Debug log

sample207.0ms

Algorithm
intervals
Results
115.0ms603×body80nan
65.0ms392×body80valid

simplify36.0ms

Counts
1 → 1
Calls
1 calls:
Slowest
36.0ms
(* (/ 1 (sqrt k)) (pow (* (* 2 PI) n) (/ (- 1 k) 2)))

prune19.0ms

Pruning

2 alts after pruning (2 fresh and 0 done)

Merged error: 0.2b

localize31.0ms

Local error

Found 3 expressions with local error:

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

rewrite10.0ms

Algorithm
rewrite-expression-head
Rules
18×times-frac
15×add-cube-cbrt
15×*-un-lft-identity
15×add-sqr-sqrt
10×unpow-prod-up
9×add-exp-log
6×associate-/l*
5×add-cbrt-cube
5×unpow-prod-down
5×sub-neg
4×pow1
3×add-log-exp
3×div-exp
3×associate-/r*
3×pow-unpow
2×pow-sub
2×pow-exp
2×pow-to-exp
1×associate-/l/
1×div-inv
1×*-commutative
1×associate-*r*
1×frac-2neg
1×pow-pow
1×clear-num
1×cbrt-undiv
Counts
3 → 68
Calls
3 calls:
Slowest
6.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))

series409.0ms

Counts
3 → 9
Calls
3 calls:
Slowest
209.0ms
(/ (pow (* PI (* n 2)) (- 1/2 (/ k 2))) (sqrt k))
132.0ms
(pow (* PI (* n 2)) (- 1/2 (/ k 2)))
68.0ms
(* PI (* n 2))

simplify3.1s

Counts
67 → 77
Calls
67 calls:
Slowest
594.0ms
(- (* (log (* PI (* n 2))) (- 1/2 (/ k 2))) (log (sqrt k)))
559.0ms
(- (* (log (* PI (* n 2))) (- 1/2 (/ k 2))) (log (sqrt k)))
265.0ms
(* (log (* PI (* n 2))) (- 1/2 (/ k 2)))
263.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)))
225.0ms
(- (log (pow (* PI (* n 2)) (- 1/2 (/ k 2)))) (log (sqrt k)))

prune837.0ms

Pruning

8 alts after pruning (8 fresh and 0 done)

Merged error: 0.0b

localize29.0ms

Local error

Found 4 expressions with local error:

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

rewrite31.0ms

Algorithm
rewrite-expression-head
Rules
241×times-frac
101×add-cube-cbrt
101×*-un-lft-identity
101×add-sqr-sqrt
96×unpow-prod-up
48×sub-neg
29×associate-/r*
26×add-exp-log
23×unpow-prod-down
17×div-exp
11×add-cbrt-cube
10×pow-sub
9×pow-to-exp
8×associate-/l*
8×div-inv
8×associate-/r/
6×pow-unpow
5×pow1
4×add-log-exp
4×pow-exp
4×cbrt-undiv
2×frac-2neg
2×clear-num
1×associate-/l/
1×pow-pow
Counts
4 → 238
Calls
4 calls:
Slowest
12.0ms
(/ (pow PI (- 1/2 (/ k 2))) (/ (sqrt k) (pow (* n 2) (- 1/2 (/ k 2)))))
7.0ms
(/ (sqrt k) (pow (* n 2) (- 1/2 (/ k 2))))
2.0ms
(pow (* n 2) (- 1/2 (/ k 2)))
1.0ms
(pow PI (- 1/2 (/ k 2)))

series538.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
267.0ms
(/ (pow PI (- 1/2 (/ k 2))) (/ (sqrt k) (pow (* n 2) (- 1/2 (/ k 2)))))
119.0ms
(pow (* n 2) (- 1/2 (/ k 2)))
112.0ms
(/ (sqrt k) (pow (* n 2) (- 1/2 (/ k 2))))
38.0ms
(pow PI (- 1/2 (/ k 2)))

simplify12.6s

Counts
345 → 250
Calls
345 calls:
Slowest
700.0ms
(/ (* (cbrt (pow PI (- 1/2 (/ k 2)))) (cbrt (pow PI (- 1/2 (/ k 2))))) (/ (* (cbrt (sqrt k)) (cbrt (sqrt k))) (* (cbrt (pow (* n 2) (- 1/2 (/ k 2)))) (cbrt (pow (* n 2) (- 1/2 (/ k 2)))))))
487.0ms
(- (log (pow PI (- 1/2 (/ k 2)))) (log (/ (sqrt k) (pow (* n 2) (- 1/2 (/ k 2))))))
477.0ms
(/ (* (* (pow PI (- 1/2 (/ k 2))) (pow PI (- 1/2 (/ k 2)))) (pow PI (- 1/2 (/ k 2)))) (/ (* (* (sqrt k) (sqrt k)) (sqrt k)) (* (* (pow (* n 2) (- 1/2 (/ k 2))) (pow (* n 2) (- 1/2 (/ k 2)))) (pow (* n 2) (- 1/2 (/ k 2))))))
432.0ms
(- (* (log PI) (- 1/2 (/ k 2))) (- (log (sqrt k)) (log (pow (* n 2) (- 1/2 (/ k 2))))))
397.0ms
(- (log (pow PI (- 1/2 (/ k 2)))) (- (log (sqrt k)) (log (pow (* n 2) (- 1/2 (/ k 2))))))

prune3.1s

Pruning

7 alts after pruning (6 fresh and 1 done)

Merged error: 0b

localize22.0ms

Local error

Found 4 expressions with local error:

1.1b
(pow (* n 2) (- 1/2 (/ k 2)))
0.2b
(/ (exp (* (log PI) (- 1/2 (/ k 2)))) (/ (sqrt k) (pow (* n 2) (- 1/2 (/ k 2)))))
0.2b
(/ (sqrt k) (pow (* n 2) (- 1/2 (/ k 2))))
0.0b
(* (log PI) (- 1/2 (/ k 2)))

rewrite34.0ms

Algorithm
rewrite-expression-head
Rules
323×times-frac
122×add-cube-cbrt
122×*-un-lft-identity
122×add-sqr-sqrt
96×exp-sum
79×sub-neg
58×unpow-prod-up
50×distribute-lft-in
50×distribute-rgt-in
29×unpow-prod-down
29×associate-/r*
18×add-exp-log
13×add-cbrt-cube
10×pow-sub
10×associate-/l*
10×div-inv
10×associate-/r/
10×div-exp
7×pow1
4×add-log-exp
4×cbrt-undiv
3×pow-exp
3×associate-*l*
3×associate-*r*
3×pow-unpow
3×pow-to-exp
2×associate-*r/
2×frac-2neg
2×clear-num
1×flip--
1×cbrt-unprod
1×*-commutative
1×prod-exp
1×pow-prod-down
1×flip3--
1×pow-pow
Counts
4 → 290
Calls
4 calls:
Slowest
17.0ms
(/ (exp (* (log PI) (- 1/2 (/ k 2)))) (/ (sqrt k) (pow (* n 2) (- 1/2 (/ k 2)))))
4.0ms
(* (log PI) (- 1/2 (/ k 2)))
4.0ms
(/ (sqrt k) (pow (* n 2) (- 1/2 (/ k 2))))
3.0ms
(pow (* n 2) (- 1/2 (/ k 2)))

series804.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
323.0ms
(* (log PI) (- 1/2 (/ k 2)))
268.0ms
(/ (exp (* (log PI) (- 1/2 (/ k 2)))) (/ (sqrt k) (pow (* n 2) (- 1/2 (/ k 2)))))
118.0ms
(/ (sqrt k) (pow (* n 2) (- 1/2 (/ k 2))))
96.0ms
(pow (* n 2) (- 1/2 (/ k 2)))

simplify27.2s

Counts
442 → 302
Calls
442 calls:
Slowest
435.0ms
(/ (* (* (exp (* (log PI) (- 1/2 (/ k 2)))) (exp (* (log PI) (- 1/2 (/ k 2))))) (exp (* (log PI) (- 1/2 (/ k 2))))) (/ (* (* (sqrt k) (sqrt k)) (sqrt k)) (* (* (pow (* n 2) (- 1/2 (/ k 2))) (pow (* n 2) (- 1/2 (/ k 2)))) (pow (* n 2) (- 1/2 (/ k 2))))))
379.0ms
(/ (* (cbrt (exp (* (log PI) (- 1/2 (/ k 2))))) (cbrt (exp (* (log PI) (- 1/2 (/ k 2)))))) (/ (* (cbrt (sqrt k)) (cbrt (sqrt k))) (* (cbrt (pow (* n 2) (- 1/2 (/ k 2)))) (cbrt (pow (* n 2) (- 1/2 (/ k 2)))))))
357.0ms
(* (* (* (log PI) (log PI)) (log PI)) (* (* (- 1/2 (/ k 2)) (- 1/2 (/ k 2))) (- 1/2 (/ k 2))))
323.0ms
(- (* (log PI) (- 1/2 (/ k 2))) (- (log (sqrt k)) (log (pow (* n 2) (- 1/2 (/ k 2))))))
303.0ms
(/ (* (cbrt (exp (* (log PI) (- 1/2 (/ k 2))))) (cbrt (exp (* (log PI) (- 1/2 (/ k 2)))))) (/ 1 (* (cbrt (pow (* n 2) (- 1/2 (/ k 2)))) (cbrt (pow (* n 2) (- 1/2 (/ k 2)))))))

prune4.1s

Pruning

7 alts after pruning (5 fresh and 2 done)

Merged error: 0b

localize38.0ms

Local error

Found 4 expressions with local error:

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

rewrite21.0ms

Algorithm
rewrite-expression-head
Rules
40×times-frac
28×add-cube-cbrt
28×*-un-lft-identity
28×add-sqr-sqrt
24×unpow-prod-up
14×add-exp-log
12×sub-neg
10×pow-sub
10×associate-/r*
8×associate-/l*
8×div-inv
8×associate-/r/
8×div-exp
8×add-cbrt-cube
8×pow-to-exp
6×pow-unpow
4×add-log-exp
4×rec-exp
4×pow1
2×frac-2neg
2×clear-num
2×cbrt-undiv
1×associate-/l/
Counts
4 → 119
Calls
4 calls:
Slowest
8.0ms
(/ (pow PI (- 1/2 (/ k 2))) (/ 1 (pow n (- 1/2 (/ k 2)))))
7.0ms
(/ (sqrt k) (pow 2 (- 1/2 (/ k 2))))
1.0ms
(pow n (- 1/2 (/ k 2)))
1.0ms
(pow PI (- 1/2 (/ k 2)))

series410.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
145.0ms
(/ (pow PI (- 1/2 (/ k 2))) (/ 1 (pow n (- 1/2 (/ k 2)))))
127.0ms
(/ (sqrt k) (pow 2 (- 1/2 (/ k 2))))
102.0ms
(pow n (- 1/2 (/ k 2)))
35.0ms
(pow PI (- 1/2 (/ k 2)))

simplify5.3s

Counts
133 → 131
Calls
133 calls:
Slowest
603.0ms
(- (log (pow PI (- 1/2 (/ k 2)))) (- (log (pow n (- 1/2 (/ k 2))))))
471.0ms
(* (exp (* -1 (* (- 1/2 (* 1/2 k)) (log (/ 1 n))))) (pow PI (- 1/2 (* 1/2 k))))
466.0ms
(* (exp (* (- 1/2 (* 1/2 k)) (- (log -1) (log (/ -1 n))))) (pow PI (- 1/2 (* 1/2 k))))
451.0ms
(exp (* (- 1/2 (* 1/2 k)) (- (log -1) (log (/ -1 n)))))
278.0ms
(/ (* (* (pow PI (- 1/2 (/ k 2))) (pow PI (- 1/2 (/ k 2)))) (pow PI (- 1/2 (/ k 2)))) (* (* (/ 1 (pow n (- 1/2 (/ k 2)))) (/ 1 (pow n (- 1/2 (/ k 2))))) (/ 1 (pow n (- 1/2 (/ k 2))))))

prune1.6s

Pruning

7 alts after pruning (4 fresh and 3 done)

Merged error: 0b

regimes123.0ms

Accuracy

0% (0.4b remaining)

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

bsearch1.0ms

end0.0ms

sample7.1s

Algorithm
intervals
Results
3.6s20818×body80nan
2.7s11839×body80valid
0.0ms1×body160nan