Tangent space embedding

Different connectivity measures, including tangent space measure
same pb np.diag(np.diag(tril_mask)) allocates a full dense matrix. Then

off_diag_mask.flat[::p + 1] -= np.diag(tril_mask)