| /petsc/src/dm/impls/plex/tutorials/ |
| H A D | dmplexgetrestoreclosureindices.F90 | 17 PetscInt :: cStart, cEnd 36 PetscCallA(DMPlexGetHeightStratum(dm, zero, cStart, cEnd, ierr)) 43 allocate (idxMatrix(nIdx, cEnd - cStart)) 49 …PetscCallA(DMPlexGetClosureIndices(cdm, gS, gS, cEnd - 1, PETSC_TRUE, nIdx, indices, offsets, PETS… 50 idxMatrix(1:nIdx, cEnd - Cstart) = indices 51 …PetscCallA(DMPlexRestoreClosureIndices(cdm, gS, gS, cEnd - 1, PETSC_TRUE, nIdx, indices, offsets, … 56 cnt = cnt + count(idxMatrix(idx, 1) == idxMatrix(1:nIdx, cEnd))
|
| H A D | ex10.c | 23 PetscInt cStart, cEnd, c; in CreateDomainLabel() local 29 PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); in CreateDomainLabel() 30 for (c = cStart; c < cEnd; ++c) { in CreateDomainLabel() 59 PetscInt dim, Nv, v, cStart, cEnd, c; in AdaptMesh() local 109 PetscCall(DMPlexGetHeightStratum(dmCur, 0, &cStart, &cEnd)); in AdaptMesh() 110 for (c = cStart; c < cEnd; ++c) { in AdaptMesh()
|
| H A D | ex18.c | 12 PetscInt dim, cStart, cEnd, fStart, fEnd; in surfArea() local 27 PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); in surfArea() 28 for (PetscInt ii = cStart; ii < cEnd; ++ii) { in surfArea() 47 PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); in surfArea() 48 for (PetscInt ii = cStart; ii < cEnd; ++ii) { in surfArea()
|
| H A D | ex16.c | 50 PetscInt dim, cStart, cEnd; in ReadCGNSDM() local 56 PetscCall(DMPlexGetHeightStratum(*dm, 0, &cStart, &cEnd)); in ReadCGNSDM() 58 for (PetscInt i = cStart + 1; i < cEnd; i++) { in ReadCGNSDM()
|
| /petsc/src/dm/dt/dualspace/impls/refined/ |
| H A D | dualspacerefined.c | 37 PetscInt cStart, cEnd, c; in PetscDualSpaceRefinedSetCellSpaces_Refined() local 44 PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); in PetscDualSpaceRefinedSetCellSpaces_Refined() 45 for (c = 0; c < cEnd - cStart; c++) { in PetscDualSpaceRefinedSetCellSpaces_Refined() 66 PetscInt cStart, cEnd, c, spdim; in PetscDualSpaceSetUp_Refined() local 75 PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); in PetscDualSpaceSetUp_Refined() 76 for (c = cStart; c < cEnd; c++) { in PetscDualSpaceSetUp_Refined() 85 for (c = cStart; c < cEnd; c++) { in PetscDualSpaceSetUp_Refined() 94 if (c < cEnd) sp->uniform = PETSC_FALSE; in PetscDualSpaceSetUp_Refined() 134 PetscInt cStart, cEnd, c; in PetscDualSpaceRefinedView_Ascii() local 137 PetscCall(DMPlexGetHeightStratum(sp->dm, 0, &cStart, &cEnd)); in PetscDualSpaceRefinedView_Ascii() [all …]
|
| /petsc/src/dm/impls/plex/tests/ |
| H A D | ex34.c | 112 PetscInt dim, pStart, pEnd, p, cStart, cEnd, c, Nhyb = 0, off[2]; in OrderHybridMesh() local 120 PetscCall(DMPlexGetHeightStratum(*dm, 0, &cStart, &cEnd)); in OrderHybridMesh() 121 for (c = cStart; c < cEnd; ++c) { in OrderHybridMesh() 128 off[1] = cEnd - Nhyb; in OrderHybridMesh() 129 for (c = cStart; c < cEnd; ++c) { in OrderHybridMesh() 136 …ck(off[0] == cEnd - Nhyb, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Number of normal cells %" PetscInt_FMT… in OrderHybridMesh() 137 …PetscCheck(off[1] == cEnd, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Number of hybrid cells %" PetscInt_FM… in OrderHybridMesh()
|
| H A D | ex16.c | 18 PetscInt cStart, cEnd, cStartSub, cEndSub; in CreateHalfCellsLabel() local 24 PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); in CreateHalfCellsLabel() 27 cEndSub = cEnd / 2; in CreateHalfCellsLabel() 29 cStartSub = cEnd / 2; in CreateHalfCellsLabel() 30 cEndSub = cEnd; in CreateHalfCellsLabel() 41 PetscInt cStart, cEnd, cdim; in CreateHalfDomainLabel() local 48 PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); in CreateHalfDomainLabel() 50 for (PetscInt c = cStart; c < cEnd; ++c) { in CreateHalfDomainLabel()
|
| H A D | ex75.c | 8 PetscInt dim = 1, d, cStart, cEnd, c, q, degree = 2, coordSize, offset; in main() local 36 PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); in main() 37 PetscCall(PetscSectionSetChart(coordSec, cStart, cEnd)); in main() 38 for (c = cStart; c < cEnd; ++c) { in main() 50 for (c = cStart; c < cEnd; ++c) { in main()
|
| H A D | ex20.c | 10 PetscInt cStart, cEnd; in main() local 20 PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); in main() 23 if (cEnd > cStart) PetscCall(DMLabelSetValue(adaptLabel, cStart, DM_ADAPT_REFINE)); in main()
|
| H A D | ex17.c | 46 PetscInt cStart, cEnd, c; in TestCentroidLocation() local 52 PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); in TestCentroidLocation() 54 PetscCall(VecCreateSeq(PETSC_COMM_SELF, (cEnd - cStart) * cdim, &points)); in TestCentroidLocation() 57 for (c = cStart; c < cEnd; ++c) { in TestCentroidLocation() 68 if (n != (cEnd - cStart)) { in TestCentroidLocation() 72 …F, PETSC_ERR_PLIB, "Located %" PetscInt_FMT " points instead of %" PetscInt_FMT, n, cEnd - cStart); in TestCentroidLocation() 74 …for (c = cStart; c < cEnd; ++c) PetscCheck(cells[c - cStart].index == c, PETSC_COMM_SELF, PETSC_ER… in TestCentroidLocation()
|
| H A D | ex32.c | 9 PetscInt dim, depth, cStart, cEnd, c; in CheckMesh() local 14 PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); in CheckMesh() 15 for (c = cStart; c < cEnd; ++c) { in CheckMesh()
|
| H A D | ex9.c | 221 PetscInt numRuns, cStart, cEnd, c, i; in TestCone() local 235 PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); in TestCone() 238 for (c = cStart; c < cEnd; ++c) { in TestCone() 248 numRuns = (cEnd - cStart) * user->iterations; in TestCone() 264 PetscInt numRuns, cStart, cEnd, c, i; in TestTransitiveClosure() local 278 PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); in TestTransitiveClosure() 281 for (c = cStart; c < cEnd; ++c) { in TestTransitiveClosure() 293 numRuns = (cEnd - cStart) * user->iterations; in TestTransitiveClosure() 311 PetscInt numRuns, cStart, cEnd, c, i; in TestVecClosure() local 348 PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); in TestVecClosure() [all …]
|
| H A D | ex13.c | 102 PetscInt h, cStart, cEnd, c; in ScrambleOrientation() local 106 PetscCall(DMPlexGetHeightStratum(dm, h, &cStart, &cEnd)); in ScrambleOrientation() 107 for (c = cStart; c < cEnd; ++c) { in ScrambleOrientation()
|
| H A D | ex41.c | 55 PetscInt cStart, cEnd; in ConstructRefineTree() local 61 PetscCall(DMPlexGetHeightStratum(odm, 0, &cStart, &cEnd)); in ConstructRefineTree() 62 for (PetscInt c = cStart; c < cEnd; ++c) { in ConstructRefineTree()
|
| H A D | ex100.c | 68 PetscInt cStart, cEnd; in main() local 78 PetscCall(DMPlexGetHeightStratum(dm_read, 0, &cStart, &cEnd)); in main() 79 for (PetscInt cell = cStart; cell < cEnd; cell++) { in main()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexorient.c | 92 PetscInt cStart, cEnd, fStart, fEnd; in DMPlexCheckFace_Internal() local 96 PetscCall(ISGetPointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexCheckFace_Internal() 104 indC[Ns] = GetPointIndex(supp[s], cStart, cEnd, cells); in DMPlexCheckFace_Internal() 254 PetscInt cStart = 0, cEnd = 0, fStart = 0, fEnd = 0; in DMPlexOrient_Serial() local 257 if (cellIS) PetscCall(ISGetPointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexOrient_Serial() 259 PetscCall(PetscBTCreate(cEnd - cStart, &seenCells)); in DMPlexOrient_Serial() 260 PetscCall(PetscBTMemzero(cEnd - cStart, seenCells)); in DMPlexOrient_Serial() 265 for (PetscInt c = 0; c < cEnd - cStart; ++c) cellComp[c] = -1; in DMPlexOrient_Serial() 270 for (cc = cStart; cc < cEnd; ++cc) in DMPlexOrient_Serial() 272 if (cc >= cEnd) break; in DMPlexOrient_Serial() [all …]
|
| H A D | plexadapt.c | 3 …elToVolumeConstraint(DM dm, DMLabel adaptLabel, PetscInt cStart, PetscInt cEnd, PetscReal refRatio… in DMPlexLabelToVolumeConstraint() argument 10 for (c = cStart; c < cEnd; c++) { in DMPlexLabelToVolumeConstraint() 55 …elToMetricConstraint(DM dm, DMLabel adaptLabel, PetscInt cStart, PetscInt cEnd, PetscInt vStart, P… in DMPlexLabelToMetricConstraint() argument 79 for (c = cStart; c < cEnd; ++c) { in DMPlexLabelToMetricConstraint() 173 PetscInt dim, cStart, cEnd, c; in DMPlexRefine_Internal() local 182 PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); in DMPlexRefine_Internal() 222 PetscCall(PetscMalloc1(cEnd - cStart, &maxVolumes)); in DMPlexRefine_Internal() 224 … PetscCall(DMPlexLabelToVolumeConstraint(dm, adaptLabel, cStart, cEnd, PETSC_DEFAULT, maxVolumes)); in DMPlexRefine_Internal() 226 for (c = cStart; c < cEnd; ++c) { in DMPlexRefine_Internal() 233 for (c = 0; c < cEnd - cStart; ++c) maxVolumes[c] = refinementLimit; in DMPlexRefine_Internal() [all …]
|
| H A D | plexvtu.c | 56 …PetscInt dim, vStart, vEnd, cStart, cEnd, pStart, pEnd, cellHeight, numLabelCells, hasLabel, … in DMPlexGetVTKConnectivity() local 65 PetscCall(DMPlexGetHeightStratum(dm, cellHeight, &cStart, &cEnd)); in DMPlexGetVTKConnectivity() 72 for (c = cStart; c < cEnd; ++c) { in DMPlexGetVTKConnectivity() 128 PetscInt cStart, cEnd, cellHeight; in DMPlexGetNonEmptyComm_Private() local 133 PetscCall(DMPlexGetHeightStratum(dm, cellHeight, &cStart, &cEnd)); in DMPlexGetNonEmptyComm_Private() 134 color = (cStart < cEnd) ? 0 : 1; in DMPlexGetNonEmptyComm_Private() 175 …PetscInt dimEmbed, cellHeight, cStart, cEnd, vStart, vEnd, numLabelCells, hasLabel… in DMPlexVTKWriteAll_VTU() local 185 PetscCall(DMPlexGetHeightStratum(dm, cellHeight, &cStart, &cEnd)); in DMPlexVTKWriteAll_VTU() 208 for (c = cStart; c < cEnd; ++c) { in DMPlexVTKWriteAll_VTU() 282 if (cEnd > cStart) PetscCall(PetscSectionGetDof(section, cStart, &bs)); in DMPlexVTKWriteAll_VTU() [all …]
|
| /petsc/src/ts/tests/ |
| H A D | ex28.c | 65 PetscInt dim, d, cStart, cEnd, c, Np, p; in SetInitialCoordinates() local 77 PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); in SetInitialCoordinates() 84 for (c = cStart; c < cEnd; ++c) { in SetInitialCoordinates() 126 PetscInt dim, d, cStart, cEnd, c, Np, p, n; in SetInitialConditions() local 134 PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); in SetInitialConditions() 135 … == (cEnd - cStart) * Np, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "TS solution local size %" PetscInt_… in SetInitialConditions() 138 for (c = cStart; c < cEnd; ++c) { in SetInitialConditions() 153 PetscInt dim, cStart, cEnd, c, Np = user->particlesPerCell, p; in CreateParticles() local 168 PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); in CreateParticles() 169 PetscCall(DMSwarmSetLocalSizes(*sw, (cEnd - cStart) * Np, 0)); in CreateParticles() [all …]
|
| H A D | ex27.c | 68 PetscInt dim, d, cStart, cEnd, c, Np, p; in SetInitialCoordinates() local 83 PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); in SetInitialCoordinates() 91 for (c = cStart; c < cEnd; ++c) { in SetInitialCoordinates() 114 for (c = cStart; c < cEnd; ++c) { in SetInitialCoordinates() 140 PetscInt dim, d, cStart, cEnd, c, Np, p, n; in SetInitialConditions() local 148 PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); in SetInitialConditions() 151 for (c = cStart; c < cEnd; ++c) { in SetInitialConditions() 166 PetscInt dim, cStart, cEnd, c, Np = user->N, p; in CreateParticles() local 188 PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); in CreateParticles() 189 PetscCall(DMSwarmSetLocalSizes(*sw, (cEnd - cStart) * Np, 0)); in CreateParticles() [all …]
|
| /petsc/src/dm/impls/swarm/tests/ |
| H A D | ex3.c | 36 PetscInt dim, cStart, cEnd, c, Np = user->particlesPerCell, p; in CreateParticles() local 50 PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); in CreateParticles() 51 PetscCall(DMSwarmSetLocalSizes(*sw, (cEnd - cStart) * Np, 0)); in CreateParticles() 54 for (c = cStart; c < cEnd; ++c) { in CreateParticles()
|
| /petsc/src/dm/dt/fe/tests/ |
| H A D | ex1.c | 188 PetscInt depth, cStart, cEnd, cell, chunkSize = cbs, Nch = 0, Nf, f, totDim, i, k; in TestIntegration() local 198 PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); in TestIntegration() 211 for (cell = cStart; cell < cEnd; cell += chunkSize, ++Nch) { in TestIntegration() 212 const PetscInt cS = cell, cE = PetscMin(cS + chunkSize, cEnd), Ne = cE - cS; in TestIntegration() 233 PetscCall(PetscFEGeomRestoreChunk(affineGeom, cStart, cEnd, &chunkGeom)); in TestIntegration() 241 PetscInt N = (cEnd - cStart) * Nf * its; in TestIntegration() 271 PetscInt cStart, cEnd, Nf, N; in TestIntegration2() local 274 PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); in TestIntegration2() 278 N = (cEnd - cStart) * Nf * eventInfo.count; in TestIntegration2()
|
| /petsc/src/dm/field/tutorials/ |
| H A D | ex1.c | 61 static PetscErrorCode TestEvaluateFE(DMField field, PetscInt n, PetscInt cStart, PetscInt cEnd, Pet… in TestEvaluateFE() argument 78 PetscCall(PetscRandomSetInterval(rand, (PetscScalar)cStart, (PetscScalar)cEnd)); in TestEvaluateFE() 105 static PetscErrorCode TestEvaluateFV(DMField field, PetscInt n, PetscInt cStart, PetscInt cEnd, Pet… in TestEvaluateFV() argument 122 PetscCall(PetscRandomSetInterval(rand, (PetscScalar)cStart, (PetscScalar)cEnd)); in TestEvaluateFV() 228 PetscInt cStart = -1, cEnd = -1; in main() local 277 PetscCall(DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd)); in main() 326 PetscCall(DMDAGetHeightStratum(dm, 0, &cStart, &cEnd)); in main() 345 if (numFE) PetscCall(TestEvaluateFE(field, numFE, cStart, cEnd, quad, rand)); in main() 346 if (numFV) PetscCall(TestEvaluateFV(field, numFV, cStart, cEnd, rand)); in main()
|
| /petsc/src/ts/tutorials/ |
| H A D | ex52.c | 82 PetscInt cell, cStart, cEnd; in FormInitialSolution() local 89 PetscCall(DMPlexGetHeightStratum(da, 0, &cStart, &cEnd)); in FormInitialSolution() 92 for (cell = cStart; cell < cEnd; cell++) { in FormInitialSolution() 149 PetscInt cell, cStart, cEnd, nC; in FormFunction() local 176 PetscCall(DMPlexGetHeightStratum(da, 0, &cStart, &cEnd)); in FormFunction() 197 for (cell = cStart; cell < cEnd; cell++) { in FormFunction()
|
| /petsc/src/binding/petsc4py/src/petsc4py/lib/_pytypes/viewer/ |
| H A D | petscpyvista.py | 70 cStart, cEnd = plex.getHeightStratum(0) 74 for c in range(cStart, cEnd): 82 for c in range(cStart, cEnd): 94 celltypes = np.zeros((cEnd - cStart), dtype=np.uint32) 95 for c in range(cStart, cEnd):
|