[petsc-users] Shell matrix in BVOrthogonalize