MatDenseGetColumnVec() locks the whole matrix
We can lock individuals columns with a PetscBT, using the global lock as a count of the existing locks.
I need this to implement a dual basis algorithm
We can lock individuals columns with a PetscBT, using the global lock as a count of the existing locks.
I need this to implement a dual basis algorithm