Thanks!
centers = np.array([ [3, 3], [4, 4] ]) point = np.array([3.5, 3.5]) vals = centers - point np.linalg.norm(vals, axis=1)