def computeLikelihood(self, ctx): """ call the native code to compute log likelihood :param ctx: an instance of a ChainContext :returns: the log likelihood :raises InvalidLikelihoodException: in case the log likelihood from the wmapWrapper is smaller than 0.0001 """ cl_tt =ctx.get(CL_TT_KEY) cl_te =ctx.get(CL_TE_KEY) cl_ee =ctx.get(CL_EE_KEY) cl_bb =ctx.get(CL_BB_KEY) loglike = -wmapWrapperManager.computewmaplikelihood(cl_tt,cl_te,cl_ee,cl_bb) if(abs(loglike)<0.0001): raise InvalidLikelihoodException() return loglike
from wmap3Wrapper import wmapWrapperManager import numpy as np def copyArray(src, trg, index): lenght= min(len(src),len(trg)) for i in xrange(0,lenght): trg[i]=src[i,index] data = np.loadtxt("lensedcls.dat") size = 1199 cl_tt = np.zeros(size) cl_te = np.zeros(size) cl_ee = np.zeros(size) cl_bb = np.zeros(size) # copyArray(data, cl_tt, 1) copyArray(data, cl_te, 4) copyArray(data, cl_ee, 2) copyArray(data, cl_bb, 3) wmapWrapperManager.setup() likelihood = wmapWrapperManager.computewmaplikelihood(cl_tt,cl_te,cl_ee,cl_bb) print likelihood