virial_radius_analytic

cosmic.sample.cmc.elson.virial_radius_analytic(gamma, r_max)

Virial radius is best calculated directly, since rmax may be pretty far from infinity. Directly integrate 4*pi*r*rho*m_enclosed from 0 to rMax to get the binding energy, then just divide 0.5 by that.