Looks almost like COO, `coo_matrix((major_index, np.ones(len(major_index),dtype=np.intc)), value))`