code: def _check_indices(self, pair_inds, pair_n=2):
          
  
          pair_inds = np.array(pair_inds).astype(dtype=np.int, casting='safe')
  
          if pair_inds.ndim != 2:
              raise ValueError("pair indices has to be a matrix.")
  
          if pair_inds.shape[1] != pair_n:
              raise ValueError("pair indices shape has to be (x, %i)." % pair_n)
  
          if pair_inds.max() > self.topology.n_atoms:
              raise ValueError("index out of bounds: %i."
                               " Maximum atom index available: %i"
                               % (pair_inds.max(), self.topology.n_atoms))
  
          return pair_inds
docs: ensure pairs are valid (shapes, all atom indices available?, etc.)