def __init__(self,
varymquin=False,
varymphan=False,
varyiniphi=False,
varybeta=False):
## two parameters: Om and h
"""Is better to start the chains at masses equal one, othewise
may take much longer"""
self.varymquin = varymquin
self.varymphan = varymphan
self.varybeta = varybeta
self.varyiniphi = varyiniphi
self.mquin = 0 if (varymphan and (not varymquin)) else mquin_par.value
self.mphan = 0 if (varymquin and (not varymphan)) else mphan_par.value
self.beta = 0 if (not varybeta) else beta_par.value
self.iniphi = iniphi_par.value
LCDMCosmology.__init__(self, mnu=0)
self.lna = np.linspace(-10, 0, 500)
self.z = np.exp(-self.lna) - 1.
self.zvals = np.linspace(0, 3, 100)
self.cte = 3.0 * self.h**2
self.n = 2.
self.m = 2
self.chatty = True
self.qp = 1
self.updateParams([])
def __init__(self,
varyw=False,
varywa=False,
varyOk=False,
varyq0=True,
varyj0=True,
varyOmega0=True):
# three parameters: w, wa, Ok
self.varyw = varyw
self.varywa = varywa
self.varyOk = varyOk
self.varyq0 = varyq0
self.varyj0 = varyj0
self.varyOmega0 = varyOmega0
self.Ok = Ok_par.value
self.w0 = w_par.value
self.wa = wa_par.value
self.q0 = q0.value
self.j0 = j0.value
self.Omega0 = Omega0.value
LCDMCosmology.__init__(self)
def __init__(self, varyalpha=True, varyOk=False):
## two parameters: Om and h
"""Is better to start the chains at masses equal one, othewise
may take much longer"""
self.varyalpha = varyalpha
self.varyOk = varyOk
self.alpha = alpha_par.value
self.Ok = Ok_par.value
#Compute Vcap either analytically or find it
self.fixVcap = False
#Pick one for initial conditions
#self.inic = 'Matter'
#self.a_in = 5e-3
self.inic = 'Radiation'
self.a_in = 5e-5
self.ndots = 1000
self.zvals = np.linspace(0, 4, 50)
self.ascal = np.linspace(self.a_in, 1, self.ndots)
LCDMCosmology.__init__(self, mnu=0)
self.updateParams([])
def __init__(self, varyalpha=False, varybeta=False, varyilam=False,\
varymu=False, varyOk=False,
alpha=1, beta=1, mu=1, ilam=1, eps=1, curv=0):
"""Is better to start the chains at masses equal one, othewise
may take much longer"""
#eps => 1=Quintes, -1=Panthom
self.varyOk = varyOk
self.varymu = varymu
self.varyilam = varyilam
self.varybeta = varybeta
self.varyalpha = varyalpha
self.Ok = curv
self.alpha = alpha
self.beta = beta
self.mu = mu
self.ilam = ilam
self.eps = eps
self.lna = np.linspace(-6, 0, 500)
self.z = np.exp(-self.lna) - 1.
self.zvals = np.linspace(0, 5, 200)
LCDMCosmology.__init__(self, mnu=0)
self.updateParams([])
def __init__(self, varyw=True, varyOk=True):
# two parameters: Om and h
self.varyw = varyw
self.varyOk = varyOk
self.Ok = Ok_par.value
self.dw0 = dw_par.value
LCDMCosmology.__init__(self)
def __init__(self, varyw=True, varywa=True, varyOk=True):
# three parameters: w, wa, Ok
self.varyw = varyw
self.varywa = varywa
self.varyOk = varyOk
self.Ok = Ok_par.value
self.w0 = w_par.value
self.wa = wa_par.value
LCDMCosmology.__init__(self)
def __init__(self, polyvary=['Om1', 'Om2', 'Ok'], Ok_prior=0.1):
# Ok, LCDM has Omega_m, we also have Omega_1 and Omega_2
# and Lambda is then what remains
##
self.Ok = Ok_par.value
self.Om1 = Om1_par.value
self.Om2 = Om2_par.value
self.varyOm1 = 'Om1' in polyvary
self.varyOm2 = 'Om2' in polyvary
self.varyOk = 'Ok' in polyvary
self.Ok_prior = Ok_prior
LCDMCosmology.__init__(self, mnu=0)
def __init__(self, zeroDE=False, kwargs_LCDM={}):
# two parameters: Om and h
self.Ok = Ok_par.value
self.zeroDE = zeroDE
LCDMCosmology.__init__(self, **kwargs_LCDM)