draw_r_vr_vt

cosmic.sample.cmc.elson.draw_r_vr_vt(N=100000, r_max=300, gamma=4)

Draw random velocities and positions from the Elson profile.

N = number of stars r_max = maximum number of virial radii for the farthest star gamma = steepness function of the profile

Note that gamma=4 is the Plummer profile

returns (vr,vt,r) in G=M_cluster=1 units

N.B. if you’re getting a “Expect x to not have duplicates” error with a large gamma, use a smaller r_max.