Omegav = 0.
Omegak = 1. - Omegam - Omegav
print 'Omegab = ',omegab/h2,'Omegam = ',(omegac+omegab)/h2,'Omegav = ',Omegav, 'Omegak = ', Omegak
params = {'H0':H0,'omegab':Omegab,'omegac':Omegac,'omegak':Omegak,'scalar_index':0.9624,
'reion__use_optical_depth':True,'reion__optical_depth':0.0925}
nz = 100
zmin = 0
zmax = 1100
# BUILD IT ###############################################################################
#pklib = ln.build_pklib(params, zmin, zmax, nz, extra_name='nolambda_',log=True)
# RESTORE IT ##############################################################################
pklib = ln.read_pklib(zmin, zmax, nz, extra_name = 'nolambda_')
def profile(x,y,range=None,nbins=10,fmt=None,plot=True, dispersion=True, median=False):
if range == None:
mini = np.min(x)
maxi = np.max(x)
else:
mini = range[0]
maxi = range[1]
clf()
subplot(2,1,1)
errorbar(nbtry, allmeans, yerr=allerrmeans)
xscale('log')
subplot(2,1,2)
errorbar(nbtry, allvars, yerr=allerrvars)
xscale('log')
#### P(K) library at various redshifts
nz = 100
zmin = 0
zmax = 1100
# BUILD IT ###############################################################################
#pklib = ln.build_pklib(params, zmin, zmax, nz)
# RESTORE IT ##############################################################################
pklib = ln.read_pklib(zmin, zmax, nz)
#### Skewers library with the same seed at different reshifts
# set the x range so that it extends well beyond the distance at maximum redshift
zvals = np.linspace(0,zmax,1000)
lcdm = cd.set_omega_k_0({'omega_M_0' : 0.3, 'omega_lambda_0' : 0.7, 'h' : 0.7})
d_lcdm = cd.comoving_distance_transverse(zvals, **lcdm)
maxd = np.max(d_lcdm)*1.2
theseed = 1
nn = 2**18
lightcone = ln.lightcone_1d(maxd, nn, pklib, seed=theseed)
thezinit = linspace(0,zmax, 10000)
dinit = cd.comoving_distance_transverse(thezinit, **lcdm)
