Working on SAR MNP gibbs sampler, need to finish linear constraint matrix function