code: def get_ir_reciprocal_mesh(self, mesh=(10, 10, 10), is_shift=(0, 0, 0)):
          
          shift = np.array([1 if i else 0 for i in is_shift])
          mapping, grid = spglib.get_ir_reciprocal_mesh(
              np.array(mesh), self._cell, is_shift=shift, symprec=self._symprec)
  
          results = []
          for i, count in zip(*np.unique(mapping, return_counts=True)):
              results.append(((grid[i] + shift * (0.5, 0.5, 0.5)) / mesh,
                              count))
          return results
docs: k-point mesh of the Brillouin zone generated taken into account
          symmetry.The method returns the irreducible kpoints of the mesh
          and their weights
  
          Args:
              mesh (3x1 array): The number of kpoint for the mesh needed in
                  each direction
              is_shift (3x1 array): Whether to shift the kpoint grid. (1, 1,
              1) means all points are shifted by 0.5, 0.5, 0.5.
  
          Returns:
              A list of irreducible kpoints and their weights as a list of
              tuples [(ir_kpoint, weight)], with ir_kpoint given
              in fractional coordinates