changing mapping from positions to cell indices to ensure consistency in floating point arrithmetic for neighboring patches
adding checks to ensure that the expected number of IB points are initialized