Use new submatrix
rewrite of submatrix to use csc, and delete subsample