code: def intersect_ranges(self, starts, stops):
          
  
          loc = self.locate_ranges(starts, stops, strict=False)
          return self.compress(loc, axis=0)
docs: Intersect with a set of ranges.
  
          Parameters
          ----------
          starts : array_like, int
              Range start values.
          stops : array_like, int
              Range stop values.
  
          Returns
          -------
          idx : SortedIndex
  
          Examples
          --------
  
          >>> import allel
          >>> import numpy as np
          >>> idx = allel.SortedIndex([3, 6, 11, 20, 35])
          >>> ranges = np.array([[0, 2], [6, 17], [12, 15], [31, 35],
          ...                    [100, 120]])
          >>> starts = ranges[:, 0]
          >>> stops = ranges[:, 1]
          >>> idx.intersect_ranges(starts, stops)
          <SortedIndex shape=(3,) dtype=int64>
          [6, 11, 35]