maybe a way to check the subset treatment if `efficient` is True is to use repeated subsets in the original data (`np.tile`), or try n_sub = 1 (in default).  But I'm not sure it doesn't break if there is no variation across subsets, or there is only one.