Following on from #77 and #80, sampling positions for multiple populations now returns lists of individual arrays. The functions in the galaxies module should be update to allow easy composition.
Assuming ngal_arr is an array of shape (N,), the galaxy counts gal_cnt would be of shape (N,), and this should be broadcast against potential leading dimensions of dndz to produce a list of gal_z samples.
The existing gal_pop output, which has always been a little awkward, will then no longer be necessary.
Following on from #77 and #80, sampling positions for multiple populations now returns lists of individual arrays. The functions in the
galaxiesmodule should be update to allow easy composition.Assuming
ngal_arris an array of shape(N,), the galaxy countsgal_cntwould be of shape(N,), and this should be broadcast against potential leading dimensions ofdndzto produce a list ofgal_zsamples.The existing
gal_popoutput, which has always been a little awkward, will then no longer be necessary.