direct-solver: remove gather_block_interaction_points

We need compute the interactions with the proxy points and neighboring
points differently, so there's no point in having a function that
makes a big array out of all of them.