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)