Remove the rows and columns corresponding to boundary condition

It is not a suggested approach, as it requires a complete restructuring of the sparsity pattern of the matrix, see for instance: https://scicomp.stackexchange.com/questions/35282/how-to-delete-nth-row-and-nth-column-of-a-matrix-k-in-petsc-and-restru