Closed
Description
The old version of fastMP had this code:
# Re-estimate centers using the selected stars
if len(st_idx) > self.N_clust_min:
xy_c, vpd_c, plx_c = self.get_5D_center(
lon[st_idx], lat[st_idx], pmRA[st_idx], pmDE[st_idx],
plx[st_idx])
idx_selected += list(st_idx)
N_runs += 1
which stored the probable members in the idx_selected list.
By default N_clust_min=25 and the estimate_nmembs function returned N_survived=25 as the minimum value. This meant that if the estimate_nmembs function returned N_survived=25 no stars would be stored in idx_selected and the probabilities would all be 0.
There are 3531 entries with N_50=0: entries_N50_0.txt