only compute triangluar matrix for pairwise marginals