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: }