add interface for PETSc submatrix and test