Skip to content

Commit

Permalink
Merge branch 'master' of github.com:jeremy-baier/enterprise_extension…
Browse files Browse the repository at this point in the history
…s into jgb-15yr-cnm
  • Loading branch information
jeremy-baier committed Sep 8, 2024
2 parents 4d48771 + 6887f81 commit 0b30f05
Show file tree
Hide file tree
Showing 3 changed files with 13 additions and 4 deletions.
2 changes: 2 additions & 0 deletions .flake8
Original file line number Diff line number Diff line change
Expand Up @@ -2,6 +2,8 @@
max-line-length = 120
max-complexity = 45
ignore =
C901,
# function is too complex
E203,
# missing whitespace around operator
E225,
Expand Down
13 changes: 10 additions & 3 deletions enterprise_extensions/sampler.py
Original file line number Diff line number Diff line change
Expand Up @@ -166,6 +166,15 @@ def __init__(self, pta, snames=None, empirical_distr=None, f_stat_file=None, sav
self.ndim = sum(p.size or 1 for p in pta.params)
self.plist = [p.name for p in pta.params]

# parameter dictionary
self.params_dict = {}
for p in self.params:
if p.size:
for ii in range(0, p.size):
self.params_dict.update({p.name + "_{}".format(ii): p})
else:
self.params_dict.update({p.name: p})

# parameter map
self.pmap = {}
ct = 0
Expand Down Expand Up @@ -616,9 +625,7 @@ def draw_from_gwb_log_uniform_distribution(self, x, iter, beta):
signal_name = [par for par in self.pnames
if ('gw' in par and 'log10_A' in par)][0]

param_names = [par.name for par in self.params]
idx = list(param_names).index(signal_name)
param = self.params[idx]
param = self.params_dict[signal_name]

q[self.pmap[str(param)]] = np.random.uniform(param.prior._defaults['pmin'], param.prior._defaults['pmax'])

Expand Down
2 changes: 1 addition & 1 deletion requirements.txt
Original file line number Diff line number Diff line change
Expand Up @@ -6,7 +6,7 @@ scikit-sparse>=0.4.5
pint-pulsar>=0.8.2
libstempo>=2.4.0
enterprise-pulsar>=3.3.0
scikit-learn==0.24
scikit-learn==1.0.1
emcee
ptmcmcsampler
numdifftools

0 comments on commit 0b30f05

Please sign in to comment.