def test_mean_to_mean_nfw():
mdef = md.SOMean(Planck15, z=0, overdensity=200)
mdef2 = md.SOMean(Planck15, z=0, overdensity=300)
cduffy = mdef._duffy_concentration(1e12)
mnew, rnew, cnew = mdef.change_definition(1e12, mdef2)
mnew_, rnew_, cnew_ = changeMassDefinition(1e12, cduffy, 0, "200m", "300m",
"nfw")
assert np.isclose(mnew, mnew_, 1e-2)
assert np.isclose(rnew * 1e3, rnew_, 1e-2)
assert np.isclose(cnew, cnew_, 1e-2)
def test_from_colossus_name(colossus_cosmo):
assert md.from_colossus_name("200c") == md.SOCritical()
assert md.from_colossus_name("200m") == md.SOMean()
assert md.from_colossus_name("fof") == md.FOF()
assert md.from_colossus_name("800c") == md.SOCritical(overdensity=800)
assert md.from_colossus_name("vir") == md.SOVirial()
with pytest.raises(ValueError):
md.from_colossus_name("derp")
def test_mean_to_crit_nfw(colossus_cosmo):
mdef = md.SOMean(overdensity=200)
mdef2 = md.SOCritical(overdensity=300)
cduffy = mdef._duffy_concentration(1e12)
mnew, rnew, cnew = mdef.change_definition(1e12, mdef2)
mnew_, rnew_, cnew_ = changeMassDefinition(1e12, cduffy, 0, "200m", "300c",
"nfw")
assert np.isclose(mnew, mnew_, rtol=1e-2)
assert np.isclose(rnew * 1e3, rnew_, rtol=1e-2)
assert np.isclose(cnew, cnew_, rtol=1e-2)
def test_mean_to_vir_nfw(colossus_cosmo):
mdef = md.SOMean()
mdef2 = md.SOVirial()
cduffy = mdef._duffy_concentration(1e12)
mnew, rnew, cnew = mdef.change_definition(1e12, mdef2)
mnew_, rnew_, cnew_ = changeMassDefinition(1e12, cduffy, 0, "200m", "vir",
"nfw")
print(mnew, mnew_)
assert np.isclose(mnew, mnew_, rtol=1e-2)
assert np.isclose(rnew * 1e3, rnew_, rtol=1e-2)
assert np.isclose(cnew, cnew_, rtol=1e-2)
def test_colossus_name(colossus_cosmo):
assert md.SOMean().colossus_name == "200m"
assert md.SOCritical().colossus_name == "200c"
assert md.SOVirial().colossus_name == "vir"
assert md.FOF().colossus_name == "fof"