def plot_Oede(i,
param,
ylabel=None,
ytick=None,
ymax=1.01,
quint=False,
ylim=None):
ax = fig.add_subplot(4, 3, i)
for i in range(rr):
sc = ax.scatter(Oede_cmb[i], param[i], color=Usf.colour(i + 1))
ax.grid(True)
if quint:
plt.xticks(np.arange(0.02, 0.0601, 0.01))
plt.xlim([0.02, 0.0601])
else:
plt.xticks(np.arange(0.0, 0.09, 0.02))
plt.xlim(xmin=0)
if ylim is not None:
yin, yf = ylim
plt.ylim([yin, yf])
if ytick is not None:
yin, yf, yd = ytick
plt.yticks(np.arange(yin, yf, yd))
plt.xlabel("$ \Omega_{ede}(z_{drag})$")
plt.ylabel(ylabel)
def plot_Oede(i, param, ylabel= None, ytick=None, ymax=1.01, quint= False, ylim=None):
ax = fig.add_subplot(4,3,i)
for i in range(rr):
sc = ax.scatter(Oede_cmb[i], param[i], color=Usf.colour(i+1))
ax.grid(True)
if quint:
plt.xticks(np.arange(0.02, 0.0601, 0.01))
plt.xlim([0.02,0.0601])
else:
plt.xticks(np.arange(0.0, 0.09, 0.02))
plt.xlim(xmin=0)
if ylim is not None:
yin, yf = ylim
plt.ylim([yin, yf])
if ytick is not None:
yin, yf, yd= ytick
plt.yticks(np.arange(yin, yf, yd))
plt.xlabel("$ \Omega_{ede}(z_{drag})$")
plt.ylabel(ylabel)
if ytick is not None:
yin, yf, yd= ytick
plt.yticks(np.arange(yin, yf, yd))
plt.xlabel("$ \Omega_{ede}(z_{drag})$")
plt.ylabel(ylabel)
#---------------Plot Files----------------------#
if True:
fig = pylab.figure(figsize=(10,12))
ax8 = fig.add_subplot(4, 3, 1)
for i in range(rr):
z = np.exp(-lna[i])
ax8.plot(z, Oede[i], color=Usf.colour(i+1))
ax8.grid(True)
ax8.set_xscale('log')
plt.xlabel("$z+1$")
plt.yticks(np.arange(0,1,0.2))
plt.ylabel("$\\Omega_{ede}$")
plt.xlim([1,1.0e6])
ax9 = fig.add_subplot(4, 3, 2)
for i in range(rr):
z = np.exp(-lna[i])
ax9.plot(z, wede[i], color=Usf.colour(i+1),
label = '$\Omega_{ede} = $%1.3f'%(Oede_cmb[i]))
ax9.legend(loc='center left', bbox_to_anchor=(1.5, 0.5))
ax9.grid(True)
T.updateParams([lam_, V0_, A_, lB_])
sol = T.Ini_phi()
xx, yy = sol.T
pts.append(T.O_phi(lna))
pts3.append(T.w_ede(lna))
Oma.append(Om0/(exp(lna)**3)*(1.0/T.hub((lna), sol.T))**2)
Ora.append(Or0/(exp(lna)**4)*(1.0/T.hub((lna), sol.T))**2)
pts4.append(xx)
if True:
fig = plt.figure(figsize=(16,5))
ax=fig.add_subplot(1,3,1)
for i in range(4):
ax.plot(z, pts[i], color = Usf.colour(i+1))
ax.plot(z, Oma[i], 'g--')
ax.plot(z, Ora[i], 'y--')
plt.ylim([-0.01,1.01])
ax.set_xscale('log')
plt.xlim([1.e-1, 1e8])
plt.xlabel("$z$")
plt.ylabel("$\\Omega_{ede}$")
ax.grid(True)
ax2=fig.add_subplot(1,3,2)
for i in range(4):
ax2.plot(z, pts3[i], color = Usf.colour(i+1))
plt.ylim(ymin=-1.01)
ax2.set_xscale('log')
plt.xlim([1.e-1, 1e8])
if ytick is not None:
yin, yf, yd = ytick
plt.yticks(np.arange(yin, yf, yd))
plt.xlabel("$ \Omega_{ede}(z_{drag})$")
plt.ylabel(ylabel)
#---------------Plot Files----------------------#
if True:
fig = pylab.figure(figsize=(10, 12))
ax8 = fig.add_subplot(4, 3, 1)
for i in range(rr):
z = np.exp(-lna[i])
ax8.plot(z, Oede[i], color=Usf.colour(i + 1))
ax8.grid(True)
ax8.set_xscale('log')
plt.xlabel("$z+1$")
plt.yticks(np.arange(0, 1, 0.2))
plt.ylabel("$\\Omega_{ede}$")
plt.xlim([1, 1.0e6])
ax9 = fig.add_subplot(4, 3, 2)
for i in range(rr):
z = np.exp(-lna[i])
ax9.plot(z,
wede[i],
color=Usf.colour(i + 1),
label='$\Omega_{ede} = $%1.3f' % (Oede_cmb[i]))
ax9.legend(loc='center left', bbox_to_anchor=(1.5, 0.5))
T.updateParams([lam_, V0_, A_, lB_])
sol = T.Ini_phi()
xx, yy = sol.T
pts.append(T.O_phi(lna))
pts3.append(T.w_ede(lna))
Oma.append(Om0 / (exp(lna)**3) * (1.0 / T.hub((lna), sol.T))**2)
Ora.append(Or0 / (exp(lna)**4) * (1.0 / T.hub((lna), sol.T))**2)
pts4.append(xx)
if True:
fig = plt.figure(figsize=(16, 5))
ax = fig.add_subplot(1, 3, 1)
for i in range(4):
ax.plot(z, pts[i], color=Usf.colour(i + 1))
ax.plot(z, Oma[i], 'g--')
ax.plot(z, Ora[i], 'y--')
plt.ylim([-0.01, 1.01])
ax.set_xscale('log')
plt.xlim([1.e-1, 1e8])
plt.xlabel("$z$")
plt.ylabel("$\\Omega_{ede}$")
ax.grid(True)
ax2 = fig.add_subplot(1, 3, 2)
for i in range(4):
ax2.plot(z, pts3[i], color=Usf.colour(i + 1))
plt.ylim(ymin=-1.01)
ax2.set_xscale('log')
plt.xlim([1.e-1, 1e8])
