Home
last modified time | relevance | path

Searched refs:closureSize (Results 1 – 25 of 32) sorted by relevance

12

/petsc/src/dm/impls/plex/adaptors/pragmatic/
H A Dpragmaticadapt.c110 PetscInt closureSize, cl; in DMAdaptMetric_Pragmatic_Plex() local
112 PetscCall(DMPlexGetTransitiveClosure(dm, bdFacesFull[f], PETSC_TRUE, &closureSize, &closure)); in DMAdaptMetric_Pragmatic_Plex()
113 for (cl = 0; cl < closureSize * 2; cl += 2) { in DMAdaptMetric_Pragmatic_Plex()
116 … PetscCall(DMPlexRestoreTransitiveClosure(dm, bdFacesFull[f], PETSC_TRUE, &closureSize, &closure)); in DMAdaptMetric_Pragmatic_Plex()
121 PetscInt closureSize, cl; in DMAdaptMetric_Pragmatic_Plex() local
123 PetscCall(DMPlexGetTransitiveClosure(dm, bdFacesFull[f], PETSC_TRUE, &closureSize, &closure)); in DMAdaptMetric_Pragmatic_Plex()
124 for (cl = 0; cl < closureSize * 2; cl += 2) { in DMAdaptMetric_Pragmatic_Plex()
127 … PetscCall(DMPlexRestoreTransitiveClosure(dm, bdFacesFull[f], PETSC_TRUE, &closureSize, &closure)); in DMAdaptMetric_Pragmatic_Plex()
/petsc/src/dm/impls/plex/
H A Dplexvtu.c85 PetscInt closureSize; in DMPlexGetVTKConnectivity() local
87 PetscCall(DMPlexGetTransitiveClosure(dm, c, PETSC_TRUE, &closureSize, &closure)); in DMPlexGetVTKConnectivity()
88 for (v = 0; v < closureSize * 2; v += 2) { in DMPlexGetVTKConnectivity()
95 PetscCall(DMPlexRestoreTransitiveClosure(dm, c, PETSC_TRUE, &closureSize, &closure)); in DMPlexGetVTKConnectivity()
219 PetscInt closureSize; in DMPlexVTKWriteAll_VTU() local
221 PetscCall(DMPlexGetTransitiveClosure(dm, c, PETSC_TRUE, &closureSize, &closure)); in DMPlexVTKWriteAll_VTU()
222 for (v = 0; v < closureSize * 2; v += 2) { in DMPlexVTKWriteAll_VTU()
228 PetscCall(DMPlexRestoreTransitiveClosure(dm, c, PETSC_TRUE, &closureSize, &closure)); in DMPlexVTKWriteAll_VTU()
466 PetscInt closureSize; in DMPlexVTKWriteAll_VTU() local
468 PetscCall(DMPlexGetTransitiveClosure(dm, c, PETSC_TRUE, &closureSize, &closure)); in DMPlexVTKWriteAll_VTU()
[all …]
H A Dplextree.c700 PetscInt closureSize, *closure = NULL; in DMPlexCreateAnchors_Tree() local
702 PetscCall(DMPlexGetTransitiveClosure(dm, ancestor, PETSC_TRUE, &closureSize, &closure)); in DMPlexCreateAnchors_Tree()
703 PetscCall(PetscSectionSetDof(aSec, p, closureSize)); in DMPlexCreateAnchors_Tree()
704 PetscCall(DMPlexRestoreTransitiveClosure(dm, ancestor, PETSC_TRUE, &closureSize, &closure)); in DMPlexCreateAnchors_Tree()
725 PetscInt j, closureSize, *closure = NULL, aOff; in DMPlexCreateAnchors_Tree() local
729 PetscCall(DMPlexGetTransitiveClosure(dm, ancestor, PETSC_TRUE, &closureSize, &closure)); in DMPlexCreateAnchors_Tree()
730 for (j = 0; j < closureSize; j++) anchors[aOff + j] = closure[2 * j]; in DMPlexCreateAnchors_Tree()
731 PetscCall(DMPlexRestoreTransitiveClosure(dm, ancestor, PETSC_TRUE, &closureSize, &closure)); in DMPlexCreateAnchors_Tree()
1238 PetscInt closureSize, closureSizeP, *closure = NULL, *closureP = NULL; in DMPlexComputeAnchorMatrix_Tree_Direct() local
1243 PetscCall(DMPlexGetTransitiveClosure(dm, c, PETSC_TRUE, &closureSize, &closure)); in DMPlexComputeAnchorMatrix_Tree_Direct()
[all …]
H A Dplexsubmesh.c228 PetscInt closureSize, c; in DMPlexLabelComplete_Internal() local
233 PetscCall(DMPlexGetTransitiveClosure(dm, q, PETSC_TRUE, &closureSize, &closure)); in DMPlexLabelComplete_Internal()
234 … for (c = 0; c < closureSize * 2; c += 2) PetscCall(DMLabelSetValue(label, closure[c], values[v])); in DMPlexLabelComplete_Internal()
235 PetscCall(DMPlexRestoreTransitiveClosure(dm, q, PETSC_TRUE, &closureSize, &closure)); in DMPlexLabelComplete_Internal()
345 PetscInt closureSize, cl, h, cStart, cEnd; in DMPlexLabelAddCells() local
359 PetscCall(DMPlexGetTransitiveClosure(dm, points[p], PETSC_FALSE, &closureSize, &closure)); in DMPlexLabelAddCells()
360 for (cl = closureSize - 1; cl > 0; --cl) { in DMPlexLabelAddCells()
367 PetscCall(DMPlexRestoreTransitiveClosure(dm, points[p], PETSC_FALSE, &closureSize, &closure)); in DMPlexLabelAddCells()
417 PetscInt closureSize, cl; in DMPlexLabelAddFaceCells() local
420 PetscCall(DMPlexGetTransitiveClosure(dm, face, PETSC_FALSE, &closureSize, &closure)); in DMPlexLabelAddFaceCells()
[all …]
H A Dplexadapt.c12 PetscInt closureSize = 0, cl; in DMPlexLabelToVolumeConstraint() local
20 PetscCall(DMPlexGetTransitiveClosure(dm, c, PETSC_TRUE, &closureSize, &closure)); in DMPlexLabelToVolumeConstraint()
21 for (cl = 0; cl < closureSize * 2; cl += 2) { in DMPlexLabelToVolumeConstraint()
43 PetscCall(DMPlexRestoreTransitiveClosure(dm, c, PETSC_TRUE, &closureSize, &closure)); in DMPlexLabelToVolumeConstraint()
H A Dplex.c931 PetscInt closureSize, cl, ct; in DMPlexView_Ascii_Geometry() local
935 PetscCall(DMPlexGetTransitiveClosure(dm, c, PETSC_TRUE, &closureSize, &closure)); in DMPlexView_Ascii_Geometry()
937 for (cl = 0; cl < closureSize * 2; cl += 2) { in DMPlexView_Ascii_Geometry()
956 PetscCall(DMPlexRestoreTransitiveClosure(dm, c, PETSC_TRUE, &closureSize, &closure)); in DMPlexView_Ascii_Geometry()
1214 PetscInt closureSize; in DMPlexView_Ascii() local
1219 PetscCall(DMPlexGetTransitiveClosure(dm, c, PETSC_TRUE, &closureSize, &closure)); in DMPlexView_Ascii()
1220 for (p = 0; p < closureSize * 2; p += 2) PetscCall(PetscBTSet(wp, closure[p] - pStart)); in DMPlexView_Ascii()
1221 PetscCall(DMPlexRestoreTransitiveClosure(dm, c, PETSC_TRUE, &closureSize, &closure)); in DMPlexView_Ascii()
1395 PetscInt closureSize, Nv = 0, v; in DMPlexView_Ascii() local
1397 PetscCall(DMPlexGetTransitiveClosure(dm, c, PETSC_TRUE, &closureSize, &closure)); in DMPlexView_Ascii()
[all …]
H A Dplexinterpolate.c1819 PetscInt *closure = NULL, closureSize, cl, coneSize = 0; in DMPlexUninterpolate() local
1821 PetscCall(DMPlexGetTransitiveClosure(dm, c, PETSC_TRUE, &closureSize, &closure)); in DMPlexUninterpolate()
1822 for (cl = 0; cl < closureSize * 2; cl += 2) { in DMPlexUninterpolate()
1827 PetscCall(DMPlexRestoreTransitiveClosure(dm, c, PETSC_TRUE, &closureSize, &closure)); in DMPlexUninterpolate()
1834 PetscInt *closure = NULL, closureSize, cl, coneSize = 0; in DMPlexUninterpolate() local
1836 PetscCall(DMPlexGetTransitiveClosure(dm, c, PETSC_TRUE, &closureSize, &closure)); in DMPlexUninterpolate()
1837 for (cl = 0; cl < closureSize * 2; cl += 2) { in DMPlexUninterpolate()
1842 PetscCall(DMPlexRestoreTransitiveClosure(dm, c, PETSC_TRUE, &closureSize, &closure)); in DMPlexUninterpolate()
/petsc/src/dm/impls/plex/adaptors/mmg/
H A Dmmgadapt.c94 PetscInt *closure = NULL, closureSize, cl; in DMAdaptMetric_Mmg_Plex() local
100 PetscCall(DMPlexGetTransitiveClosure(dm, f, PETSC_TRUE, &closureSize, &closure)); in DMAdaptMetric_Mmg_Plex()
101 for (cl = 0; cl < closureSize * 2; cl += 2) { in DMAdaptMetric_Mmg_Plex()
104 PetscCall(DMPlexRestoreTransitiveClosure(dm, f, PETSC_TRUE, &closureSize, &closure)); in DMAdaptMetric_Mmg_Plex()
109 PetscInt *closure = NULL, closureSize, cl; in DMAdaptMetric_Mmg_Plex() local
114 PetscCall(DMPlexGetTransitiveClosure(dm, f, PETSC_TRUE, &closureSize, &closure)); in DMAdaptMetric_Mmg_Plex()
115 for (cl = 0; cl < closureSize * 2; cl += 2) { in DMAdaptMetric_Mmg_Plex()
118 PetscCall(DMPlexRestoreTransitiveClosure(dm, f, PETSC_TRUE, &closureSize, &closure)); in DMAdaptMetric_Mmg_Plex()
/petsc/src/dm/impls/plex/tests/
H A Dex26.c152 PetscInt numCells, cell, closureSize, *closureA = NULL; in main() local
192 PetscCall(DMPlexGetTransitiveClosure(dm, cellID[0], PETSC_TRUE, &closureSize, &closureA)); in main()
193 switch (closureSize) { in main()
231 …_COMM_SELF, PETSC_ERR_ARG_INCOMP, "Unknown element with closure size %" PetscInt_FMT, closureSize); in main()
233 PetscCall(DMPlexRestoreTransitiveClosure(dm, cellID[0], PETSC_TRUE, &closureSize, &closureA)); in main()
238 PetscCall(DMPlexGetTransitiveClosure(dm, cellID[cell], PETSC_TRUE, &closureSize, &closure)); in main()
239 for (p = 0; p < closureSize; ++p) { in main()
250 … PetscCall(DMPlexRestoreTransitiveClosure(dm, cellID[cell], PETSC_TRUE, &closureSize, &closure)); in main()
H A Dex64.c72 PetscInt numCells, cell, closureSize, *closureA = NULL; in main() local
112 PetscCall(DMPlexGetTransitiveClosure(dm, cellID[0], PETSC_TRUE, &closureSize, &closureA)); in main()
113 switch (closureSize) { in main()
151 …_COMM_SELF, PETSC_ERR_ARG_INCOMP, "Unknown element with closure size %" PetscInt_FMT, closureSize); in main()
153 PetscCall(DMPlexRestoreTransitiveClosure(dm, cellID[0], PETSC_TRUE, &closureSize, &closureA)); in main()
158 PetscCall(DMPlexGetTransitiveClosure(dm, cellID[cell], PETSC_TRUE, &closureSize, &closure)); in main()
159 for (p = 0; p < closureSize; ++p) { in main()
170 … PetscCall(DMPlexRestoreTransitiveClosure(dm, cellID[cell], PETSC_TRUE, &closureSize, &closure)); in main()
H A Dex9.c283 PetscInt closureSize; in TestTransitiveClosure() local
285 PetscCall(DMPlexGetTransitiveClosure(dm, c, PETSC_TRUE, &closureSize, &closure)); in TestTransitiveClosure()
286 PetscCall(DMPlexRestoreTransitiveClosure(dm, c, PETSC_TRUE, &closureSize, &closure)); in TestTransitiveClosure()
354 PetscInt closureSize = 64; in TestVecClosure() local
356 PetscCall(DMPlexVecGetClosure(dm, s, v, c, &closureSize, &closure)); in TestVecClosure()
357 … if (!user->reuseArray) PetscCall(DMPlexVecRestoreClosure(dm, s, v, c, &closureSize, &closure)); in TestVecClosure()
H A Dex60.c91 PetscInt *closure = NULL, closureSize, cl; in main() local
94 PetscCall(DMPlexGetTransitiveClosure(dm, f, PETSC_TRUE, &closureSize, &closure)); in main()
95 for (cl = 0; cl < closureSize * 2; cl += 2) { in main()
104 PetscCall(DMPlexRestoreTransitiveClosure(dm, f, PETSC_TRUE, &closureSize, &closure)); in main()
H A Dex26f90.F9061 PetscInt :: numCells, cell, closureSize
361 closureSize = size(xyz)
363 do j = 0, closureSize - 1, sdim
366 cval(offUA + i) = cval(offUA + i)*sdim/closureSize
H A Dex62f90.F9044 PetscInt :: numCells, cell, closureSize
345 closureSize = size(xyz)
347 do j = 0, closureSize - 1, sdim
350 cval(offUA + i) = cval(offUA + i)*sdim/closureSize
/petsc/src/dm/dt/tests/
H A Dex4.c10 …PetscInt nFunc, *ids, *idsCopy, *idsCopy2, i, closureSize, *closure = NULL, offset, de… in CheckSymmetry() local
50 PetscCall(DMPlexGetTransitiveClosure(dm, 0, PETSC_TRUE, &closureSize, &closure)); in CheckSymmetry()
52 for (i = 0, offset = 0; i < closureSize; i++, offset += numDofs[depth]) { in CheckSymmetry()
108 PetscCall(DMPlexRestoreTransitiveClosure(dm, 0, PETSC_TRUE, &closureSize, &closure)); in CheckSymmetry()
/petsc/src/dm/impls/plex/generators/triangle/
H A Dtrigenerate.c297 PetscInt closureSize; in DMPlexRefine_Triangle() local
299 PetscCall(DMPlexGetTransitiveClosure(dm, c, PETSC_TRUE, &closureSize, &closure)); in DMPlexRefine_Triangle()
300 …closureSize != 4) || !(closureSize != 7), comm, PETSC_ERR_ARG_WRONG, "Mesh has cell which is not a… in DMPlexRefine_Triangle()
301 …for (v = 0; v < 3; ++v) PetscCall(PetscCIntCast(closure[(v + closureSize - 3) * 2] - vStart, &in.t… in DMPlexRefine_Triangle()
302 PetscCall(DMPlexRestoreTransitiveClosure(dm, c, PETSC_TRUE, &closureSize, &closure)); in DMPlexRefine_Triangle()
/petsc/src/dm/impls/plex/exodusii/
H A Dplexexodusii2.c874 PetscInt csSize, closureSize; in DMView_PlexExodusII() local
879 PetscCall(DMPlexVecGetClosure(cdm, NULL, coord, cells[0], &closureSize, &xyz)); in DMView_PlexExodusII()
882 …closureSize == 2 * dim, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Number of vertices %" PetscInt… in DMView_PlexExodusII()
886 if (closureSize == 3 * dim) { in DMView_PlexExodusII()
888 } else if (closureSize == 4 * dim) { in DMView_PlexExodusII()
890 …es %" PetscInt_FMT " in dimension %" PetscInt_FMT " has no ExodusII type", closureSize / dim, dim); in DMView_PlexExodusII()
893 if (closureSize == 4 * dim) { in DMView_PlexExodusII()
895 } else if (closureSize == 8 * dim) { in DMView_PlexExodusII()
897 …es %" PetscInt_FMT " in dimension %" PetscInt_FMT " has no ExodusII type", closureSize / dim, dim); in DMView_PlexExodusII()
908 PetscCall(DMPlexVecRestoreClosure(cdm, NULL, coord, cells[0], &closureSize, &xyz)); in DMView_PlexExodusII()
[all …]
/petsc/src/dm/impls/plex/adaptors/parmmg/
H A Dparmmgadapt.c142 PetscInt *closure = NULL, closureSize, cl; in DMAdaptMetric_ParMmg_Plex() local
157 PetscCall(DMPlexGetTransitiveClosure(dm, f, PETSC_TRUE, &closureSize, &closure)); in DMAdaptMetric_ParMmg_Plex()
158 for (cl = 0; cl < closureSize * 2; cl += 2) { in DMAdaptMetric_ParMmg_Plex()
161 PetscCall(DMPlexRestoreTransitiveClosure(dm, f, PETSC_TRUE, &closureSize, &closure)); in DMAdaptMetric_ParMmg_Plex()
165 PetscInt *closure = NULL, closureSize, cl; in DMAdaptMetric_ParMmg_Plex() local
169 PetscCall(DMPlexGetTransitiveClosure(dm, f, PETSC_TRUE, &closureSize, &closure)); in DMAdaptMetric_ParMmg_Plex()
170 for (cl = 0; cl < closureSize * 2; cl += 2) { in DMAdaptMetric_ParMmg_Plex()
173 PetscCall(DMPlexRestoreTransitiveClosure(dm, f, PETSC_TRUE, &closureSize, &closure)); in DMAdaptMetric_ParMmg_Plex()
/petsc/src/dm/impls/plex/tutorials/
H A Dex6.c127 PetscInt closureSize, ki, kj, f, c, foff = 0; in ReadData2D() local
129 PetscCall(DMPlexVecGetClosure(dm, NULL, u, cell, &closureSize, &closure)); in ReadData2D()
146 PetscCall(DMPlexVecRestoreClosure(dm, NULL, u, cell, &closureSize, &closure)); in ReadData2D()
160 PetscInt closureSize, ki, kj, kk, f, c, foff = 0; in ReadData3D() local
162 PetscCall(DMPlexVecGetClosure(dm, NULL, u, cell, &closureSize, &closure)); in ReadData3D()
182 PetscCall(DMPlexVecRestoreClosure(dm, NULL, u, cell, &closureSize, &closure)); in ReadData3D()
/petsc/src/dm/dt/fe/impls/composite/
H A Dfecomposite.c46 PetscInt closureSize; in PetscFESetUp_Composite() local
49 PetscCall(DMPlexGetTransitiveClosure(K, s, PETSC_TRUE, &closureSize, &closure)); in PetscFESetUp_Composite()
50 for (j = 0; j < closureSize; j++) { in PetscFESetUp_Composite()
58 PetscCall(DMPlexRestoreTransitiveClosure(K, s, PETSC_TRUE, &closureSize, &closure)); in PetscFESetUp_Composite()
/petsc/src/dm/impls/plex/generators/ctetgen/
H A Dctetgenerate.c391 PetscInt closureSize; in DMPlexRefine_CTetgen() local
393 PetscCall(DMPlexGetTransitiveClosure(dm, c, PETSC_TRUE, &closureSize, &closure)); in DMPlexRefine_CTetgen()
394 …closureSize == 5) || (closureSize == 15), comm, PETSC_ERR_ARG_WRONG, "Mesh has cell which is not a… in DMPlexRefine_CTetgen()
395 …++v) in->tetrahedronlist[idx * in->numberofcorners + v] = closure[(v + closureSize - 4) * 2] - vSt… in DMPlexRefine_CTetgen()
396 PetscCall(DMPlexRestoreTransitiveClosure(dm, c, PETSC_TRUE, &closureSize, &closure)); in DMPlexRefine_CTetgen()
/petsc/src/dm/impls/plex/generators/tetgen/
H A Dtetgenerate.cxx420 PetscInt closureSize; in DMPlexRefine_Tetgen() local
422 PetscCall(DMPlexGetTransitiveClosure(dm, c, PETSC_TRUE, &closureSize, &closure)); in DMPlexRefine_Tetgen()
423 …closureSize != 5) || !(closureSize != 15), comm, PETSC_ERR_ARG_WRONG, "Mesh has cell which is not … in DMPlexRefine_Tetgen()
424 …; ++v) in.tetrahedronlist[idx * in.numberofcorners + v] = closure[(v + closureSize - 4) * 2] - vSt… in DMPlexRefine_Tetgen()
425 PetscCall(DMPlexRestoreTransitiveClosure(dm, c, PETSC_TRUE, &closureSize, &closure)); in DMPlexRefine_Tetgen()
/petsc/src/dm/impls/forest/tests/
H A Dex2.c97 PetscInt cl, closureSize, *closure = NULL; in IdentifyBadPoints() local
114 PetscCall(DMPlexGetTransitiveClosure(dmplex, p, PETSC_TRUE, &closureSize, &closure)); in IdentifyBadPoints()
115 for (cl = 0; cl < closureSize; cl++) { in IdentifyBadPoints()
122 PetscCall(DMPlexRestoreTransitiveClosure(dmplex, p, PETSC_TRUE, &closureSize, &closure)); in IdentifyBadPoints()
/petsc/src/dm/impls/forest/p4est/
H A Dpforest.h1516 PetscInt *star = NULL, *closure = NULL, closureSize, starSize, cttSize; in DMPlexCreateConnectivity_pforest() local
1538 PetscCall(DMPlexGetTransitiveClosure(dm, p, PETSC_TRUE, &closureSize, &closure)); in DMPlexCreateConnectivity_pforest()
1539 …closureSize == P4EST_INSUL, PETSC_COMM_SELF, PETSC_ERR_ARG_INCOMP, "Cell %" PetscInt_FMT " with wr… in DMPlexCreateConnectivity_pforest()
1546 PetscCall(DMPlexRestoreTransitiveClosure(dm, p, PETSC_TRUE, &closureSize, &closure)); in DMPlexCreateConnectivity_pforest()
1569 PetscCall(DMPlexGetTransitiveClosure(dm, p, PETSC_TRUE, &closureSize, &closure)); in DMPlexCreateConnectivity_pforest()
1570 …PetscCheck(closureSize == P4EST_INSUL, PETSC_COMM_SELF, PETSC_ERR_ARG_INCOMP, "Cell with wrong clo… in DMPlexCreateConnectivity_pforest()
1577 PetscCall(DMPlexRestoreTransitiveClosure(dm, p, PETSC_TRUE, &closureSize, &closure)); in DMPlexCreateConnectivity_pforest()
1698 PetscCall(DMPlexGetTransitiveClosure(dm, p, PETSC_TRUE, &closureSize, &closure)); in DMPlexCreateConnectivity_pforest()
1699 … PetscCheck(closureSize == P4EST_INSUL, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Non-standard closure"); in DMPlexCreateConnectivity_pforest()
1722 PetscCall(DMPlexRestoreTransitiveClosure(dm, p, PETSC_TRUE, &closureSize, &closure)); in DMPlexCreateConnectivity_pforest()
[all …]
/petsc/src/dm/dt/dualspace/impls/lagrange/
H A Ddspacelagrange.c384 PetscInt closureSize = 0; in PetscLagNodeIndicesComputeVertexOrder() local
432 PetscCall(DMPlexGetTransitiveClosure(dm, 0, PETSC_TRUE, &closureSize, &closure)); in PetscLagNodeIndicesComputeVertexOrder()
433 c = closureSize - nVerts; in PetscLagNodeIndicesComputeVertexOrder()
436 PetscCall(DMPlexRestoreTransitiveClosure(dm, 0, PETSC_TRUE, &closureSize, &closure)); in PetscLagNodeIndicesComputeVertexOrder()
545 PetscInt closureSize = 0; in PetscLagNodeIndicesPushForward() local
571 PetscCall(DMPlexGetTransitiveClosure_Internal(dm, p, ornt, PETSC_TRUE, &closureSize, &closure)); in PetscLagNodeIndicesPushForward()
572 c = closureSize - nSubVert; in PetscLagNodeIndicesPushForward()
575 PetscCall(DMPlexRestoreTransitiveClosure(dm, p, PETSC_TRUE, &closureSize, &closure)); in PetscLagNodeIndicesPushForward()
2629 PetscInt closureSize = 0; in PetscDualSpaceGetBoundarySymmetries_Internal() local
2634 PetscCall(DMPlexGetTransitiveClosure(sp->dm, 0, PETSC_TRUE, &closureSize, &closure)); in PetscDualSpaceGetBoundarySymmetries_Internal()
[all …]

12