code: def smear(idx, factor):
      
  
  
      s = [idx]
      for i in range(factor+1):
          a = i - factor/2
          s += [idx + a]
      return numpy.unique(numpy.concatenate(s))
docs: This function will take as input an array of indexes and return every
      unique index within the specified factor of the inputs.
  
      E.g.: smear([5,7,100],2) = [3,4,5,6,7,8,9,98,99,100,101,102]
  
      Parameters
      -----------
      idx : numpy.array of ints
          The indexes to be smeared.
      factor : idx
          The factor by which to smear out the input array.
  
      Returns
      --------
      new_idx : numpy.array of ints
          The smeared array of indexes.