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