virial_radius_numerical¶ cosmic.sample.cmc.king.virial_radius_numerical(r, rho_r, M_enclosed)[source]¶ Virial radius is best calculated directly. Directly integrate 4*pi*r*rho*m_enclosed over the samples binding energy, then just divide 0.5 by that.