def all_col_dist(m): D = m.shape[1] d = np.zeros((D,D)) for i in xrange(D): div = m[:,i] for j in xrange(D): djv = m[:,j] d[j][i] = cosine(div,djv) return d