[ENH] md.compute_neighbors

This looks like it's working. Needs code review.
Am I missing a filtering step later in the code or is the haystack meant to include the query indices by default?