int_HMF2.M.copy(), int_HMF2.zarr[zbins], param_vals2)[-1, -1])
#pl = Plotter()
#pl.add(np.log10(int_HMF.M200[:,zbins]),np.log10(int_HMF.dn_dM(int_HMF.M200,200)[:,zbins]*int_HMF.M200[:,zbins]),color='b',alpha=0.9)
#pl.add(np.log10(int_HMF2.M200[:,zbins]),np.log10(int_HMF2.dn_dM(int_HMF2.M200,200)[:,zbins]*int_HMF2.M200[:,zbins]),color='r',alpha=0.9)
#pl.add(np.log10(int_HMF.M200[:,zbins]),np.log10(dn_dzdm_int(int_HMF.zarr[zbins],np.log10(int_HMF.M.copy()))[:,0]*int_HMF.M200[:,zbins]),linestyle='--',color='g',alpha=0.9)
#pl.add(np.log10(int_HMF2.M200[:,zbins]),np.log10(dn_dzdm_int2(int_HMF2.zarr[zbins],np.log10(int_HMF2.M.copy()))[:,0]*int_HMF2.M200[:,zbins]),linestyle='--',color='c',alpha=0.9)
#pl.done("test_MF.png")
#print np.log10(dn_dzdm_int(int_HMF.zarr[zbins],np.log10(int_HMF.M.copy()))[:,0]*int_HMF.M200[:,zbins])
print(
"MF interp check",
np.sum(
np.log10(
int_HMF.dn_dM(int_HMF.M200, 200)[:, zbins] *
int_HMF.M200[:, zbins]) - np.log10(
dn_dzdm_int(int_HMF.zarr[zbins], np.log10(
int_HMF.M.copy()))[:, 0] * int_HMF.M200[:, zbins])))
print(int_HMF.zarr[zbins])
start = time.time()
print(CL.lnlike(parvals, parlist)) #,priorvals,priorlist)
print("like call", time.time() - start)
print("prior", CL.lnprior(parvals, parlist, priorvals, priorlist))
#print "Prob",CL.lnprob(parvals,parlist,priorvals,priorlist)
start = time.time()
print(CL.lnlike(parvals2, parlist)) #,priorvals,priorlist)
print("like call", time.time() - start)
print(old_div(Y[-1,-1],Ytilde[-1,-1]))
print(np.log(old_div(Y[-1,-1],Ytilde[-1,-1])))
#print -1.*(np.log(Y[-1,30:35]/Ytilde[-1,30:35]))
print("P of Y", CL.P_Yo(np.outer(np.ones(len(int_HMF.M.copy())),CL.LgY),int_HMF2.M.copy(),int_HMF2.zarr[zbins],param_vals2)[-1,-1])
#pl = Plotter()
#pl.add(np.log10(int_HMF.M200[:,zbins]),np.log10(int_HMF.dn_dM(int_HMF.M200,200)[:,zbins]*int_HMF.M200[:,zbins]),color='b',alpha=0.9)
#pl.add(np.log10(int_HMF2.M200[:,zbins]),np.log10(int_HMF2.dn_dM(int_HMF2.M200,200)[:,zbins]*int_HMF2.M200[:,zbins]),color='r',alpha=0.9)
#pl.add(np.log10(int_HMF.M200[:,zbins]),np.log10(dn_dzdm_int(int_HMF.zarr[zbins],np.log10(int_HMF.M.copy()))[:,0]*int_HMF.M200[:,zbins]),linestyle='--',color='g',alpha=0.9)
#pl.add(np.log10(int_HMF2.M200[:,zbins]),np.log10(dn_dzdm_int2(int_HMF2.zarr[zbins],np.log10(int_HMF2.M.copy()))[:,0]*int_HMF2.M200[:,zbins]),linestyle='--',color='c',alpha=0.9)
#pl.done("test_MF.png")
#print np.log10(dn_dzdm_int(int_HMF.zarr[zbins],np.log10(int_HMF.M.copy()))[:,0]*int_HMF.M200[:,zbins])
print("MF interp check", np.sum(np.log10(int_HMF.dn_dM(int_HMF.M200,200)[:,zbins]*int_HMF.M200[:,zbins]) - np.log10(dn_dzdm_int(int_HMF.zarr[zbins],np.log10(int_HMF.M.copy()))[:,0]*int_HMF.M200[:,zbins])))
print(int_HMF.zarr[zbins])
start = time.time()
print(CL.lnlike(parvals,parlist))#,priorvals,priorlist)
print("like call", time.time() - start)
print("prior",CL.lnprior(parvals,parlist,priorvals,priorlist))
#print "Prob",CL.lnprob(parvals,parlist,priorvals,priorlist)
start = time.time()
print(CL.lnlike(parvals2,parlist))#,priorvals,priorlist)
print("like call", time.time() - start)
print("prior",CL.lnprior(parvals2,parlist,priorvals,priorlist))
#print "Prob",CL.lnprob(parvals2,parlist,priorvals,priorlist)
