let is_subset = function
 (i1, i2) -> lower_bound (i2) <= lower_bound (i1) & upper_bound (i2) >= upper_bound (i1)