code: def _unique_coords(self, coords_in, magmoms_in=None, lattice=None):
          
          coords = []
          if magmoms_in:
              magmoms = []
              if len(magmoms_in) != len(coords_in):
                  raise ValueError
              for tmp_coord, tmp_magmom in zip(coords_in, magmoms_in):
                  for op in self.symmetry_operations:
                      coord = op.operate(tmp_coord)
                      coord = np.array([i - math.floor(i) for i in coord])
                      if isinstance(op, MagSymmOp):
                          # Up to this point, magmoms have been defined relative
                          # to crystal axis. Now convert to Cartesian and into
                          # a Magmom object.
                          magmom = Magmom.from_moment_relative_to_crystal_axes(
                              op.operate_magmom(tmp_magmom),
                              lattice=lattice
                          )
                      else:
                          magmom = Magmom(tmp_magmom)
                      if not in_coord_list_pbc(coords, coord,
                                               atol=self._site_tolerance):
                          coords.append(coord)
                          magmoms.append(magmom)
              return coords, magmoms
          else:
              for tmp_coord in coords_in:
                  for op in self.symmetry_operations:
                      coord = op.operate(tmp_coord)
                      coord = np.array([i - math.floor(i) for i in coord])
                      if not in_coord_list_pbc(coords, coord,
                                               atol=self._site_tolerance):
                          coords.append(coord)
              return coords, [Magmom(0)] * len(coords)
docs: Generate unique coordinates using coord and symmetry positions
          and also their corresponding magnetic moments, if supplied.