def createChainContext(self, p):
"""
Returns a new instance of a chain context
"""
try:
p = Params(*zip(self.params.keys, p))
except Exception:
# no params or params has no keys
pass
return ChainContext(self, p)
wmap9Likelihood = wmap9.WmapLikelihoodModule()
#setting up wrappers
cambCore.setup()
pyCambCore.setup()
wmap3Likelihood.setup()
wmap5Likelihood.setup()
wmap7Likelihood.setup()
wmap9Likelihood.setup()
#camb wrapper
print("setup done. calling CAMB power spectrum computation")
paramValues = [70, 0.0226, 0.122, 2.1E-009, 0.96, 0.09]
ctx = ChainContext(None, paramValues)
cambCore(ctx)
print("cls computation done. calling WMAP likelihood")
print("Likelihood WMAP 3: %s" % wmap3Likelihood.computeLikelihood(ctx))
print("Likelihood WMAP 5: %s" % wmap5Likelihood.computeLikelihood(ctx))
print("Likelihood WMAP 7: %s" % wmap7Likelihood.computeLikelihood(ctx))
print("Likelihood WMAP 9: %s" % wmap9Likelihood.computeLikelihood(ctx))
#pycamb
ctx = ChainContext(None, paramValues)
pyCambCore(ctx)
print("cls computation done. calling WMAP likelihood")
print("Likelihood WMAP 3: %s" % wmap3Likelihood.computeLikelihood(ctx))
print("Likelihood WMAP 5: %s" % wmap5Likelihood.computeLikelihood(ctx))
def setup(self):
self.ctx=ChainContext(self, self.params)
def createChainContext(self, p):
"""
Returns a new instance of a chain context
"""
return ChainContext(self, p)
camb = PyCambCoreModule(lmax=2700)
camb.setup()
# Path to clik files
path = ''
clik_files = [
join(path, 'commander_v4.1_lm49.clik'),
join(path, 'CAMspec_v6.2TN_2013_02_26_dist.clik'),
join(path, 'lowlike_v222.clik')
]
# clik_files = join(path,'commander_v4.1_lm49.clik')
params = array([
67.04, 0.022032, 0.12038,
exp(3.0980) * 1e-10, 0.9619, 0.0925, 152, 63.3, 117.0, 0.0, 27.2, 6.80,
0.916, 0.406, 0.601, 1.0, 1.0, 0.03, 0.9, 0.0
])
nuisance_indices = [[], range(6, 6 + 14), []]
# nuisance_indices = [[]]
like = PlanckLikelihoodModule(clik_files, nuisance_indices)
like.setup()
ctx = ChainContext(None, params)
camb(ctx)
print like.computeLikelihood(ctx)
