get_velocities¶

cosmic.sample.cmc.king.get_velocities(r, r_profile, psi_profile, M_enclosed_profile)[source]¶

The correct way to generate velocities: start from the distribution function and use rejection sampling.

returns (vr,vt) with same length as r