docs: def is_consecutive(bins, rtol: float = 1.e-5, atol: float = 1.e-8) -> bool:
      """Check whether the bins are consecutive (edges match).
  
      Does not check if the bins are in rising order.
      """
      bins = np.asarray(bins)
      if bins.ndim == 1:
          return True
      else:
          bins = make_bin_array(bins)
          return np.allclose(bins[1:, 0], bins[:-1, 1], rtol, atol)