Add CMBlensing mg factors
@@ -814,9 +814,10 @@ static double transfer_cmblens(int l,double k,ccl_cosmology *cosmo,CCL_ClTracer
double w=1-chi/clt->chi_source;
// If muSigma parameterisation of gravity is in effect and
// Sigma0>0, add the relevant factor here.
- //if (fabs(cosmo->params.sigma_0)>1e-15){
- // w = w * (1. + ccl_Sig_MG(cosmo,ccl_scale_factor_of_chi(cosmo,chi, status), status));
- //return clt->prefac_lensing*l*(l+1.)*w/(a*chi*k*k);
+ if (fabs(cosmo->params.sigma_0)>1e-15){
+ w = w * (1. + ccl_Sig_MG(cosmo,ccl_scale_factor_of_chi(cosmo,chi, status), status));
+ }
+ return clt->prefac_lensing*l*(l+1.)*w/(a*chi*k*k);
}
return 0;
}
