Home
last modified time | relevance | path

Searched refs:cellIS (Results 1 – 23 of 23) sorted by relevance

/petsc/src/ts/utils/
H A Ddmplexts.c52 IS cellIS; in DMPlexTSComputeRHSFunctionFVM() local
60 PetscCall(DMGetStratumIS(plex, "dim", depth, &cellIS)); in DMPlexTSComputeRHSFunctionFVM()
61 if (!cellIS) PetscCall(DMGetStratumIS(plex, "depth", depth, &cellIS)); in DMPlexTSComputeRHSFunctionFVM()
64 PetscCall(DMPlexComputeResidualByKey(plex, key, cellIS, time, locX, NULL, time, locF, ctx)); in DMPlexTSComputeRHSFunctionFVM()
68 PetscCall(ISDestroy(&cellIS)); in DMPlexTSComputeRHSFunctionFVM()
146 IS cellIS; in DMPlexTSComputeIFunctionFEM() local
155 cellIS = allcellIS; in DMPlexTSComputeIFunctionFEM()
161 PetscCall(ISIntersect_Caching_Internal(allcellIS, pointIS, &cellIS)); in DMPlexTSComputeIFunctionFEM()
164 PetscCall(DMPlexComputeResidualByKey(plex, key, cellIS, time, locX, locX_t, time, locF, ctx)); in DMPlexTSComputeIFunctionFEM()
165 PetscCall(ISDestroy(&cellIS)); in DMPlexTSComputeIFunctionFEM()
[all …]
/petsc/src/dm/dt/fe/tests/
H A Dex2.c94 PetscErrorCode CellRangeGetFEGeom(IS cellIS, DMField coordField, PetscQuadrature quad, PetscFEGeomM… in CellRangeGetFEGeom() argument
103 PetscCall(PetscObjectQuery((PetscObject)cellIS, composeStr, (PetscObject *)&container)); in CellRangeGetFEGeom()
107 PetscCall(DMFieldCreateFEGeom(coordField, cellIS, quad, mode, geom)); in CellRangeGetFEGeom()
108 …PetscCall(PetscObjectContainerCompose((PetscObject)cellIS, composeStr, *geom, PetscFEGeomDestroy_V… in CellRangeGetFEGeom()
113 PetscErrorCode CellRangeRestoreFEGeom(IS cellIS, DMField coordField, PetscQuadrature quad, PetscBoo… in CellRangeRestoreFEGeom() argument
120 static PetscErrorCode CreateFEGeometry(DM dm, PetscDS ds, IS cellIS, PetscQuadrature *affineQuad, P… in CreateFEGeometry() argument
132 PetscCall(DMFieldGetDegree(coordField, cellIS, NULL, &maxDegree)); in CreateFEGeometry()
134 PetscCall(DMFieldCreateDefaultQuadrature(coordField, cellIS, affineQuad)); in CreateFEGeometry()
135 …if (*affineQuad) PetscCall(CellRangeGetFEGeom(cellIS, coordField, *affineQuad, PETSC_FEGEOM_BASIC,… in CreateFEGeometry()
144 … PetscCall(CellRangeGetFEGeom(cellIS, coordField, (*quads)[f], PETSC_FEGEOM_BASIC, &(*geoms)[f])); in CreateFEGeometry()
[all …]
H A Dex1.c103 PetscErrorCode CellRangeGetFEGeom(IS cellIS, DMField coordField, PetscQuadrature quad, PetscFEGeomM… in CellRangeGetFEGeom() argument
112 PetscCall(PetscObjectQuery((PetscObject)cellIS, composeStr, (PetscObject *)&container)); in CellRangeGetFEGeom()
116 PetscCall(DMFieldCreateFEGeom(coordField, cellIS, quad, mode, geom)); in CellRangeGetFEGeom()
117 …PetscCall(PetscObjectContainerCompose((PetscObject)cellIS, composeStr, *geom, PetscFEGeomDestroy_V… in CellRangeGetFEGeom()
122 PetscErrorCode CellRangeRestoreFEGeom(IS cellIS, DMField coordField, PetscQuadrature quad, PetscBoo… in CellRangeRestoreFEGeom() argument
129 static PetscErrorCode CreateFEGeometry(DM dm, PetscDS ds, IS cellIS, PetscQuadrature *affineQuad, P… in CreateFEGeometry() argument
141 PetscCall(DMFieldGetDegree(coordField, cellIS, NULL, &maxDegree)); in CreateFEGeometry()
143 PetscCall(DMFieldCreateDefaultQuadrature(coordField, cellIS, affineQuad)); in CreateFEGeometry()
144 …if (*affineQuad) PetscCall(CellRangeGetFEGeom(cellIS, coordField, *affineQuad, PETSC_FEGEOM_BASIC,… in CreateFEGeometry()
153 … PetscCall(CellRangeGetFEGeom(cellIS, coordField, (*quads)[f], PETSC_FEGEOM_BASIC, &(*geoms)[f])); in CreateFEGeometry()
[all …]
/petsc/src/dm/field/tutorials/
H A Dex1.c69 IS cellIS; in TestEvaluateFE() local
86 PetscCall(ISCreateGeneral(comm, n, cells, PETSC_OWN_POINTER, &cellIS)); in TestEvaluateFE()
87 PetscCall(PetscObjectSetName((PetscObject)cellIS, "FE Test Cells")); in TestEvaluateFE()
91 PetscCall(DMFieldEvaluateFE(field, cellIS, quad, PETSC_SCALAR, B, D, H)); in TestEvaluateFE()
92 PetscCall(DMFieldEvaluateFE(field, cellIS, quad, PETSC_REAL, rB, rD, rH)); in TestEvaluateFE()
97 PetscCall(ISView(cellIS, viewer)); in TestEvaluateFE()
101 PetscCall(ISDestroy(&cellIS)); in TestEvaluateFE()
113 IS cellIS; in TestEvaluateFV() local
130 PetscCall(ISCreateGeneral(comm, n, cells, PETSC_OWN_POINTER, &cellIS)); in TestEvaluateFV()
131 PetscCall(PetscObjectSetName((PetscObject)cellIS, "FV Test Cells")); in TestEvaluateFV()
[all …]
/petsc/src/dm/field/interface/
H A Ddmfield.c276 PetscErrorCode DMFieldEvaluateFE(DMField field, IS cellIS, PetscQuadrature points, PetscDataType da… in DMFieldEvaluateFE() argument
280 PetscValidHeaderSpecific(cellIS, IS_CLASSID, 2); in DMFieldEvaluateFE()
285 PetscUseTypeMethod(field, evaluateFE, cellIS, points, datatype, B, D, H); in DMFieldEvaluateFE()
316 PetscErrorCode DMFieldEvaluateFV(DMField field, IS cellIS, PetscDataType datatype, void *B, void *D… in DMFieldEvaluateFV() argument
320 PetscValidHeaderSpecific(cellIS, IS_CLASSID, 2); in DMFieldEvaluateFV()
324 PetscUseTypeMethod(field, evaluateFV, cellIS, datatype, B, D, H); in DMFieldEvaluateFV()
346 PetscErrorCode DMFieldGetDegree(DMField field, IS cellIS, PeOp PetscInt *minDegree, PeOp PetscInt *… in DMFieldGetDegree() argument
350 PetscValidHeaderSpecific(cellIS, IS_CLASSID, 2); in DMFieldGetDegree()
357 PetscTryTypeMethod(field, getDegree, cellIS, minDegree, maxDegree); in DMFieldGetDegree()
/petsc/src/dm/impls/plex/tests/
H A Dex26.c150 IS cellIS; in main() local
155 PetscCall(DMGetStratumIS(dm, "Cell Sets", csID[set], &cellIS)); in main()
191 PetscCall(ISGetIndices(cellIS, &cellID)); in main()
252 PetscCall(ISRestoreIndices(cellIS, &cellID)); in main()
253 PetscCall(ISDestroy(&cellIS)); in main()
418 IS cellIS; in main() local
424 PetscCall(DMGetStratumIS(dmS, "Cell Sets", csID[set], &cellIS)); in main()
425 PetscCall(ISGetIndices(cellIS, &cellID)); in main()
426 PetscCall(ISGetSize(cellIS, &numCells)); in main()
439 PetscCall(ISRestoreIndices(cellIS, &cellID)); in main()
[all …]
H A Dex26f90.F9059 type(tIS) :: cellIS
224 PetscCallA(DMGetStratumIS(dm, 'Cell Sets', csID(set), cellIS, ierr))
238 PetscCallA(ISGetIndices(cellIS, cellID, ierr))
295 PetscCallA(ISRestoreIndices(cellIS, cellID, ierr))
296 PetscCallA(ISDestroy(cellIS, ierr))
445 PetscCallA(DMGetStratumIS(dmS, 'Cell Sets', csID(set), cellIS, ierr))
446 PetscCallA(ISGetIndices(cellIS, cellID, ierr))
447 PetscCallA(ISGetSize(cellIS, numCells, ierr))
462 PetscCallA(ISRestoreIndices(cellIS, cellID, ierr))
463 PetscCallA(ISDestroy(cellIS, ierr))
H A Dex62f90.F9042 type(tIS) :: cellIS
207 PetscCallA(DMGetStratumIS(pdm, 'Cell Sets', csID(set), cellIS, ierr))
221 PetscCallA(ISGetIndices(cellIS, cellID, ierr))
278 PetscCallA(ISRestoreIndices(cellIS, cellID, ierr))
279 PetscCallA(ISDestroy(cellIS, ierr))
429 PetscCallA(DMGetStratumIS(dmS, 'Cell Sets', csID(set), cellIS, ierr))
430 PetscCallA(ISGetIndices(cellIS, cellID, ierr))
431 PetscCallA(ISGetSize(cellIS, numCells, ierr))
446 PetscCallA(ISRestoreIndices(cellIS, cellID, ierr))
447 PetscCallA(ISDestroy(cellIS, ierr))
H A Dex64.c70 IS cellIS; in main() local
75 PetscCall(DMGetStratumIS(dm, "Cell Sets", csID[set], &cellIS)); in main()
111 PetscCall(ISGetIndices(cellIS, &cellID)); in main()
172 PetscCall(ISRestoreIndices(cellIS, &cellID)); in main()
173 PetscCall(ISDestroy(&cellIS)); in main()
/petsc/src/dm/field/impls/da/
H A Ddmfieldda.c203 static PetscErrorCode DMFieldEvaluateFE_DA(DMField field, IS cellIS, PetscQuadrature points, PetscD… in DMFieldEvaluateFE_DA() argument
256 PetscCall(ISGetLocalSize(cellIS, &nCells)); in DMFieldEvaluateFE_DA()
257 PetscCall(PetscObjectTypeCompare((PetscObject)cellIS, ISSTRIDE, &isStride)); in DMFieldEvaluateFE_DA()
258 if (isStride) PetscCall(ISStrideGetInfo(cellIS, &sfirst, &stride)); in DMFieldEvaluateFE_DA()
259 else PetscCall(ISGetIndices(cellIS, &cells)); in DMFieldEvaluateFE_DA()
294 if (!isStride) PetscCall(ISRestoreIndices(cellIS, &cells)); in DMFieldEvaluateFE_DA()
302 static PetscErrorCode DMFieldEvaluateFV_DA(DMField field, IS cellIS, PetscDataType datatype, void *… in DMFieldEvaluateFV_DA() argument
332 PetscCall(ISGetLocalSize(cellIS, &numCells)); in DMFieldEvaluateFV_DA()
334 PetscCall(PetscObjectTypeCompare((PetscObject)cellIS, ISSTRIDE, &isStride)); in DMFieldEvaluateFV_DA()
335 if (isStride) PetscCall(ISStrideGetInfo(cellIS, &sfirst, &stride)); in DMFieldEvaluateFV_DA()
[all …]
/petsc/src/snes/utils/
H A Ddmplexsnes.c296 IS cellIS; in DMPlexSNESComputeResidualFEM() local
305 cellIS = allcellIS; in DMPlexSNESComputeResidualFEM()
311 PetscCall(ISIntersect_Caching_Internal(allcellIS, pointIS, &cellIS)); in DMPlexSNESComputeResidualFEM()
314 PetscCall(DMPlexComputeResidualByKey(plex, key, cellIS, PETSC_MIN_REAL, X, NULL, 0.0, F, ctx)); in DMPlexSNESComputeResidualFEM()
315 PetscCall(ISDestroy(&cellIS)); in DMPlexSNESComputeResidualFEM()
353 IS cellIS; in DMPlexSNESComputeResidualDS() local
387 cellIS = allcellIS; in DMPlexSNESComputeResidualDS()
392 PetscCall(ISIntersect_Caching_Internal(allcellIS, pointIS, &cellIS)); in DMPlexSNESComputeResidualDS()
395 …PetscCall(DMPlexComputeResidualByKey(plex, reskeys[k], cellIS, PETSC_MIN_REAL, X, NULL, 0.0, F, ct… in DMPlexSNESComputeResidualDS()
396 PetscCall(ISDestroy(&cellIS)); in DMPlexSNESComputeResidualDS()
[all …]
/petsc/src/dm/impls/plex/ftn-custom/
H A Dzplexfemf90.c21 PETSC_EXTERN void dmplexgetcellfields_(DM *dm, IS *cellIS, Vec *locX, Vec *locX_t, Vec *locA, F90Ar… in dmplexgetcellfields_() argument
27 *ierr = ISGetLocalSize(*cellIS, &numCells); in dmplexgetcellfields_()
29 *ierr = DMPlexGetCellFields(*dm, *cellIS, *locX, *locX_t, *locA, &u, &u_t, &a); in dmplexgetcellfields_()
53 PETSC_EXTERN void dmplexrestorecellfields_(DM *dm, IS *cellIS, Vec *locX, Vec *locX_t, Vec *locA, F… in dmplexrestorecellfields_() argument
63 …*ierr = DMPlexRestoreCellFields(*dm, *cellIS, *locX, NULL, NULL, &u, u_t ? &u_t : NULL, a ? &a : N… in dmplexrestorecellfields_()
/petsc/src/dm/impls/plex/
H A Dplexorient.c86 …nternal(DM dm, PetscInt *faceFIFO, PetscInt *fTop, PetscInt *fBottom, IS cellIS, IS faceIS, PetscB… in DMPlexCheckFace_Internal() argument
96 PetscCall(ISGetPointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexCheckFace_Internal()
249 static PetscErrorCode DMPlexOrient_Serial(DM dm, IS cellIS, IS faceIS, PetscInt *Ncomp, PetscInt ce… in DMPlexOrient_Serial() argument
257 if (cellIS) PetscCall(ISGetPointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexOrient_Serial()
293 …while (fTop < fBottom) PetscCall(DMPlexCheckFace_Internal(dm, faceFIFO, &fTop, &fBottom, cellIS, f… in DMPlexOrient_Serial()
330 IS cellIS, faceIS; in DMPlexOrient()
333 PetscCall(DMPlexGetAllCells_Internal(dm, &cellIS)); in DMPlexOrient()
335 PetscCall(DMPlexOrientCells_Internal(dm, cellIS, faceIS)); in DMPlexOrient()
336 PetscCall(ISDestroy(&cellIS)); in DMPlexOrient()
701 static PetscErrorCode CreateCellAndFaceIS_Private(DM dm, DMLabel label, IS *cellIS, IS *faceIS) in CreateCellAndFaceIS_Private() argument
[all …]
H A Dplexfem.c2423 IS cellIS; in DMPlexComputeIntegral_Internal() local
2443 PetscCall(ISCreateStride(PETSC_COMM_SELF, numCells, cStart, 1, &cellIS)); in DMPlexComputeIntegral_Internal()
2461 PetscCall(DMFieldGetDegree(coordField, cellIS, NULL, &maxDegree)); in DMPlexComputeIntegral_Internal()
2463 PetscCall(DMFieldCreateDefaultQuadrature(coordField, cellIS, &affineQuad)); in DMPlexComputeIntegral_Internal()
2464 …if (affineQuad) PetscCall(DMFieldCreateFEGeom(coordField, cellIS, affineQuad, PETSC_FEGEOM_BASIC, … in DMPlexComputeIntegral_Internal()
2541 …if (!affineQuad) PetscCall(DMFieldCreateFEGeom(coordField, cellIS, q, PETSC_FEGEOM_BASIC, &cgeomFE… in DMPlexComputeIntegral_Internal()
2581 PetscCall(ISDestroy(&cellIS)); in DMPlexComputeIntegral_Internal()
3751 PetscErrorCode DMPlexGetCellFields(DM dm, IS cellIS, Vec locX, PeOp Vec locX_t, PeOp Vec locA, Pets… in DMPlexGetCellFields() argument
3769 PetscCall(ISGetPointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexGetCellFields()
3814 PetscCall(ISRestorePointRange(cellIS, &cStart, &cEnd, &cells)); in DMPlexGetCellFields()
[all …]
H A Dplex.c5897 PetscErrorCode DMPlexGetAllCells_Internal(DM plex, IS *cellIS) in DMPlexGetAllCells_Internal() argument
5903 PetscCall(DMGetStratumIS(plex, "dim", depth, cellIS)); in DMPlexGetAllCells_Internal()
5904 if (!*cellIS) PetscCall(DMGetStratumIS(plex, "depth", depth, cellIS)); in DMPlexGetAllCells_Internal()
10379 IS cellIS; in DMCreateMassMatrixLumped_Plex() local
10393 PetscCall(DMGetStratumIS(dm, "depth", depth, &cellIS)); in DMCreateMassMatrixLumped_Plex()
10400 …PetscCall(DMPlexComputeJacobianActionByKey(dmc, key, cellIS, 0.0, 0.0, ones, NULL, ones, locmass, … in DMCreateMassMatrixLumped_Plex()
10401 PetscCall(ISDestroy(&cellIS)); in DMCreateMassMatrixLumped_Plex()
10427 IS cellIS; in DMCreateMassMatrix_Plex() local
10440 PetscCall(DMGetStratumIS(dmc, "depth", depth, &cellIS)); in DMCreateMassMatrix_Plex()
10446 PetscCall(DMPlexComputeJacobianByKey(dmc, key, cellIS, 0.0, 0.0, u, NULL, *mass, *mass, NULL)); in DMCreateMassMatrix_Plex()
[all …]
/petsc/src/dm/interface/
H A Ddmceed.c134 static PetscErrorCode DMCeedCreateGeometry(DM dm, IS cellIS, PetscInt *Nqdata, CeedElemRestriction … in DMCeedCreateGeometry() argument
151 PetscCall(ISGetPointRange(cellIS, &cStart, &cEnd, &cells)); in DMCeedCreateGeometry()
197 PetscErrorCode DMCeedCreate_Internal(DM dm, IS cellIS, PetscBool createGeometry, CeedQFunctionUser … in DMCeedCreate_Internal() argument
221 PetscCall(DMCeedCreateGeometry(cdm, cellIS, &Nqdata, &sd->erq, &sd->qd, &sd->geom)); in DMCeedCreate_Internal()
261 IS cellIS; in DMCeedCreate() local
265 PetscCall(DMPlexGetAllCells_Internal(plex, &cellIS)); in DMCeedCreate()
267 PetscCall(DMCeedCreate_Internal(dm, cellIS, createGeometry, func, func_source, &dm->dmceed)); in DMCeedCreate()
269 PetscCall(ISDestroy(&cellIS)); in DMCeedCreate()
/petsc/include/petsc/private/
H A Dpcpatchimpl.h75 IS cellIS; /* Temporary IS for each cell patch */ member
/petsc/src/ksp/pc/impls/patch/
H A Dpcpatch.c607 …on(PC pc, PetscErrorCode (*func)(PC pc, PetscInt point, Vec x, Vec f, IS cellIS, PetscInt n, const… in PCPatchSetComputeFunction() argument
1891 …omputeFunction_DMPlex_Private(PC pc, PetscInt patchNum, Vec x, Vec F, IS cellIS, PetscInt n, const… in PCPatchComputeFunction_DMPlex_Private() argument
1924 …PetscCall(DMPlexComputeResidual_Patch_Internal(dm, patch->patchSection, cellIS, 0.0, x, NULL, F, c… in PCPatchComputeFunction_DMPlex_Private()
1956 …l(ISCreateGeneral(PETSC_COMM_SELF, ncell, cellsArray + offset, PETSC_USE_POINTER, &patch->cellIS)); in PCPatchComputeFunction_Internal()
1957 …PetscCallBack("PCPatch callback", patch->usercomputef(pc, point, x, F, patch->cellIS, ncell * patc… in PCPatchComputeFunction_Internal()
1958 PetscCall(ISDestroy(&patch->cellIS)); in PCPatchComputeFunction_Internal()
1973 …omputeOperator_DMPlex_Private(PC pc, PetscInt patchNum, Vec x, Mat J, IS cellIS, PetscInt n, const… in PCPatchComputeOperator_DMPlex_Private() argument
2006 …uteJacobian_Patch_Internal(dm, patch->patchSection, patch->patchSection, cellIS, 0.0, 0.0, x, NULL… in PCPatchComputeOperator_DMPlex_Private()
2061 …l(ISCreateGeneral(PETSC_COMM_SELF, ncell, cellsArray + offset, PETSC_USE_POINTER, &patch->cellIS)); in PCPatchComputeOperator_Internal()
2063 …patch->usercomputeop(pc, point, x, mat, patch->cellIS, ncell * patch->totalDofsPerCell, dofsArray … in PCPatchComputeOperator_Internal()
[all …]
/petsc/src/dm/impls/plex/hdf5/
H A Dplexhdf5.c1086 …vate(DM dm, PetscInt cStart, PetscInt cEnd, IS globalCellNumbers, PetscInt *numCorners, IS *cellIS) in CreateConesIS_Private() argument
1179 …ISCreateGeneral(PetscObjectComm((PetscObject)dm), conesSize, vertices, PETSC_OWN_POINTER, cellIS)); in CreateConesIS_Private()
1180 PetscCall(PetscLayoutSetBlockSize((*cellIS)->map, *numCorners)); in CreateConesIS_Private()
1181 PetscCall(PetscObjectSetName((PetscObject)*cellIS, "cells")); in CreateConesIS_Private()
1189 IS cellIS; in DMPlexTopologyView_HDF5_XDMF_Private() local
1216 PetscCall(CreateConesIS_Private(dm, pStart, pEnd, globalCellNumbers, &numCorners, &cellIS)); in DMPlexTopologyView_HDF5_XDMF_Private()
1225 PetscCall(ISView(cellIS, viewer)); in DMPlexTopologyView_HDF5_XDMF_Private()
1226 …PetscCall(PetscViewerHDF5WriteObjectAttribute(viewer, (PetscObject)cellIS, "cell_corners", PETSC_I… in DMPlexTopologyView_HDF5_XDMF_Private()
1227 …PetscCall(PetscViewerHDF5WriteObjectAttribute(viewer, (PetscObject)cellIS, "cell_dim", PETSC_INT, … in DMPlexTopologyView_HDF5_XDMF_Private()
1228 PetscCall(ISDestroy(&cellIS)); in DMPlexTopologyView_HDF5_XDMF_Private()
/petsc/src/dm/field/impls/ds/
H A Ddmfieldds.c792 IS cellIS; in DMFieldComputeFaceData_DS() local
799 PetscCall(DMLabelGetStratumIS(depthLabel, dim + 1, &cellIS)); in DMFieldComputeFaceData_DS()
800 PetscCall(DMFieldGetDegree(field, cellIS, NULL, &maxDegree)); in DMFieldComputeFaceData_DS()
1127 PetscCall(ISDestroy(&cellIS)); in DMFieldComputeFaceData_DS()
/petsc/src/ts/tutorials/hamiltonian/
H A Dex3.c826 IS cellIS; in CreateMesh() local
835 PetscCall(ISCreateStride(PETSC_COMM_SELF, cEnd - cStart, cStart, 1, &cellIS)); in CreateMesh()
842 PetscCall(DMFieldCreateFEGeom(coordField, cellIS, quad, PETSC_FEGEOM_BASIC, &user->fegeom)); in CreateMesh()
844 PetscCall(ISDestroy(&cellIS)); in CreateMesh()
H A Dex2.c917 IS cellIS; in CreateMesh() local
926 PetscCall(ISCreateStride(PETSC_COMM_SELF, cEnd - cStart, cStart, 1, &cellIS)); in CreateMesh()
934 PetscCall(DMFieldCreateFEGeom(coordField, cellIS, quad, PETSC_FEGEOM_BASIC, &ctx->fegeom)); in CreateMesh()
936 PetscCall(ISDestroy(&cellIS)); in CreateMesh()
H A Dex4.c1425 IS cellIS; in CreateMesh() local
1442 PetscCall(ISCreateStride(PETSC_COMM_SELF, cEnd - cStart, cStart, 1, &cellIS)); in CreateMesh()
1450 PetscCall(DMFieldCreateFEGeom(coordField, cellIS, quad, PETSC_FEGEOM_BASIC, &ctx->fegeom)); in CreateMesh()
1452 PetscCall(ISDestroy(&cellIS)); in CreateMesh()