now have diagonal matrices that discern whether values are on the boundary or not (for dirichlet boundary conditions)