[subset] use add_array to populate feature_indexes. This is much faster then calling a bunch of individual add()'s.