diff --git a/docs/fully_coherent_search_using_MCMC_on_glitching_data.md b/docs/fully_coherent_search_using_MCMC_on_glitching_data.md index 80073253fe6f50dd3f948ae327c43e30f10d137b..4360d4f5fb3c3bf77d5e47dc02e06cd2867ea232 100644 --- a/docs/fully_coherent_search_using_MCMC_on_glitching_data.md +++ b/docs/fully_coherent_search_using_MCMC_on_glitching_data.md @@ -39,7 +39,7 @@ between the different peaks in the posterior). ntemps = 2 log10temperature_min = -0.01 nwalkers = 100 -nsteps = [5000, 10000] +nsteps = [500, 500] mcmc = MCMCSearch('fully_coherent_search_using_MCMC_on_glitching_data', 'data', sftfilepath='data/*_glitch*.sft', diff --git a/docs/img/fully_coherent_search_using_MCMC_on_glitching_data_corner.png b/docs/img/fully_coherent_search_using_MCMC_on_glitching_data_corner.png index a3812bfd1aa21d5add813c9405d37cc1b3ea5639..d5956109ad50e01f98cc941cbd0a280ff9faa898 100644 Binary files a/docs/img/fully_coherent_search_using_MCMC_on_glitching_data_corner.png and b/docs/img/fully_coherent_search_using_MCMC_on_glitching_data_corner.png differ diff --git a/docs/img/fully_coherent_search_using_MCMC_on_glitching_data_walkers.png b/docs/img/fully_coherent_search_using_MCMC_on_glitching_data_walkers.png index 6ad24a11b4e4488a071cae1f861e44bef4f15fd1..c5f3da6c366a796ae9ac54801448d0cd3bb169d5 100644 Binary files a/docs/img/fully_coherent_search_using_MCMC_on_glitching_data_walkers.png and b/docs/img/fully_coherent_search_using_MCMC_on_glitching_data_walkers.png differ diff --git a/examples/fully_coherent_search_using_MCMC_on_glitching_data.py b/examples/fully_coherent_search_using_MCMC_on_glitching_data.py index e7fbeea1251924f6416f8816a8dea8fff67d1103..75c35cbd33573b335737771964d9b172fb7aefbc 100644 --- a/examples/fully_coherent_search_using_MCMC_on_glitching_data.py +++ b/examples/fully_coherent_search_using_MCMC_on_glitching_data.py @@ -1,10 +1,12 @@ +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 = [5000, 10000] +nsteps = [500, 500] 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()