def occupancy(grid, points, spacing=<NUM_LIT>):
distances = ((grid[:,None,:] - points[None,:,:])**<NUM_LIT:2>).sum(axis=<NUM_LIT:2>)<EOL>occupied = (distances < spacing).sum(axis=<NUM_LIT:1>)<EOL>return occupied<EOL>
Return a vector with the occupancy of each grid point for 
    given array of points