poinc_dist_num.py

def num(points):
    expd=np.expand_dims(points,2) #need another dimension...
    tiled=np.tile(expd, points.shape[0]) #...to tile up the vectors
    trans=np.transpose(points) #Also need to transpose the points matrix to fit well with broadcasting
    diff=trans-tiled           #doing the difference, exploiting Numpy broadcasting capabilities
    num=np.sum(np.square(diff), axis=1) #an then obtain the squared norm of the difference
    return num