def P_Phi_p(fv, **cosmo): P_Phi_p = power_primordial(fv, **cosmo) loglog(cc.k, P_Phi_p, next(linecycler), linewidth = 2, label = '$k^{-2}$-factor: %s' %fv) legend(loc = 'best', prop = {'size':9}) axvline(x = 0.17 * 10**-2, color='r', ls = '--') # vertical line at k_eq axvline(x = 0.14 * 10**-1, color='b', ls = '-.') # vertical line at k_eq xlim((cc.k_min, cc.k_max)) xlabel('k') ylabel('$P^p_{Phi}$(k)')
def P_Phi(model, fv, **cosmo): P_Phi_p = power_primordial(fv, **cosmo); g_Delta_m = var.grwoth_factor(model,**cosmo)[0] g_Phi = var.grwoth_factor(model,**cosmo)[1] P_Delta = (g_Delta_m * g_Phi)**2 * P_Phi_p loglog(cc.k, P_Delta, next(linecycler), linewidth = 2, label = '$k^{-2}$-factor: %s' %fv) legend(loc = 'best', prop = {'size':9}) axvline(x = 0.17 * 10**-2, color='r', ls = '--') # vertical line at k_eq axvline(x = 0.14 * 10**-1, color='b', ls = '-.') # vertical line at k_eq xlim((cc.k_min, cc.k_max)) xlabel('k') ylabel('$P_m$(k, a = 1)')