2: #include <../src/mat/impls/aij/mpi/mpiaij.h> 4: PetscErrorCode MatGetDiagonalBlock_MPIAIJ(Mat A,Mat *a) 5: { 6: *a = ((Mat_MPIAIJ*)A->data)->A; 7: return 0; 8: }