I have been using BVOrthogonalize in SLEPc with a symmetric positive definite matrix (say B) in my BVSetMatrix. I was wondering if I can create a shell matrix for B and use it to overload the MatMult operations. If yes, what MATOPs like MATOP_MULT, MATOP_MULT_TRANSPOSE, MATOP_GET_DIAGONAL should I provide?