Skip to content
GitLab
Projects
Groups
Snippets
/
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
Gregory Ashton
PyFstat
Commits
90df76df
Commit
90df76df
authored
Feb 16, 2017
by
Gregory Ashton
Browse files
Update fc on glitching data example
parent
8b5d333d
Changes
4
Hide whitespace changes
Inline
Side-by-side
docs/fully_coherent_search_using_MCMC_on_glitching_data.md
View file @
90df76df
...
...
@@ -39,7 +39,7 @@ between the different peaks in the posterior).
ntemps = 2
log10temperature_min = -0.01
nwalkers = 100
nsteps = [500
0
,
100
00]
nsteps = [500,
5
00]
mcmc = MCMCSearch('fully_coherent_search_using_MCMC_on_glitching_data', 'data',
sftfilepath='data/*_glitch*.sft',
...
...
docs/img/fully_coherent_search_using_MCMC_on_glitching_data_corner.png
View replaced file @
8b5d333d
View file @
90df76df
119 KB
|
W:
|
H:
119 KB
|
W:
|
H:
2-up
Swipe
Onion skin
docs/img/fully_coherent_search_using_MCMC_on_glitching_data_walkers.png
View replaced file @
8b5d333d
View file @
90df76df
74 KB
|
W:
|
H:
193 KB
|
W:
|
H:
2-up
Swipe
Onion skin
examples/fully_coherent_search_using_MCMC_on_glitching_data.py
View file @
90df76df
import
numpy
as
np
from
pyfstat
import
MCMCSearch
F0
=
30.0
F1
=
-
1e-10
F2
=
0
Alpha
=
5e-3
Delta
=
6e-2
Alpha
=
np
.
radians
(
83.6292
)
Delta
=
np
.
radians
(
22.0144
)
tref
=
362750407.0
tstart
=
1000000000
...
...
@@ -21,11 +23,11 @@ theta_prior = {'F0': {'type': 'unif', 'lower': F0-1e-4, 'upper': F0+1e-4},
ntemps
=
2
log10temperature_min
=
-
0.01
nwalkers
=
100
nsteps
=
[
500
0
,
100
00
]
nsteps
=
[
500
,
5
00
]
mcmc
=
MCMCSearch
(
'fully_coherent_search_using_MCMC_on_glitching_data'
,
'data'
,
sftfilepath
=
'data/*_glitch*.sft'
,
theta_prior
=
theta_prior
,
tref
=
tref
,
tstart
=
tstart
,
tend
=
tend
,
theta_prior
=
theta_prior
,
tref
=
tref
,
minStartTime
=
tstart
,
maxStartTime
=
tend
,
nsteps
=
nsteps
,
nwalkers
=
nwalkers
,
ntemps
=
ntemps
,
log10temperature_min
=
log10temperature_min
)
mcmc
.
run
()
...
...
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment