docs: def find(self, groupid):
          """ return all of the indices of particles of groupid """
          return self.indices[self.offset[groupid]
                  :self.offset[groupid]+ self.length[groupid]]