sbaij/mpi/sbaijov.c

14   IS             *is_new, *is_row;
24   PetscCall(PetscMalloc1(is_max, &is_new));
26   PetscCall(ISCompressIndicesGeneral(N, C->rmap->n, bs, is_max, is, is_new));
34     for (i = 0; i < ov; ++i) PetscCall(MatIncreaseOverlap_MPISBAIJ_Once(C, is_max, is_new));
76           PetscCall(MatCreateSubMatrices_MPIBAIJ_local(C, max_no, is_row + pos, is_new + pos, MAT_INITIAL_MATRIX, submats + pos, PETSC_TRUE));
82         PetscCall(MatIncreaseOverlap_MPIBAIJ_Once(C, is_max, is_new));
89           /* put is_new obtained from MatIncreaseOverlap_MPIBAIJ() to table */
92           PetscCall(ISGetIndices(is_new[i], &idx));
93           PetscCall(ISGetLocalSize(is_new[i], &nis));
107           PetscCall(ISRestoreIndices(is_new[i], &idx));
108           PetscCall(ISDestroy(&is_new[i]));
110           /* create updated is_new */
111           PetscCall(ISCreateGeneral(PETSC_COMM_SELF, isz, nidx, PETSC_COPY_VALUES, is_new + i));
127     PetscCall(ISGetLocalSize(is_new[i], &nis));
128     PetscCall(ISGetIndices(is_new[i], &idx));
129     PetscCall(ISCreateBlock(PetscObjectComm((PetscObject)is_new[i]), bs, nis, idx, PETSC_COPY_VALUES, &is[i]));
130     PetscCall(ISDestroy(&is_new[i]));
132   PetscCall(PetscFree(is_new));