Compute the neighbours list directly as a numpy array and start using the multiprocessing utils.