This breaks this function completely. The whole key invariant is that
if `data, indices, filenames = featurize_all(featurizer, trajectories)`, then

`data[i] = featurize(md.load(filenames[i])[indices[i]])`

but this breaks that.