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)