make_summation_matrix function

Construct symmetry constraint matrix for bivariate symmetric smoothing.