Not sure I understand this correctly: we're now adding 1s to the diagonal for all boundary nodes, but we're not zero'ing the rest of the row? Is this mathematically equivalent?