Returns a set of adjacent indices. The pairs are symmetric – (x, y) and (y, x) will both be populated. Adjacency with boundary_marker doesn’t count.
