impls/plex/plex.c

53   PetscInt       cStart, cEnd;
56   PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd));
57   if (cEnd <= cStart) {
61   PetscCall(DMPlexGetCellType(dm, cStart, &ct));
74 + cStart - The first "normal" cell, pass `NULL` if not needed
84 PetscErrorCode DMPlexGetSimplexOrBoxCells(DM dm, PetscInt height, PeOp PetscInt *cStart, PeOp PetscInt *cEnd)
103       PetscCall(DMPlexGetHeightStratum(dm, PetscMax(height, 0), cStart, cEnd));
120   if (cStart) *cStart = cS;
127   PetscInt                 cdim, pStart, pEnd, vStart, vEnd, cStart, cEnd, c, depth, cellHeight, t;
154     PetscCall(DMLabelGetStratumBounds(ctLabel, ict, &cStart, &cEnd));
156       PetscCall(DMLabelGetValue(depthLabel, cStart, &dep));
160       if ((cStart >= pStart) && (cStart < pEnd)) PetscCall(PetscSectionGetFieldDof(section, cStart, field, &vcdof[c]));
162       if ((cStart >= pStart) && (cStart < pEnd)) PetscCall(PetscSectionGetDof(section, cStart, &vcdof[c]));
186       PetscCall(DMLabelGetStratumBounds(ctLabel, ict, &cStart, &cEnd));
187       sStart[t] = cStart;
220   PetscInt cdim, pStart, pEnd, vStart, vEnd, cStart, cEnd;
227   PetscCall(DMPlexGetSimplexOrBoxCells(dm, 0, &cStart, &cEnd));
231     if ((cStart >= pStart) && (cStart < pEnd)) PetscCall(PetscSectionGetFieldDof(section, cStart, field, &vcdof[1]));
234     if ((cStart >= pStart) && (cStart < pEnd)) PetscCall(PetscSectionGetDof(section, cStart, &vcdof[1]));
243     *sStart = cStart;
424   PetscInt           dim, Nf, f, Nc, comp, vStart, vEnd, cStart, cEnd, c, N, level, step, w = 0;
439   PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd));
487       for (c = cStart; c < cEnd; ++c) {
914   PetscInt           dim, pStart, pEnd, cStart, cEnd, c;
922   PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd));
929   for (c = cStart; c < cEnd; ++c) {
1147     PetscInt     numLabels, l, numColors, numLColors, dim, d, depth, cStart, cEnd, c, vStart, vEnd, v, eStart = 0, eEnd = 0, fStart = 0, fEnd = 0, e, p, n;
1201     PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd));
1244       PetscInt maxStratum = PetscMax(vEnd - vStart, PetscMax(eEnd - eStart, PetscMax(fEnd - fStart, cEnd - cStart)));
1258       PetscCall(PetscViewerASCIIPrintf(viewer, "\\newcommand{\\cStart}{%" PetscInt_FMT "}\n", cStart));
1260       PetscCall(PetscViewerASCIIPrintf(viewer, "\\newcommand{\\numCells}{%" PetscInt_FMT "}\n", cEnd - cStart));
1261       PetscCall(PetscViewerASCIIPrintf(viewer, "\\newcommand{\\cShift}{%.2f}\n", 3 + (maxStratum - (cEnd - cStart)) / 2.));
1378       for (c = cStart; c < cEnd; ++c) {
1427     for (c = cStart; c < cEnd; ++c) {
1481       PetscCall(PetscViewerASCIIPrintf(viewer, "\\foreach \\c in {\\cStart,...,\\cEnd}\n"));
1483       PetscCall(PetscViewerASCIIPrintf(viewer, "  \\node(\\c_%d) [draw,shape=circle,color=%s,minimum size = 6mm] at (\\cShift+\\c-\\cStart,%d) {\\c};\n", rank, color, height++));
1538     PetscInt               i, c, cStart, cEnd, cum, numVertices, ect, ectn, cellHeight;
1562     PetscCall(DMPlexGetHeightStratum(dm, cellHeight, &cStart, &cEnd));
1566     for (c = cStart, cum = 0; c < cEnd; c++) {
1567       if (idxs[c - cStart] < 0) continue;
1568       idxs2[cum++] = idxs[c - cStart];
1950   PetscInt     dim, vStart, vEnd, cStart, cEnd, c, cDegree, edgeDiv, lineColor = PETSC_DETERMINE, cellColor = PETSC_DETERMINE;
1966   PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd));
1977   for (c = cStart; c < cEnd; ++c) {
2094       PetscInt cStart, cEnd, c;
2096       PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd));
2097       for (c = cStart; c < cEnd; ++c) PetscCall(DMSetLabelValue(dm, "Cell Sets", c, 1));
7792   PetscInt numFields, sStart, sEnd, cStart, cEnd;
7797   PetscCall(PetscSectionGetChart(cSec, &cStart, &cEnd));
7806     if (b >= cStart && b < cEnd) PetscCall(PetscSectionGetDof(cSec, b, &cSecDof));
7845   PetscInt        cStart = -1, cEnd = -1;
7950   PetscCall(PetscSectionGetChart(cSec, &cStart, &cEnd));
9030   PetscInt cellHeight, cStart, cEnd;
9034   if (includeAll) PetscCall(DMPlexGetHeightStratum(dm, cellHeight, &cStart, &cEnd));
9035   else PetscCall(DMPlexGetSimplexOrBoxCells(dm, cellHeight, &cStart, &cEnd));
9036   PetscCall(DMPlexCreateNumbering_Plex(dm, cStart, cEnd, 0, NULL, dm->sf, globalCellNumbers));
9237   PetscInt       dim, cStart, cEnd, c;
9246   PetscCall(DMPlexGetHeightStratum(rdm, 0, &cStart, &cEnd));
9247   PetscCall(DMPlexGetCellType(dm, cStart, &ct));
9257   for (c = cStart; c < cEnd; ++c) {
9292   PetscInt       dim, pStart, pEnd, p, cStart;
9303   PetscCall(DMPlexGetHeightStratum(plex, 0, &cStart, NULL));
9304   PetscCall(DMPlexGetCellType(plex, cStart, &ct));
9505   PetscInt               vStart, vEnd, cStart, cEnd, c;
9510   PetscCall(DMPlexGetHeightStratum(dm, cellHeight, &cStart, &cEnd));
9512   for (c = cStart; c < cEnd; ++c) {
9563   PetscInt               dim, depth, vStart, vEnd, cStart, cEnd, c, h;
9579     PetscCall(DMPlexGetHeightStratum(dm, h, &cStart, &cEnd));
9580     for (c = cStart; c < cEnd; ++c) {
9653   PetscInt  dim, depth, dE, d, cStart, cEnd, c;
9661   PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd));
9665   for (c = cStart; c < cEnd; ++c) {
9768     PetscInt cellHeight, cStart, cEnd;
9771     PetscCall(DMPlexGetHeightStratum(dm, cellHeight, &cStart, &cEnd));
9775       PetscCheck(point < cStart || point >= cEnd, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Point SF contains %" PetscInt_FMT " which is a cell", point);
9923   PetscInt     cdim, cStart, cEnd, c, eStart, eEnd, count = 0;
9937   PetscCall(DMPlexGetSimplexOrBoxCells(dm, 0, &cStart, &cEnd));
9939   for (c = cStart; c < cEnd; c++) {
10070   PetscInt               nc, cellHeight, cStart, cEnd, cell, cellIter = 0;
10109   PetscCall(DMPlexGetHeightStratum(dm, cellHeight, &cStart, &cEnd));
10115   PetscCall(VecSetSizes(*OrthQual, cEnd - cStart, PETSC_DETERMINE));
10125   PetscCall(PetscMalloc5(cEnd - cStart, &idx, cEnd - cStart, &oqVals, nc, &ci, nc, &fi, nc, &Ai));
10126   for (cell = cStart; cell < cEnd; cellIter++, cell++) {
10133     idx[cellIter] = cell - cStart;
10193   PetscCall(VecSetValuesLocal(*OrthQual, cEnd - cStart, idx, oqVals, INSERT_VALUES));
10995     PetscInt           cStart, cEnd, Nf, N;
10999     PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd));
11003     N        = (cEnd - cStart) * Nf * eventInfo.count;