where hi_m gives the top-order m bits

This essentially:

• partitions each dimension into 2^m regions

• maps vector into "identifier" of region that contains it

For non-equi-width regions, use mapping array to determine partition boundaries.