Functions
def	st_gs (sigma, dc)

def	st_b1 (sigma, dc)

def	duffy_c (m, z)

def	get_halomod_pk (cp, karr, z, integrate=True, fname_out=None)

Variables
dictionary	cpar1

dictionary	cpar2

dictionary	cpar3

bool	USE_INT = False

int	LKHMIN = -4

int	LKHMAX = 1

int	NKH = 256

int	NM = 2048

float	LM0 = 5.5

int	LMMIN = 6

int	LMMAX = 16

dictionary	params

	integrate

	fname_out

Function Documentation

def halomod_bm.duffy_c	(	m,
		z
	)

Definition at line 48 of file halomod_bm.py.

Referenced by get_halomod_pk().

 def duffy_c(m,z) :
     return params['DUF_A']*(m/2E12)**params['DUF_B']*(1+z)**params['DUF_C']
 

def halomod_bm.get_halomod_pk	(	cp,
		karr,
		z,
		integrate = `True`,
		fname_out = `None`
	)

Definition at line 51 of file halomod_bm.py.

References duffy_c(), st_b1(), and st_gs().

 def get_halomod_pk(cp,karr,z,integrate=True,fname_out=None) :
     cosmo=ccl.Cosmology(Omega_c=cp['Om_m']-cp['Om_b'],Omega_b=cp['Om_b'],h=cp['h'],
                         sigma8=cp['sig8'],n_s=cp['n'],
                         transfer_function='eisenstein_hu',matter_power_spectrum='linear')
 
     lmarr=LMMIN+(LMMAX-LMMIN)*(np.arange(NM)+0.5)/NM #log(M*h/M_sum)
     dlm=(LMMAX-LMMIN)/NM
     lmarrb=np.zeros(NM+2); lmarrb[0]=lmarr[0]-dlm; lmarrb[1:-1]=lmarr; lmarrb[-1]=lmarr[-1]+dlm
     marr=10.**lmarr #M (in Msun/h) M_h Ms/h = M Ms
     marrb=10.**lmarrb #M (in Msun/h) M_h Ms/h = M Ms
     sigmarrb=ccl.sigmaM(cosmo,marrb/cp['h'],1./(1+z))
     sigmarr=sigmarrb[1:-1]
     dlsMdl10M=np.fabs((sigmarrb[2:]-sigmarrb[:-2])/(2*dlm))/sigmarr
     omega_z=cp['Om_m']*(1+z)**3/(cp['Om_m']*(1+z)**3+1-cp['Om_m'])
     delta_c=params['DELTA_C']*(1+0.012299*np.log10(omega_z))
     Delta_v=(18*np.pi**2+82*(omega_z-1)-39*(omega_z-1)**2)/omega_z
     fM=st_gs(sigmarr,delta_c)*dlsMdl10M
     bM=st_b1(sigmarr,delta_c)
     cM=duffy_c(marr,z)
     cMf=interp1d(lmarr,cM,kind='linear',bounds_error=False,fill_value=cM[0])
     rhoM=ccl.rho_x(cosmo,1.,'matter')/cp['h']**2
 
     #Compute mass function normalization
     fM0=1-np.sum(fM)*dlm
     fbM0=1-np.sum(fM*bM)*dlm
 
     rvM=(3*marr/(4*np.pi*rhoM*Delta_v))**0.333333333333
     rsM=rvM/cM
     rsM0=(3*10.**LM0/(4*np.pi*rhoM*Delta_v))**0.333333333/duffy_c(10.**LM0,z)
     rsMf=interp1d(lmarr,rsM,kind='linear',bounds_error=False,fill_value=rsM0)
 
     f1hf=interp1d(lmarr,marr*fM/rhoM,kind='linear',bounds_error=False,fill_value=0)
     f1h0=fM0*10.**LM0/rhoM
     f2hf=interp1d(lmarr,fM*bM,kind='linear',bounds_error=False,fill_value=0)
     f2h0=fbM0
 
     #NFW profile
     def u_nfw(lm,k) :
         x=k*rsMf(lm) #k*r_s
         c=cMf(lm)
         sic,cic=sici(x*(1+c))
         six,cix=sici(x)
         fsin=np.sin(x)*(sic-six)
         fcos=np.cos(x)*(cic-cix)
         fsic=np.sin(c*x)/(x*(1+c))
         fcon=np.log(1.+c)-c/(1+c)
 
         return (fsin+fcos-fsic)/fcon
 
     def p1h(k) :
         def integ(lm) :
             u=u_nfw(lm,k)
             f1h=f1hf(lm)
             return u*u*f1h
         if integrate :
             return quad(integ,lmarr[0],lmarr[-1])[0]+f1h0*(u_nfw(LM0,k))**2
         else :
             return np.sum(integ(lmarr))*dlm+f1h0*(u_nfw(LM0,k))**2
     
     def b2h(k) :
         def integ(lm) :
             u=u_nfw(lm,k)
             f2h=f2hf(lm)
             return u*f2h
         if integrate :
             return quad(integ,lmarr[0],lmarr[-1])[0]+f2h0*u_nfw(LM0,k)
         else :
             return np.sum(integ(lmarr))*dlm+f2h0*u_nfw(LM0,k)
 
 
     p1harr=np.array([p1h(kk) for kk in karr])
     b2harr=np.array([b2h(kk) for kk in karr])
     pklin=ccl.linear_matter_power(cosmo,karr*cp['h'],1./(1+z))*cp['h']**3
 
     p2harr=pklin*b2harr**2
     ptarr=p1harr+p2harr
     np.savetxt(fname_out,np.transpose([karr,pklin,p2harr,p1harr,ptarr]))
     return pklin,pklin*b2harr**2,p1harr,pklin*b2harr**2+p1harr
 
 k_arr=10.**(LKHMIN+(LKHMAX-LKHMIN)*(np.arange(NKH)+0.5)/NKH)

def halomod_bm.st_b1	(	sigma,
		dc
	)

Definition at line 42 of file halomod_bm.py.

Referenced by get_halomod_pk().

 def st_b1(sigma,dc) :
     nu=dc/sigma
     anu2=params['ST_a']*nu*nu
 
     return 1+(anu2-1+2*params['ST_p']/(1+anu2**params['ST_p']))/params['DELTA_C']
 

def halomod_bm.st_gs	(	sigma,
		dc
	)

Definition at line 35 of file halomod_bm.py.

Referenced by get_halomod_pk().

 def st_gs(sigma,dc) :
     nu=dc/sigma
     anu2=params['ST_a']*nu*nu
     expnu=np.exp(-anu2*0.5)
 
     return params['ST_A']*np.sqrt(2*params['ST_a']/np.pi)*(1+(1./anu2)**params['ST_p'])*nu*expnu
 

Variable Documentation

dictionary halomod_bm.cpar1

Initial value:

1 = {'Om_m':0.3,'Om_b':0.05,'Om_nu':0.,

2 'h':0.70,'sig8':0.8,'n':0.96}

Definition at line 9 of file halomod_bm.py.

dictionary halomod_bm.cpar2

Initial value:

1 = {'Om_m':0.272,'Om_b':0.0455,'Om_nu':0.,

2 'h':0.704,'sig8':0.81,'n':0.967}

Definition at line 11 of file halomod_bm.py.

dictionary halomod_bm.cpar3

Initial value:

1 = {'Om_m':0.3175,'Om_b':0.0490,'Om_nu':0.,

2 'h':0.6711,'sig8':0.834,'n':0.9624}

Definition at line 13 of file halomod_bm.py.

halomod_bm.fname_out

Definition at line 133 of file halomod_bm.py.

halomod_bm.integrate

Definition at line 133 of file halomod_bm.py.

Referenced by App_Var< T >::Impl< T >.integration().

int halomod_bm.LKHMAX = 1

Definition at line 19 of file halomod_bm.py.

int halomod_bm.LKHMIN = -4

Definition at line 18 of file halomod_bm.py.

float halomod_bm.LM0 = 5.5

Definition at line 23 of file halomod_bm.py.

int halomod_bm.LMMAX = 16

Definition at line 25 of file halomod_bm.py.

int halomod_bm.LMMIN = 6

Definition at line 24 of file halomod_bm.py.

int halomod_bm.NKH = 256

Definition at line 20 of file halomod_bm.py.

int halomod_bm.NM = 2048

Definition at line 22 of file halomod_bm.py.

Referenced by pwr_spec_k(), pwr_spec_k_init(), and vel_pwr_spec_k().

dictionary halomod_bm.params

Initial value:

 = {'DELTA_C':1.686470199841,
         'ST_A': 0.322,
         'ST_a': 0.707,
         'ST_p': 0.3,
         'DUF_A':7.85,
         'DUF_B':-0.081,
         'DUF_C':-0.71}

Definition at line 27 of file halomod_bm.py.

Referenced by __attribute__(), __ctest_parameters_create_general_nu_list_run(), __ctest_parameters_create_general_nu_sum_equal_run(), __ctest_parameters_create_general_nu_sum_inverted_run(), __ctest_parameters_create_general_nu_sum_run(), __ctest_parameters_read_write_run(), ccl_cosmology_compute_growth(), ccl_cosmology_create(), ccl_cosmology_create_flat_lcdm(), ccl_create_dNdz_info(), ccl_create_photoz_info(), ccl_parameters_create(), ccl_parameters_create_flat_lcdm(), ccl_parameters_read_yaml(), compare_bbks(), compare_bcm(), compare_cls(), compare_corr(), compare_correlation_3d(), compare_distances(), compare_distances_hiz(), compare_distances_hiz_mnu(), compare_distances_mnu(), compare_eh(), compare_emu(), compare_emu_nu(), compare_growth(), compare_growth_hiz(), compare_halomod(), compare_massfunc(), compare_power_nu(), compare_power_nu_nl(), compare_sigmam(), growth_ode_system(), main(), anonymous_namespace{core_power.cpp}::Integr_obj_qag.operator()(), anonymous_namespace{core_power.cpp}::Integr_obj_qagiu.operator()(), anonymous_namespace{core_power.cpp}::Integr_obj_qawo.operator()(), anonymous_namespace{core_power.cpp}::Integr_obj_qawf.operator()(), and write_Pk_CLASS().

halomod_bm.USE_INT = False

Definition at line 16 of file halomod_bm.py.

Functions

Variables

Function Documentation

Variable Documentation