732   const PetscInt          *is_indices;  in PCBDDCGraphSetUp()  local
857     PetscCall(ISGetIndices(ISForDofs[i], &is_indices));  in PCBDDCGraphSetUp()
864 …if (is_indices[jj] > -1 && is_indices[jj] < nvtxs) { /* out of bounds indices (if any) are skipped…  in PCBDDCGraphSetUp()
865           graph->nodes[is_indices[jj]].which_dof = k + b;  in PCBDDCGraphSetUp()
869     PetscCall(ISRestoreIndices(ISForDofs[i], &is_indices));  in PCBDDCGraphSetUp()
876     PetscCall(ISGetIndices(neumann_is, &is_indices));  in PCBDDCGraphSetUp()
878 …if (is_indices[i] > -1 && is_indices[i] < nvtxs) { /* out of bounds indices (if any) are skipped */  in PCBDDCGraphSetUp()
879         graph->nodes[is_indices[i]].special_dof = PCBDDCGRAPH_NEUMANN_MARK;  in PCBDDCGraphSetUp()
882     PetscCall(ISRestoreIndices(neumann_is, &is_indices));  in PCBDDCGraphSetUp()
888     PetscCall(ISGetIndices(dirichlet_is, &is_indices));  in PCBDDCGraphSetUp()
890 …if (is_indices[i] > -1 && is_indices[i] < nvtxs) { /* out of bounds indices (if any) are skipped */  in PCBDDCGraphSetUp()
892           graph->nodes[is_indices[i]].touched = PETSC_TRUE;  in PCBDDCGraphSetUp()
893           graph->nodes[is_indices[i]].subset  = 0;  in PCBDDCGraphSetUp()
895         graph->nodes[is_indices[i]].special_dof = PCBDDCGRAPH_DIRICHLET_MARK;  in PCBDDCGraphSetUp()
898     PetscCall(ISRestoreIndices(dirichlet_is, &is_indices));  in PCBDDCGraphSetUp()
904     PetscCall(ISGetIndices(custom_primal_vertices, &is_indices));  in PCBDDCGraphSetUp()
906 …if (is_indices[i] > -1 && is_indices[i] < nvtxs && graph->nodes[is_indices[i]].special_dof != PCBD…  in PCBDDCGraphSetUp()
907         graph->nodes[is_indices[i]].special_dof = PCBDDCGRAPH_SPECIAL_MARK - j;  in PCBDDCGraphSetUp()
911     PetscCall(ISRestoreIndices(custom_primal_vertices, &is_indices));  in PCBDDCGraphSetUp()