I\'m wondering if there is a way to speed up the following piece of code using numpy. The only value I\'m interested in is the points_within_distance, the numpy array can be
