Hi all, Do you know what type of the resulting vector y by MatMult of a MatNest A and a Vec x, i.e. y = Ax? Because after that I want to modify a SubVec of the result y. If it is a VectNest, I think VecNestGetSubVec should be used but not VecGetSubVector. Thank you for always helps! Hui