code: def get_neighbors(self, site, r, include_index=False, include_image=False):
          
          nn = self.get_sites_in_sphere(site.coords, r,
                                        include_index=include_index,
                                        include_image=include_image)
          return [d for d in nn if site != d[0]]
docs: Get all neighbors to a site within a sphere of radius r.  Excludes the
          site itself.
  
          Args:
              site (Site): Which is the center of the sphere.
              r (float): Radius of sphere.
              include_index (bool): Whether the non-supercell site index
                  is included in the returned data
              include_image (bool): Whether to include the supercell image
                  is included in the returned data
  
          Returns:
              [(site, dist) ...] since most of the time, subsequent processing
              requires the distance.
              If include_index == True, the tuple for each neighbor also includes
              the index of the neighbor.
              If include_supercell == True, the tuple for each neighbor also includes
              the index of supercell.