def delta_bias(model, f_NL, **cosmo): delta_bias = var.galaxy_bias(model, f_NL, **cosmo )[1] semilogx(cc.k, delta_bias, next(linecycler), linewidth = 2, label = model + "-$f_{NL}$_%s" %f_NL) legend(loc = 'best') xlim((0.3 * 10**-3, 0.3* 10**-2)) xlabel('k') ylabel('$\Delta b_g$(k, a = 1)')
def bias(model, f_NL, **cosmo): bias = var.galaxy_bias(model, f_NL, **cosmo)[0] loglog(cc.k, bias, next(linecycler), linewidth = 2, label = model+ "-$f_{NL}$_%s" %f_NL) legend(loc = 'best') xlim((0.3 * 10**-3, 0.3* 10**-2)) xlabel('k') ylabel('$b_g$(k, a = 1)')