code: def _get_symmetry(self):
          
          d = spglib.get_symmetry(self._cell, symprec=self._symprec,
                                  angle_tolerance=self._angle_tol)
          # Sometimes spglib returns small translation vectors, e.g.
          # [1e-4, 2e-4, 1e-4]
          # (these are in fractional coordinates, so should be small denominator
          # fractions)
          trans = []
          for t in d["translations"]:
              trans.append([float(Fraction.from_float(c).limit_denominator(1000))
                            for c in t])
          trans = np.array(trans)
  
          # fractional translations of 1 are more simply 0
          trans[np.abs(trans) == 1] = 0
          return d["rotations"], trans
docs: Get the symmetry operations associated with the structure.
  
          Returns:
              Symmetry operations as a tuple of two equal length sequences.
              (rotations, translations). "rotations" is the numpy integer array
              of the rotation matrices for scaled positions
              "translations" gives the numpy float64 array of the translation
              vectors in scaled positions.