[ENH] md.compute_neighbors

This looks like it's working. Needs code review.