code: def get_equivalent_atoms(self, tolerance=0.3):
          
          PA = self._get_point_group_analyzer(tolerance=tolerance)
          eq = PA.get_equivalent_atoms()
          self._convert_eq(eq)
          return eq
docs: Returns sets of equivalent atoms with symmetry operations
  
          Args:
              tolerance (float): Tolerance to generate the full set of symmetry
                  operations.
  
          Returns:
              dict: The returned dictionary has two possible keys:
  
              ``eq_sets``:
              A dictionary of indices mapping to sets of indices,
              each key maps to indices of all equivalent atoms.
              The keys are guaranteed to be not equivalent.
  
              ``sym_ops``:
              Twofold nested dictionary.
              ``operations[i][j]`` gives the symmetry operation
              that maps atom ``i`` unto ``j``.