| /petsc/src/dm/impls/plex/transform/impls/refine/bl/ |
| H A D | plexrefbl.c | 21 PetscInt Nc, No, coff, i, ict; in DMPlexTransformSetUp_BL() local
166 coff = 8 * n; in DMPlexTransformSetUp_BL()
167 bl->cone[ct][coff + 0] = DM_POLYTOPE_SEGMENT; in DMPlexTransformSetUp_BL()
168 bl->cone[ct][coff + 1] = 1; in DMPlexTransformSetUp_BL()
169 bl->cone[ct][coff + 2] = 0; in DMPlexTransformSetUp_BL()
170 bl->cone[ct][coff + 3] = 0; in DMPlexTransformSetUp_BL()
171 bl->cone[ct][coff + 4] = DM_POLYTOPE_SEGMENT; in DMPlexTransformSetUp_BL()
172 bl->cone[ct][coff + 5] = 0; in DMPlexTransformSetUp_BL()
173 bl->cone[ct][coff + 6] = 0; in DMPlexTransformSetUp_BL()
174 bl->cone[ct][coff + 7] = DM_POLYTOPE_POINT_PRISM_TENSOR; in DMPlexTransformSetUp_BL()
[all …]
|
| /petsc/src/dm/impls/plex/transform/impls/extrude/ |
| H A D | plextrcohesive.c | 194 PetscInt rt, Nc, No, coff, ooff; in DMPlexTransformCohesiveExtrudeSetUp_Segment() local
218 coff = 8; in DMPlexTransformCohesiveExtrudeSetUp_Segment()
221 ex->cone[rt][coff + 0] = DM_POLYTOPE_SEGMENT; in DMPlexTransformCohesiveExtrudeSetUp_Segment()
222 ex->cone[rt][coff + 1] = 0; in DMPlexTransformCohesiveExtrudeSetUp_Segment()
223 ex->cone[rt][coff + 2] = 0; in DMPlexTransformCohesiveExtrudeSetUp_Segment()
224 ex->cone[rt][coff + 3] = DM_POLYTOPE_SEGMENT; in DMPlexTransformCohesiveExtrudeSetUp_Segment()
225 ex->cone[rt][coff + 4] = 0; in DMPlexTransformCohesiveExtrudeSetUp_Segment()
226 ex->cone[rt][coff + 5] = 0; in DMPlexTransformCohesiveExtrudeSetUp_Segment()
227 ex->cone[rt][coff + 6] = DM_POLYTOPE_POINT_PRISM_TENSOR; in DMPlexTransformCohesiveExtrudeSetUp_Segment()
228 ex->cone[rt][coff + 7] = 1; in DMPlexTransformCohesiveExtrudeSetUp_Segment()
[all …]
|
| H A D | plextrextrude.c | 181 PetscInt Nc, No, coff, ooff; in DMPlexTransformExtrudeSetUp_Segment() local
205 coff = 8 * Np; in DMPlexTransformExtrudeSetUp_Segment()
209 ex->cone[ct][coff + 14 * i + 0] = DM_POLYTOPE_SEGMENT; in DMPlexTransformExtrudeSetUp_Segment()
210 ex->cone[ct][coff + 14 * i + 1] = 0; in DMPlexTransformExtrudeSetUp_Segment()
211 ex->cone[ct][coff + 14 * i + 2] = i; in DMPlexTransformExtrudeSetUp_Segment()
212 ex->cone[ct][coff + 14 * i + 3] = DM_POLYTOPE_SEGMENT; in DMPlexTransformExtrudeSetUp_Segment()
213 ex->cone[ct][coff + 14 * i + 4] = 0; in DMPlexTransformExtrudeSetUp_Segment()
214 ex->cone[ct][coff + 14 * i + 5] = (i + 1) % Np; in DMPlexTransformExtrudeSetUp_Segment()
215 ex->cone[ct][coff + 14 * i + 6] = DM_POLYTOPE_POINT_PRISM_TENSOR; in DMPlexTransformExtrudeSetUp_Segment()
216 ex->cone[ct][coff + 14 * i + 7] = 1; in DMPlexTransformExtrudeSetUp_Segment()
[all …]
|
| /petsc/src/dm/impls/plex/adaptors/pragmatic/ |
| H A D | pragmaticadapt.c | 22 …PetscInt dim, cStart, cEnd, numCells, c, coff, vStart, vEnd, numVertices, numLocVertices… in DMAdaptMetric_Pragmatic_Plex() local
69 for (c = 0, coff = 0; c < numCells; ++c) { in DMAdaptMetric_Pragmatic_Plex()
75 for (cl = 0; cl < coneSize; ++cl) cells[coff++] = cone[cl] - vStart; in DMAdaptMetric_Pragmatic_Plex()
209 coff = (c - cStart) * (dim + 1); in DMAdaptMetric_Pragmatic_Plex()
213 if (bdTags[coff + d]) { in DMAdaptMetric_Pragmatic_Plex()
224 PetscCall(DMLabelSetValue(bdLabelNew, cone[perm[dim][d]], bdTags[coff + d])); in DMAdaptMetric_Pragmatic_Plex()
|
| /petsc/src/dm/impls/plex/transform/interface/ |
| H A D | plextransform.c | 1446 PetscInt c, coff = *coneoff, ooff = *orntoff; in DMPlexTransformGetCone_Internal() local
1465 …const DMPolytopeType ft = (DMPolytopeType)rcone[coff++]; /* Cell type for new cone point of pNe… in DMPlexTransformGetCone_Internal()
1466 …const PetscInt fn = rcone[coff++]; /* Number of cones of p that need to b… in DMPlexTransformGetCone_Internal()
1473 const PetscInt acp = rcone[coff++]; in DMPlexTransformGetCone_Internal()
1489 pr = rcone[coff++]; in DMPlexTransformGetCone_Internal()
1521 *coneoff = coff; in DMPlexTransformGetCone_Internal()
1541 PetscInt coff, ooff; in DMPlexTransformSetCones() local
1548 for (n = 0, coff = 0, ooff = 0; n < Nct; ++n) { in DMPlexTransformSetCones()
1553 …PetscCall(DMPlexTransformGetCone_Internal(tr, p, 0, ct, ctNew, rcone, &coff, rornt, &ooff, coneNew… in DMPlexTransformSetCones()
1575 PetscInt maxConeSize = 0, Nct, p, r, n, nr, coff = 0, ooff = 0; in DMPlexTransformGetConeOriented() local
[all …]
|
| /petsc/src/dm/impls/plex/generators/ctetgen/ |
| H A D | ctetgenerate.c | 12 PetscInt bound = numCells * numCorners, coff; in DMPlexInvertCells_CTetgen() local
21 for (coff = 0; coff < bound; coff += numCorners) SWAP(cells[coff], cells[coff + 1]); in DMPlexInvertCells_CTetgen()
|
| /petsc/src/dm/impls/plex/generators/tetgen/ |
| H A D | tetgenerate.cxx | 29 PetscInt bound = numCells * numCorners, coff; in DMPlexInvertCells_Tetgen() local
38 for (coff = 0; coff < bound; coff += numCorners) SWAP(cells[coff], cells[coff + 1]); in DMPlexInvertCells_Tetgen()
|
| /petsc/src/dm/impls/plex/adaptors/parmmg/ |
| H A D | parmmgadapt.c | 35 PetscInt dim, off, coff, maxConeSize, bdSize, i, j, k, Neq, verbosity, numIter; in DMAdaptMetric_ParMmg_Plex() local
104 for (c = 0, coff = 0; c < numCells; ++c) { in DMAdaptMetric_ParMmg_Plex()
113 cells[coff++] = vertexNumber[cone[cl] - vStart]; in DMAdaptMetric_ParMmg_Plex()
181 for (c = cStart, coff = 0; c < cEnd; ++c) { in DMAdaptMetric_ParMmg_Plex()
182 if (!cIsLeaf[c - cStart]) PetscCall(DMLabelGetValue(rgLabel, c, &cellTags[coff++])); in DMAdaptMetric_ParMmg_Plex()
|
| /petsc/src/snes/utils/dm/ |
| H A D | dminterpolatesnes.c | 893 PetscInt coff = 0, foff = 0, clSize; in DMInterpolationEvaluate() local
916 interpolant[p * ctx->dof + coff + fc] = 0.0; in DMInterpolationEvaluate()
917 …for (PetscInt f = 0; f < Nb; ++f) interpolant[p * ctx->dof + coff + fc] += xa[foff + f] * basis[(0… in DMInterpolationEvaluate()
919 coff += Nc; in DMInterpolationEvaluate()
929 … for (PetscInt fc = 0; fc < Nc; ++fc) interpolant[p * ctx->dof + coff + fc] = xa[foff + fc]; in DMInterpolationEvaluate()
930 coff += Nc; in DMInterpolationEvaluate()
935 …ck(coff == ctx->dof, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Total components %" PetscInt_FMT " != %" Pe… in DMInterpolationEvaluate()
|
| /petsc/src/dm/impls/plex/adaptors/mmg/ |
| H A D | mmgadapt.c | 30 PetscInt dim, off, coff, maxConeSize, bdSize, i, j, k, Neq, verbosity, pStart, pEnd; in DMAdaptMetric_Mmg_Plex() local
63 for (c = 0, coff = 0; c < numCells; ++c) { in DMAdaptMetric_Mmg_Plex()
69 for (cl = 0; cl < coneSize; ++cl) cells[coff++] = cone[cl] - vStart + 1; in DMAdaptMetric_Mmg_Plex()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexfem.c | 6898 const PetscInt cind = c - cS, coff = cind * totDim * totDim; in DMPlexComputeJacobianHybridByKey() local
6913 …) elemMatCoh[coff + (off + i) * totDim + j] += s[0][soff + i] * elemMatNeg[coff + (off + i) * totD… in DMPlexComputeJacobianHybridByKey()
6915 …elemMatCohP[coff + (off + i) * totDim + j] += s[0][soff + i] * elemMatNegP[coff + (off + i) * totD… in DMPlexComputeJacobianHybridByKey()
6922 …; j < totDim; ++j) elemMatCoh[coff + (off + i) * totDim + j] += elemMatNeg[coff + (off + i) * totD… in DMPlexComputeJacobianHybridByKey()
6924 …j < totDim; ++j) elemMatCohP[coff + (off + i) * totDim + j] += elemMatNegP[coff + (off + i) * totD… in DMPlexComputeJacobianHybridByKey()
6930 …for (i = 0; i < totDim * totDim; ++i) elemMatCoh[coff + i] += elemMatNeg[coff + i] + elemMatPos[co… in DMPlexComputeJacobianHybridByKey()
6932 …for (i = 0; i < totDim * totDim; ++i) elemMatCohP[coff + i] += elemMatNegP[coff + i] + elemMatPosP… in DMPlexComputeJacobianHybridByKey()
|
| H A D | plex.c | 7236 PetscInt cdof, coff, cind = 0, k; in DMPlexVecSetClosure_Depth1_Static() local
7239 PetscCall(PetscSectionGetOffset(section, cp, &coff)); in DMPlexVecSetClosure_Depth1_Static()
7240 a = &array[coff]; in DMPlexVecSetClosure_Depth1_Static()
|
| /petsc/src/ksp/pc/impls/patch/ |
| H A D | pcpatch.c | 1064 …PetscInt dof, off, cdof, coff, efdof, efoff, ifdof, ifoff, pdof, n = 0, cn = 0, ifn = 0, efn … in PCPatchCreateCellPatches() local
1069 PetscCall(PetscSectionGetOffset(cellCounts, v, &coff)); in PCPatchCreateCellPatches()
1111 if (point >= cStart && point < cEnd) cellsArray[coff + cn++] = point; in PCPatchCreateCellPatches()
1124 if (!found0 && cell0 == cellsArray[coff + n]) { in PCPatchCreateCellPatches()
1128 if (!found1 && cell1 == cellsArray[coff + n]) { in PCPatchCreateCellPatches()
|