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)')
def g_Delta_m(model, BP, **cosmo): g_Delta_m = var.grwoth_factor(model, BP,**cosmo)[0] loglog(cc.k, g_Delta_m, next(linecycler), linewidth = 2, label = model) legend(loc = 'best', prop = {'size':9}) xlim((0.3 * 10**-3, 0.3)) xlabel('k') ylabel('$g_m$(k, a = 1)')