Vectorizing loops with NumPy: nearest neighbors of two thousand points

Afterwards you can time code with timeit, replace Python loops by NumPy array operations and broadcasting, and split into chunks what would not fit in memory.