code: def get_symmetry_operations(self, cartesian=False):
          
          rotation, translation = self._get_symmetry()
          symmops = []
          mat = self._structure.lattice.matrix.T
          invmat = np.linalg.inv(mat)
          for rot, trans in zip(rotation, translation):
              if cartesian:
                  rot = np.dot(mat, np.dot(rot, invmat))
                  trans = np.dot(trans, self._structure.lattice.matrix)
              op = SymmOp.from_rotation_and_translation(rot, trans)
              symmops.append(op)
          return symmops
docs: Return symmetry operations as a list of SymmOp objects.
          By default returns fractional coord symmops.
          But cartesian can be returned too.
  
          Returns:
              ([SymmOp]): List of symmetry operations.