Commit e725a15f authored by Gregory Ashton's avatar Gregory Ashton
Browse files

Update examples to match paper

parent 8d979b1b
...@@ -7,7 +7,7 @@ plt.style.use('paper') ...@@ -7,7 +7,7 @@ plt.style.use('paper')
label = 'semicoherent_glitch_robust_directed_MCMC_search_on_1_glitch' label = 'semicoherent_glitch_robust_directed_MCMC_search_on_1_glitch'
Nstar = 100 Nstar = 1000
F0_width = np.sqrt(Nstar)*np.sqrt(12)/(np.pi*duration) F0_width = np.sqrt(Nstar)*np.sqrt(12)/(np.pi*duration)
F1_width = np.sqrt(Nstar)*np.sqrt(180)/(np.pi*duration**2) F1_width = np.sqrt(Nstar)*np.sqrt(180)/(np.pi*duration**2)
...@@ -30,7 +30,7 @@ theta_prior = { ...@@ -30,7 +30,7 @@ theta_prior = {
'Delta': Delta, 'Delta': Delta,
} }
ntemps = 2 ntemps = 3
log10beta_min = -0.5 log10beta_min = -0.5
nwalkers = 100 nwalkers = 100
nsteps = [500, 500] nsteps = [500, 500]
......
...@@ -14,10 +14,10 @@ label = 'semicoherent_glitch_robust_directed_grid_search_on_1_glitch' ...@@ -14,10 +14,10 @@ label = 'semicoherent_glitch_robust_directed_grid_search_on_1_glitch'
plt.style.use('paper') plt.style.use('paper')
Nstar = 150 Nstar = 1000
F0_width = np.sqrt(Nstar)*np.sqrt(12)/(np.pi*duration) F0_width = np.sqrt(Nstar)*np.sqrt(12)/(np.pi*duration)
F1_width = np.sqrt(Nstar)*np.sqrt(180)/(np.pi*duration**2) F1_width = np.sqrt(Nstar)*np.sqrt(180)/(np.pi*duration**2)
N = 61 N = 30
F0s = [F0-F0_width/2., F0+F0_width/2., F0_width/N] F0s = [F0-F0_width/2., F0+F0_width/2., F0_width/N]
F1s = [F1-F1_width/2., F1+F1_width/2., F1_width/N] F1s = [F1-F1_width/2., F1+F1_width/2., F1_width/N]
F2s = [F2] F2s = [F2]
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment