x = mesh.points[edges[:, (0)]]
y = mesh.points[edges[:, (1)]]
dists = numpy.sum((x - y) ** 2, 1)
idx = numpy.argmin(dists)
return edges[idx]
dim = 3
N = 1000 * dim
pts = numpy.random.random(N).reshape(N / dim, dim)
Identifying points with the smallest Euclidean distance