28 …PetscInt        n, N, i, j, k, l, nsub, Jlow[2] = {0, 1}, *jlow, Jhigh[2] = {3, 4}, *jhigh, row, c…  in main()  local
98   for (l = 0; l < nsub; ++l) {  in main()
102       for (j = jlow[l]; j < jhigh[l]; ++j) {  in main()
107     PetscCall(ISCreateGeneral(PETSC_COMM_SELF, 12, subindices, PETSC_OWN_POINTER, rowis + l));  in main()
109       PetscCall(ISDuplicate(rowis[l], colis + l));  in main()
111       PetscCall(PetscObjectReference((PetscObject)rowis[l]));  in main()
112       colis[l] = rowis[l];  in main()
114     if (sort_rows) PetscCall(ISSort(rowis[l]));  in main()
115     if (sort_cols) PetscCall(ISSort(colis[l]));  in main()
128       for (l = 0; l < nsub; ++l) {  in main()
130 …F, "[%" PetscInt_FMT ":%" PetscInt_FMT "]: Subdomain row IS %" PetscInt_FMT ":\n", rank, size, l));  in main()
131         PetscCall(ISView(rowis[l], PETSC_VIEWER_STDOUT_SELF));  in main()
132 …F, "[%" PetscInt_FMT ":%" PetscInt_FMT "]: Subdomain col IS %" PetscInt_FMT ":\n", rank, size, l));  in main()
133         PetscCall(ISView(colis[l], PETSC_VIEWER_STDOUT_SELF));  in main()
134 …OMM_SELF, "[%" PetscInt_FMT ":%" PetscInt_FMT "]: Submatrix %" PetscInt_FMT ":\n", rank, size, l));  in main()
135         PetscCall(MatView(S[l], PETSC_VIEWER_STDOUT_SELF));  in main()
137           PetscCall(MatGetOwnershipRange(S[l], &i0, &i1));  in main()
139             PetscCall(MatGetRow(S[l], i, &ncols, &cols, NULL));  in main()
146             PetscCall(MatRestoreRow(S[l], i, &ncols, &cols, NULL));  in main()
159   for (l = 0; l < nsub; ++l) {  in main()
160     PetscCall(ISDestroy(&rowis[l]));  in main()
161     PetscCall(ISDestroy(&colis[l]));  in main()