def _mdelta_to_mdelta_via_m200(self, m_in, func, overdensity_in,
overdensity_out, z):
m200 = mass_concentration.mdelta_to_m200(m_in, func, overdensity_in,
(z, self._cosmo))
nfw = NFW(m200, func(m200, z, self._cosmo), z)
m_out = nfw.mass_Delta(overdensity_out)
return m_out
def test_mdelta_to_m200(self):
m_in = u.Quantity(2e14, u.solMass)
z = 0.2
func = mass_concentration.duffy_concentration
delta_in = 450
# consistency with mdelta_to_mdelta
md1 = mass_concentration.mdelta_to_m200(m_in, func, delta_in,
(z, self._cosmo))
md2 = mass_concentration.mdelta_to_mdelta(m_in, func,
delta_in, 200,
(z, self._cosmo))
assert_quantity_allclose(md1, md2)
# consistency with mass_Delta in NFW
nfw = NFW(md1, func(md1, z, self._cosmo), z)
m_out = nfw.mass_Delta(450)
assert_quantity_allclose(m_in, m_out)
