def setup(self): """ Sets up the cmb likelihood wrapper """ wmapWrapperManager.setup(self.path)
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