Hi, I am thinking of using MatMPIAIJGetSeqAIJ to have a reference to the local SeqAIJ Matrix (without actually creating a new matrix). Since I don’t see any description of that function on petsc website, I want to confirm if this functions works similar to MatDenseGetLocalMatrix. Many thanks, Ian