def ParseModel(model):
# if model=="LCDM_2":
if model == "LCDM" or "LCDM_" in model:
T = LCDMCosmology()
elif model == "LCDMmasslessnu" or "LCDMmasslessnu_" in model:
T = LCDMCosmology(mnu=0)
elif model == "nuLCDM" or "nuLCDM_" in model:
T = LCDMCosmology()
T.setVaryMnu()
elif model == "NnuLCDM" or "NnuLCDM_" in model:
LCDMCosmology.rd_approx = "tabulated_Nnu"
T = LCDMCosmology()
T.setVaryNnu()
elif model == "noradLCDM" or "noradLCDM_" in model:
T = LCDMCosmology(disable_radiation=True)
elif model == "oLCDM" or "oLCDM_" in model:
T = oLCDMCosmology()
elif model == "wCDM" or "wCDM_" in model:
T = wLCDMCosmology()
elif model == "waCDM" or "waCDM_" in model:
T = owa0CDMCosmology(varyOk=False)
elif model == "owCDM" or "owCDM_" in model:
T = owa0CDMCosmology(varywa=False)
elif model == "owaCDM" or "owaCDM_" in model:
T = owa0CDMCosmology()
elif model == "JordiCDM" or "JordiCDM_" in model:
T = JordiCDMCosmology()
elif model == "WeirdCDM" or "WeirdCDM_" in model:
T = WeirdCDMCosmology()
elif model == "TLight" or "TLight_" in model:
T = TLightCosmology()
elif model == "StepCDM" or "StepCDM_" in model:
T = StepCDMCosmology()
elif model == "Spline" or "Spline_" in model:
T = SplineLCDMCosmology()
elif model == "Decay" or "Decay_" in model:
T = DecayLCDMCosmology()
elif model == "PolyCDM" or "PolyCDM_" in model:
T = PolyCDMCosmology()
elif model == "EarlyDE" or "EarlyDE_" in model:
T = EarlyDECosmology()
else:
print "Cannot recognize model", model
sys.exit(1)
return T
def ParseModel(model):
if model=="LCDM":
T=LCDMCosmology()
elif model=="LCDMmasslessnu":
T=LCDMCosmology(mnu=0)
elif model=="nuLCDM":
T=LCDMCosmology()
T.setVaryMnu()
elif model=="NeffLCDM":
LCDMCosmology.rd_approx="CuestaNeff"
T=LCDMCosmology()
T.setVaryNnu()
elif model=="noradLCDM":
T=LCDMCosmology(disable_radiation=True)
elif model=="oLCDM":
T=oLCDMCosmology()
elif model=="nuoLCDM":
T=oLCDMCosmology()
T.setVaryMnu()
elif model=="wCDM":
T=wLCDMCosmology()
elif model=="nuwCDM":
T=wLCDMCosmology()
T.setVaryMnu()
elif model=="waCDM":
T=owa0CDMCosmology(varyOk=False)
elif model=="owCDM":
T=owa0CDMCosmology(varywa=False)
elif model=="owaCDM":
T=owa0CDMCosmology()
elif model=='FCDM':
T=OzcosCosmology()
elif model=="JordiCDM":
T=JordiCDMCosmology()
elif model=="WeirdCDM":
T=WeirdCDMCosmology()
elif model=="TLight":
T=TiredLightDecorator(PolyCDMCosmology())
elif model=="StepCDM":
T=StepCDMCosmology()
elif model=="Spline":
T=SplineLCDMCosmology()
elif model=="DecayFrac":
T=DecayLCDMCosmology()
elif model=="Decay":
T=DecayLCDMCosmology(varyxfrac=False)
elif model=="Decay01":
T=DecayLCDMCosmology(varyxfrac=False,xfrac=0.1)
elif model=="Decay05":
T=DecayLCDMCosmology(varyxfrac=False,xfrac=0.5)
elif model=="PolyCDM":
T=PolyCDMCosmology()
elif model=="fPolyCDM":
T=PolyCDMCosmology(varyOk=False)
elif model=="EarlyDE":
T=EarlyDECosmology(varyw=False,userd_DE=False)
elif model=="EarlyDE_rd_DE":
T=EarlyDECosmology(varyw=False)
elif model=="SlowRDE":
T=SlowRDECosmology(varyOk=False)
elif model=="Quint":
T=QuintCosmology()
elif model=='wDM':
T=wDMCosmology()
else:
print "Cannot recognize model", model
sys.exit(1)
return T
def ParseModel(model):
if model == "LCDM":
T = LCDMCosmology()
elif model == "LCDMmasslessnu":
T = LCDMCosmology(mnu=0)
elif model == "nuLCDM":
T = LCDMCosmology()
T.setVaryMnu()
elif model == "NeffLCDM":
LCDMCosmology.rd_approx = "CuestaNeff"
T = LCDMCosmology()
T.setVaryNnu()
elif model == "noradLCDM":
T = LCDMCosmology(disable_radiation=True)
elif model == "oLCDM":
T = oLCDMCosmology()
elif model == "nuoLCDM":
T = oLCDMCosmology()
T.setVaryMnu()
elif model == "wCDM":
T = wLCDMCosmology()
elif model == "nuwCDM":
T = wLCDMCosmology()
T.setVaryMnu()
elif model == "waCDM":
T = owa0CDMCosmology(varyOk=False)
elif model == "owCDM":
T = owa0CDMCosmology(varywa=False)
elif model == "owaCDM":
T = owa0CDMCosmology()
elif model == "JordiCDM":
T = JordiCDMCosmology()
elif model == "WeirdCDM":
T = WeirdCDMCosmology()
elif model == "TLight":
T = TiredLightDecorator(PolyCDMCosmology())
elif model == "StepCDM":
T = StepCDMCosmology()
elif model == "Spline":
T = SplineLCDMCosmology()
elif model == "DecayFrac":
T = DecayLCDMCosmology()
elif model == "Decay":
T = DecayLCDMCosmology(varyxfrac=False, xfrac=1.0)
elif model == "Decay01":
T = DecayLCDMCosmology(varyxfrac=False, xfrac=0.1)
elif model == "Decay05":
T = DecayLCDMCosmology(varyxfrac=False, xfrac=0.5)
elif model == "PolyCDM":
T = PolyCDMCosmology()
elif model == "fPolyCDM":
T = PolyCDMCosmology(varyOk=False)
elif model == "EarlyDE":
T = EarlyDECosmology(varyw=False, userd_DE=False)
elif model == "EarlyDE_rd_DE":
T = EarlyDECosmology(varyw=False)
elif model == "SlowRDE":
T = SlowRDECosmology(varyOk=False)
elif model == "Quint_last":
T = QuintCosmology()
elif model == 'wDM':
T = wDMCosmology()
else:
print(("Cannot recognize model", model))
sys.exit(1)
return T
def ParseModel(model):
if model == "LCDM":
T = LCDMCosmology()
elif model == "LCDMmasslessnu":
T = LCDMCosmology(mnu=0)
elif model == "nuLCDM":
T = LCDMCosmology()
T.setVaryMnu()
elif model == "NeffLCDM":
LCDMCosmology.rd_approx = "CuestaNeff"
T = LCDMCosmology()
T.setVaryNnu()
elif model == "noradLCDM":
T = LCDMCosmology(disable_radiation=True)
elif model == "oLCDM":
T = oLCDMCosmology()
elif model == "nuoLCDM":
T = oLCDMCosmology()
T.setVaryMnu()
elif model == "wCDM":
T = wLCDMCosmology()
elif model == "nuwCDM":
T = wLCDMCosmology()
T.setVaryMnu()
elif model == "waCDM":
T = owa0CDMCosmology(varyOk=False)
elif model == "owCDM":
T = owa0CDMCosmology(varywa=False)
elif model == "owaCDM":
T = Lowa0CDMCosmology()
elif model == 'FCDM':
T = OzcosCosmology()
elif model == 'HornFco_c_cubic':
T = HornFcoCosmology(varyc=True)
elif model == 'HornFco_f_cubic':
T = HornFcoCosmology(varyf=True)
elif model == "JordiCDM":
T = JordiCDMCosmology()
elif model == "WeirdCDM":
T = WeirdCDMCosmology()
elif model == "TLight":
T = TiredLightDecorator(PolyCDMCosmology())
elif model == "StepCDM":
T = StepCDMCosmology()
elif model == "Spline":
T = SplineLCDMCosmology()
elif model == "DecayFrac":
T = DecayLCDMCosmology()
elif model == "Decay":
T = DecayLCDMCosmology(varyxfrac=False)
elif model == "Decay01":
T = DecayLCDMCosmology(varyxfrac=False, xfrac=0.1)
elif model == "Decay05":
T = DecayLCDMCosmology(varyxfrac=False, xfrac=0.5)
elif model == "PolyCDM":
T = PolyCDMCosmology()
elif model == "fPolyCDM":
T = PolyCDMCosmology(varyOk=False)
elif model == "EarlyDE":
T = EarlyDECosmology(varyw=False, userd_DE=False)
elif model == "EarlyDE_rd_DE":
T = EarlyDECosmology(varyw=False)
elif model == "SlowRDE":
T = SlowRDECosmology(varyOk=False)
elif model == "Quint":
T = QuintCosmology_try()
elif model == "Quintess":
T = QuintomCosmology(varymquin=True)
elif model == "Phantom":
T = QuintomCosmology(varymphan=True)
elif model == "Quintom":
T = QuintomCosmology(varymquin=True, varymphan=True, varyiniphi=True)
elif model == "Anisotropic":
T = AnisotropicCosmology()
elif model == "Fourier":
T = FourierCosmology()
elif model == 'wDM':
T = wDMCosmology()
else:
print "Cannot recognize model", model
sys.exit(1)
return T
def ParseModel(model):
if model == "LCDM":
T = LCDMCosmology()
elif model == "LCDMmasslessnu":
T = LCDMCosmology(mnu=0)
elif model == "nuLCDM":
T = LCDMCosmology()
T.setVaryMnu()
elif model == "NeffLCDM":
LCDMCosmology.rd_approx = "CuestaNeff"
T = LCDMCosmology()
T.setVaryNnu()
elif model == "noradLCDM":
T = LCDMCosmology(disable_radiation=True)
elif model == "oLCDM":
T = oLCDMCosmology()
elif model == "nuoLCDM":
T = oLCDMCosmology()
T.setVaryMnu()
elif model == "wCDM":
T = wLCDMCosmology()
elif model == "nuwCDM":
T = wLCDMCosmology()
T.setVaryMnu()
elif model == "waCDM":
T = owa0CDMCosmology(varyOk=False)
elif model == "owCDM":
T = owa0CDMCosmology(varywa=False)
elif model == "owaCDM":
T = owa0CDMCosmology()
elif model == "JordiCDM":
T = JordiCDMCosmology()
elif model == "WeirdCDM":
T = WeirdCDMCosmology()
elif model == "TLight":
T = TiredLightDecorator(PolyCDMCosmology())
elif model == "StepCDM":
T = StepCDMCosmology()
elif model == "Spline":
T = SplineLCDMCosmology()
elif model == "DecayFrac":
T = DecayLCDMCosmology()
elif model == "Decay":
T = DecayLCDMCosmology(varyxfrac=False, xfrac=1.0)
elif model == "Decay01":
T = DecayLCDMCosmology(varyxfrac=False, xfrac=0.1)
elif model == "Decay05":
T = DecayLCDMCosmology(varyxfrac=False, xfrac=0.5)
elif model == "PolyCDM":
T = PolyCDMCosmology()
elif model == "fPolyCDM":
T = PolyCDMCosmology(polyvary=['Om1','Om2'])
elif model == "PolyOk": ## polycdm for OK
T = PolyCDMCosmology(Ok_prior=10.)
elif model == "PolyOkc": ## polycdm sans Om2 term to couple two
T = PolyCDMCosmology(polyvary=['Om1','Ok'],Ok_prior=10.)
elif model == "PolyOkf": ## polycdm sans Om2 term to couple two
T = PolyCDMCosmology(polyvary=['Om1','Om2'])
elif model == "EarlyDE":
T = EarlyDECosmology(varyw=False, userd_DE=False)
elif model == "EarlyDE_rd_DE":
T = EarlyDECosmology(varyw=False)
elif model == "SlowRDE":
T = SlowRDECosmology(varyOk=False)
elif model == "Quint_last":
T = QuintCosmology()
elif model == 'wDM':
T = wDMCosmology()
else:
print("Cannot recognize model", model)
sys.exit(1)
return T
