Home
last modified time | relevance | path

Searched refs:numCells (Results 1 – 25 of 43) sorted by relevance

12

/petsc/src/dm/impls/plex/hdf5/
H A Dplexhdf5xdmf.c67 PetscInt spatialDim, topoDim = -1, numCells, numVertices, NVertices, numCorners; in DMPlexLoad_HDF5_Xdmf_Internal() local
94 PetscCall(PetscViewerHDF5ReadSizes(viewer, topo_name, NULL, &numCells)); in DMPlexLoad_HDF5_Xdmf_Internal()
95 PetscCall(PetscLayoutSetSize(cells->map, numCells)); in DMPlexLoad_HDF5_Xdmf_Internal()
96 numCells = rank == 0 ? numCells : 0; in DMPlexLoad_HDF5_Xdmf_Internal()
97 PetscCall(PetscLayoutSetLocalSize(cells->map, numCells)); in DMPlexLoad_HDF5_Xdmf_Internal()
100 PetscCall(ISGetLocalSize(cells, &numCells)); in DMPlexLoad_HDF5_Xdmf_Internal()
104 numCells /= numCorners; in DMPlexLoad_HDF5_Xdmf_Internal()
124 …nt_FMT " numVertices %" PetscInt_FMT " spatialDim %" PetscInt_FMT "\n", numCells, numCorners, numV… in DMPlexLoad_HDF5_Xdmf_Internal()
146 …PetscCall(DMPlexBuildFromCellListParallel(dm, numCells, numVertices, NVertices, numCorners, cells_… in DMPlexLoad_HDF5_Xdmf_Internal()
/petsc/src/dm/impls/plex/tests/
H A Dex18.c320 PetscInt numCells, numVertices, network; in CreateMesh_1D() local
329 numCells = 3; in CreateMesh_1D()
330 PetscCall(PetscOptionsGetInt(NULL, NULL, "-ncells", &numCells, NULL)); in CreateMesh_1D()
331 …PetscCheck(numCells >= 3, comm, PETSC_ERR_ARG_OUTOFRANGE, "Test ncells %" PetscInt_FMT " must >=3"… in CreateMesh_1D()
334 numVertices = numCells + 1; in CreateMesh_1D()
335 PetscCall(PetscMalloc2(2 * numCells, &cells, 2 * numVertices, &coords)); in CreateMesh_1D()
336 for (i = 0; i < numCells; i++) { in CreateMesh_1D()
343 …PetscCall(DMPlexCreateFromCellListPetsc(comm, user->dim, numCells, numVertices, numCorners, PETSC_… in CreateMesh_1D()
353 numCells = 2; in CreateMesh_1D()
354 numVertices = numCells; in CreateMesh_1D()
[all …]
H A Dex1f90.F9011 PetscInt c, firstCell, numCells
23 numCells = 2
25 numPoints = numCells + numVertices
27 do c = firstCell, numCells - 1
H A Dex8.c49 PetscInt dim, cStart, cEnd, numCells, n; in ProcessOptions() local
55 numCells = cEnd - cStart; in ProcessOptions()
56 …tscCall(PetscMalloc4(numCells * dim, &options->v0, numCells * dim * dim, &options->J, numCells * d… in ProcessOptions()
57 PetscCall(PetscMalloc1(numCells * dim, &options->centroid)); in ProcessOptions()
58 PetscCall(PetscMalloc1(numCells * dim, &options->normal)); in ProcessOptions()
59 PetscCall(PetscMalloc1(numCells, &options->vol)); in ProcessOptions()
60 n = numCells * dim; in ProcessOptions()
62 …Flg || n == numCells * dim, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Invalid size of v0 %" PetscInt_FM… in ProcessOptions()
63 n = numCells * dim * dim; in ProcessOptions()
65 …g || n == numCells * dim * dim, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Invalid size of J %" PetscInt… in ProcessOptions()
[all …]
H A Ddmplexcomputecellgeometryfem.F908 PetscInt :: numCells = 1
39 …PetscCallA(DMPlexCreateFromCellListPetsc(PETSC_COMM_WORLD, numDim, numCells, numVertices, numCorne…
69 numCells = 1
73 …PetscCallA(DMPlexCreateFromCellListPetsc(PETSC_COMM_WORLD, numDim, numCells, numVertices, numCorne…
H A Dex26.c152 PetscInt numCells, cell, closureSize, *closureA = NULL; in main() local
154 PetscCall(DMGetStratumSize(dm, "Cell Sets", csID[set], &numCells)); in main()
156 if (numCells > 0) { in main()
235 for (cell = 0; cell < numCells; cell++) { in main()
420 PetscInt numCells, cell; in main() local
426 PetscCall(ISGetSize(cellIS, &numCells)); in main()
427 for (cell = 0; cell < numCells; cell++) { in main()
H A Dex26f90.F9061 PetscInt :: numCells, cell, closureSize
223 PetscCallA(DMGetStratumSize(dm, 'Cell Sets', csID(set), numCells, ierr))
225 if (numCells > 0) then
279 do cell = 1, numCells!
447 PetscCallA(ISGetSize(cellIS, numCells, ierr))
448 do cell = 1, numCells
H A Dex62f90.F9044 PetscInt :: numCells, cell, closureSize
206 PetscCallA(DMGetStratumSize(pdm, 'Cell Sets', csID(set), numCells, ierr))
208 if (numCells > 0) then
262 do cell = 1, numCells!
431 PetscCallA(ISGetSize(cellIS, numCells, ierr))
432 do cell = 1, numCells
/petsc/lib/petsc/bin/
H A Dpetsc_gen_xdmf.py23 def writeCells(self, fp, topologyPath, numCells, numCorners, cellsName = "cells"): argument
32 ''' % (cellsName, numCells, numCorners, topologyPath))
74 def writeSpaceGridHeader(self, fp, numCells, numCorners, cellDim, spaceDim, cellsName = "cells"): argument
89 ''' % (self.cellMap[cellDim][numCorners], numCells, cellsName, "XYZ" if spaceDim > 2 else "XY"))
205 …def write(self, hdfFilename, topologyPath, numCells, numCorners, cellDim, htopologyPath, numHCells… argument
212 self.writeCells(fp, topologyPath, numCells, numCorners)
222 self.writeSpaceGridHeader(fp, numCells, numCorners, cellDim, spaceDim)
268 numCells = cells.shape[0]
294 …Xdmf(xdmfFilename).write(hdfFilename, topoPath, numCells, numCorners, cellDim, htopoPath, numHCell…
/petsc/src/dm/impls/plex/generators/ctetgen/
H A Dctetgenerate.c10 static PetscErrorCode DMPlexInvertCells_CTetgen(PetscInt numCells, PetscInt numCorners, PetscInt ce… in DMPlexInvertCells_CTetgen() argument
12 PetscInt bound = numCells * numCorners, coff; in DMPlexInvertCells_CTetgen()
147 const PetscInt numCells = out->numberoftetrahedra; in DMPlexGenerate_CTetgen() local
165 PetscCall(PetscMalloc1(numCells * numCorners, &cells)); in DMPlexGenerate_CTetgen()
166 for (i = 0; i < numCells * numCorners; i++) cells[i] = (PetscInt)out->tetrahedronlist[i]; in DMPlexGenerate_CTetgen()
169 PetscCall(DMPlexInvertCells_CTetgen(numCells, numCorners, cells)); in DMPlexGenerate_CTetgen()
170 …PetscCall(DMPlexCreateFromCellListPetsc(comm, dim, numCells, numVertices, numCorners, interpolate,… in DMPlexGenerate_CTetgen()
177 …PetscCall(DMUniversalLabelSetLabelValue(universal, *dm, PETSC_TRUE, v + numCells, out->pointmarker… in DMPlexGenerate_CTetgen()
184 …t PetscInt vertices[2] = {out->edgelist[e * 2 + 0] + numCells, out->edgelist[e * 2 + 1] + numCell… in DMPlexGenerate_CTetgen()
196 … {out->trifacelist[f * 3 + 0] + numCells, out->trifacelist[f * 3 + 1] + numCells, out->trifacelist… in DMPlexGenerate_CTetgen()
[all …]
/petsc/src/dm/impls/plex/generators/tetgen/
H A Dtetgenerate.cxx27 static PetscErrorCode DMPlexInvertCells_Tetgen(PetscInt numCells, PetscInt numCorners, PetscInt cel… in DMPlexInvertCells_Tetgen() argument
29 PetscInt bound = numCells * numCorners, coff; in DMPlexInvertCells_Tetgen()
162 const PetscInt numCells = out.numberoftetrahedra; in DMPlexGenerate_Tetgen() local
180 cells = new PetscInt[numCells * numCorners]; in DMPlexGenerate_Tetgen()
181 for (i = 0; i < numCells * numCorners; i++) cells[i] = (PetscInt)out.tetrahedronlist[i]; in DMPlexGenerate_Tetgen()
184 PetscCall(DMPlexInvertCells_Tetgen(numCells, numCorners, cells)); in DMPlexGenerate_Tetgen()
185 …PetscCall(DMPlexCreateFromCellListPetsc(comm, dim, numCells, numVertices, numCorners, interpolate,… in DMPlexGenerate_Tetgen()
190 …PetscCall(DMUniversalLabelSetLabelValue(universal, *dm, PETSC_TRUE, v + numCells, out.pointmarkerl… in DMPlexGenerate_Tetgen()
197 …st PetscInt vertices[2] = {out.edgelist[e * 2 + 0] + numCells, out.edgelist[e * 2 + 1] + numCells… in DMPlexGenerate_Tetgen()
209 … = {out.trifacelist[f * 3 + 0] + numCells, out.trifacelist[f * 3 + 1] + numCells, out.trifacelist[… in DMPlexGenerate_Tetgen()
[all …]
/petsc/src/dm/impls/plex/generators/triangle/
H A Dtrigenerate.c169 const PetscInt numCells = out.numberoftriangles; in DMPlexGenerate_Triangle() local
187 PetscCall(PetscMalloc1(numCells * numCorners, &cells)); in DMPlexGenerate_Triangle()
188 for (i = 0; i < numCells * numCorners; i++) cells[i] = (PetscInt)out.trianglelist[i]; in DMPlexGenerate_Triangle()
190 …PetscCall(DMPlexCreateFromCellListPetsc(comm, dim, numCells, numVertices, numCorners, interpolate,… in DMPlexGenerate_Triangle()
204 if (glabel) PetscCall(DMLabelSetValue(glabel, v + numCells, out.pointmarkerlist[v])); in DMPlexGenerate_Triangle()
210 …st PetscInt vertices[2] = {out.edgelist[e * 2 + 0] + numCells, out.edgelist[e * 2 + 1] + numCells… in DMPlexGenerate_Triangle()
334 const PetscInt numCells = out.numberoftriangles; in DMPlexRefine_Triangle() local
353 PetscCall(PetscMalloc1(numCells * numCorners, &cells)); in DMPlexRefine_Triangle()
354 for (i = 0; i < numCells * numCorners; i++) cells[i] = (PetscInt)out.trianglelist[i]; in DMPlexRefine_Triangle()
357 …PetscCall(DMPlexCreateFromCellListPetsc(comm, dim, numCells, numVertices, numCorners, interpolate,… in DMPlexRefine_Triangle()
[all …]
/petsc/src/dm/impls/plex/
H A Dplexfluent.c172 PetscInt numCells = s->last - s->first + 1; in DMPlexCreateFluent_ReadSection() local
174 PetscCall(PetscMalloc1(numCells, (PetscInt **)&s->data)); in DMPlexCreateFluent_ReadSection()
175 …PetscCall(DMPlexCreateFluent_ReadValues(viewer, s->data, numCells, PETSC_INT, s->index == 2012 ? P… in DMPlexCreateFluent_ReadSection()
684 PetscInt numCells = 0; in DMPlexCreateFluent() local
730 numCells = s.last; in DMPlexCreateFluent()
731 …cCall(PetscInfo((PetscObject)viewer, "CASE: Found number of cells %" PetscInt_FMT "\n", numCells)); in DMPlexCreateFluent()
733 PetscCall(PetscMalloc2(numCells, &cellSizes, numCells, &cellTypes)); in DMPlexCreateFluent()
734 …for (PetscInt c = 0; c < numCells; ++c) PetscCall(GetNumCellFaces(s.nd ? s.nd : (int)((PetscInt *)… in DMPlexCreateFluent()
736 …cObject)viewer, "CASE: Found number of cell faces %" PetscInt_FMT "\n", numCells && s.nd ? cellSiz… in DMPlexCreateFluent()
813 PetscCall(DMPlexSetChart(*dm, 0, numCells + numFaces + numVertices)); in DMPlexCreateFluent()
[all …]
H A Dplexreorder.c71 PetscInt numCells = 0; in DMPlexGetOrdering() local
77 PetscCall(DMPlexCreateNeighborCSR(dm, 0, &numCells, &start, &adjacency)); in DMPlexGetOrdering()
78 PetscCall(PetscMalloc3(numCells, &cperm, numCells, &mask, numCells * 2, &xls)); in DMPlexGetOrdering()
79 if (numCells) { in DMPlexGetOrdering()
81 for (i = 0; i < start[numCells]; ++i) ++adjacency[i]; in DMPlexGetOrdering()
82 for (i = 0; i <= numCells; ++i) ++start[i]; in DMPlexGetOrdering()
83 PetscCall(SPARSEPACKgenrcm(&numCells, start, adjacency, cperm, mask, xls)); in DMPlexGetOrdering()
88 for (c = 0; c < numCells; ++c) --cperm[c]; in DMPlexGetOrdering()
101 …PetscCall(PetscCalloc4(numCells, &sperm, numValues, &values, numValues, &vsize, numValues + 1, &vo… in DMPlexGetOrdering()
111 … numCells, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Label only covers %" PetscInt_FMT " cells !… in DMPlexGetOrdering()
[all …]
H A Dplexcreate.c947 PetscInt i, fStart, fEnd, numCells = 0, numVerts = 0; in DMPlexCreateLineMesh_Internal() local
965 if (rank == 0) numCells = segments; in DMPlexCreateLineMesh_Internal()
969 numPoints[1] = numCells; in DMPlexCreateLineMesh_Internal()
970 …PetscCall(PetscMalloc4(numCells + numVerts, &coneSize, numCells * 2, &cones, numCells + numVerts, … in DMPlexCreateLineMesh_Internal()
971 PetscCall(PetscArrayzero(coneOrientations, numCells + numVerts)); in DMPlexCreateLineMesh_Internal()
972 for (i = 0; i < numCells; ++i) coneSize[i] = 2; in DMPlexCreateLineMesh_Internal()
973 for (i = 0; i < numVerts; ++i) coneSize[numCells + i] = 0; in DMPlexCreateLineMesh_Internal()
974 for (i = 0; i < numCells; ++i) { in DMPlexCreateLineMesh_Internal()
975 cones[2 * i] = numCells + i % numVerts; in DMPlexCreateLineMesh_Internal()
976 cones[2 * i + 1] = numCells + (i + 1) % numVerts; in DMPlexCreateLineMesh_Internal()
[all …]
H A Dplexfem.c2407 PetscInt dim, numCells = cEnd - cStart, c, f; in DMPlexComputeIntegral_Internal() local
2443 PetscCall(ISCreateStride(PETSC_COMM_SELF, numCells, cStart, 1, &cellIS)); in DMPlexComputeIntegral_Internal()
2457 PetscCall(PetscCalloc1(numCells * totDim, &u)); in DMPlexComputeIntegral_Internal()
2458 if (dmAux) PetscCall(PetscMalloc1(numCells * totDimAux, &a)); in DMPlexComputeIntegral_Internal()
2537 numChunks = numCells / (numBatches * batchSize); in DMPlexComputeIntegral_Internal()
2539 Nr = numCells % (numBatches * batchSize); in DMPlexComputeIntegral_Internal()
2540 offset = numCells - Nr; in DMPlexComputeIntegral_Internal()
2544 PetscCall(PetscFEGeomGetChunk(cgeomFEM, offset, numCells, &chunkGeom)); in DMPlexComputeIntegral_Internal()
2546 PetscCall(PetscFEGeomRestoreChunk(cgeomFEM, offset, numCells, &chunkGeom)); in DMPlexComputeIntegral_Internal()
2555 for (c = 0; c < numCells; ++c) { in DMPlexComputeIntegral_Internal()
[all …]
H A Dplexfvm.c91 PetscInt ghost, c, pd, d, numChildren, numCells; in DMPlexReconstructGradients_Internal() local
97 PetscCall(DMPlexGetSupportSize(dm, face, &numCells)); in DMPlexReconstructGradients_Internal()
98 …cCheck(numCells == 2, PETSC_COMM_SELF, PETSC_ERR_PLIB, "facet %" PetscInt_FMT " has %" PetscInt_FM… in DMPlexReconstructGradients_Internal()
/petsc/src/dm/impls/plex/adaptors/pragmatic/
H A Dpragmaticadapt.c22 …PetscInt dim, cStart, cEnd, numCells, c, coff, vStart, vEnd, numVertices, numLocVertices… in DMAdaptMetric_Pragmatic_Plex() local
58 numCells = cEnd - cStart; in DMAdaptMetric_Pragmatic_Plex()
59 if (numCells == 0) { in DMAdaptMetric_Pragmatic_Plex()
66 …numVertices, &y, numVertices, &z, numVertices * PetscSqr(dim), &metric, numCells * maxConeSize, &c… in DMAdaptMetric_Pragmatic_Plex()
69 for (c = 0, coff = 0; c < numCells; ++c) { in DMAdaptMetric_Pragmatic_Plex()
162 …ernalVoid("pragmatic_2d_mpi_init", pragmatic_2d_mpi_init(&numVertices, &numCells, cells, x, y, l2g… in DMAdaptMetric_Pragmatic_Plex()
165 …ernalVoid("pragmatic_3d_mpi_init", pragmatic_3d_mpi_init(&numVertices, &numCells, cells, x, y, z, … in DMAdaptMetric_Pragmatic_Plex()
/petsc/src/dm/dt/fe/interface/
H A Dfegeom.c19 PetscErrorCode PetscFEGeomCreate(PetscQuadrature quad, PetscInt numCells, PetscInt dimEmbed, PetscF… in PetscFEGeomCreate() argument
30 g->numCells = numCells; in PetscFEGeomCreate()
34 N = numCells * Nq; in PetscFEGeomCreate()
37 PetscCall(PetscCalloc2(numCells, &g->face, N * dimEmbed, &g->n)); in PetscFEGeomCreate()
100 (*chunkGeom)->numCells = cEnd - cStart; in PetscFEGeomGetChunk()
265 N = geom->numPoints * geom->numCells; in PetscFEGeomComplete()
/petsc/src/dm/impls/plex/ftn-custom/
H A Dzplexfemf90.c25 PetscInt numCells, totDim, totDimAux = 0; in dmplexgetcellfields_() local
27 *ierr = ISGetLocalSize(*cellIS, &numCells); in dmplexgetcellfields_()
46 …*ierr = F90Array1dCreate((void *)u, MPIU_SCALAR, 1, numCells * totDim, uPtr PETSC_F90_2PTR_PARAM(u… in dmplexgetcellfields_()
48 …*ierr = F90Array1dCreate((void *)u_t, MPIU_SCALAR, 1, locX_t ? numCells * totDim : 0, utPtr PETSC_… in dmplexgetcellfields_()
50 …*ierr = F90Array1dCreate((void *)a, MPIU_SCALAR, 1, locA ? numCells * totDimAux : 0, aPtr PETSC_F9… in dmplexgetcellfields_()
/petsc/src/tao/tutorials/
H A Dex3.c81 PetscInt dim = 2, numCells; in CreateMesh() local
108 PetscCall(VecGetSize(topology, &numCells)); in CreateMesh()
109 numCells = numCells / vertices_per_cell; in CreateMesh()
116 PetscCall(PetscMalloc1(numCells * vertices_per_cell, &cells)); in CreateMesh()
117 for (j = 0; j < numCells * vertices_per_cell; j++) cells[j] = (PetscInt)topo_f[j]; in CreateMesh()
120 …PetscCall(DMPlexCreateFromCellListPetsc(comm, dim, numCells, numVertices, vertices_per_cell, PETSC… in CreateMesh()
/petsc/src/dm/impls/plex/adaptors/mmg/
H A Dmmgadapt.c29 …PetscInt cStart, cEnd, c, numCells, fStart, fEnd, numFaceTags, f, vStart, vEnd, v, numVe… in DMAdaptMetric_Mmg_Plex() local
58 numCells = cEnd - cStart; in DMAdaptMetric_Mmg_Plex()
62 PetscCall(PetscMalloc1(numCells * maxConeSize, &cells)); in DMAdaptMetric_Mmg_Plex()
63 for (c = 0, coff = 0; c < numCells; ++c) { in DMAdaptMetric_Mmg_Plex()
124 PetscCall(PetscCalloc2(numVertices, &verTags, numCells, &cellTags)); in DMAdaptMetric_Mmg_Plex()
169 PetscCallMMG_NONSTANDARD(MMG2D_Set_meshSize, mmg_mesh, numVertices, numCells, 0, numFaceTags); in DMAdaptMetric_Mmg_Plex()
186 …PetscCallMMG_NONSTANDARD(MMG3D_Set_meshSize, mmg_mesh, numVertices, numCells, 0, numFaceTags, 0, 0… in DMAdaptMetric_Mmg_Plex()
/petsc/src/dm/field/impls/ds/
H A Ddmfieldds.c151 PetscInt nq, nc, dim, meshDim, numCells, feDim, i, K = H ? 2 : (D ? 1 : (B ? 0 : -1)); in DMFieldEvaluateFE_DS() local
171 PetscCall(ISGetLocalSize(pointIS, &numCells)); in DMFieldEvaluateFE_DS()
178 for (i = 0; i < numCells; i++) { in DMFieldEvaluateFE_DS()
232 PetscInt c, nFound, numCells, feDim, nc; in DMFieldEvaluate_DS() local
260 PetscCall(PetscSFGetGraph(cellSF, &numCells, &nFound, NULL, &cells)); in DMFieldEvaluate_DS()
264 for (c = 0, gatherSize = 0, gatherMax = 0; c < numCells; c++) { in DMFieldEvaluate_DS()
280 for (c = 0, offset = 0; c < numCells; c++) { in DMFieldEvaluate_DS()
802 numFaces = geom->numCells; in DMFieldComputeFaceData_DS()
837 PetscInt numCells, offset, *cells; in DMFieldComputeFaceData_DS() local
862 for (p = 0, numCells = 0; p < numFaces; p++) { in DMFieldComputeFaceData_DS()
[all …]
/petsc/src/dm/impls/plex/exodusii/
H A Dplexexodusii2.c778 …PetscInt pStart, pEnd, p, cStart, cEnd, numCells, vStart, vEnd, numVertices, eStart, eEn… in DMView_PlexExodusII() local
814 numCells = cEnd - cStart; in DMView_PlexExodusII()
817 …PetscCheck(!(rank && (numCells || numEdges || numVertices)), PETSC_COMM_SELF, PETSC_ERR_SUP, "Writ… in DMView_PlexExodusII()
869 cellsNotInConnectivity = numCells; in DMView_PlexExodusII()
932 …if (num_cs) PetscCallExternal(ex_put_init, exo->exoid, dmName, dim, numNodes, numCells, num_cs, nu… in DMView_PlexExodusII()
1675 int dim = 0, dimEmbed = 0, numVertices = 0, numCells = 0; in DMPlexCreateExodus() local
1686 …PetscCallExternal(ex_get_init, exoid, title, &dimEmbed, &numVertices, &numCells, &num_cs, &num_vs,… in DMPlexCreateExodus()
1692 PetscCall(DMPlexSetChart(*dm, 0, numCells + numVertices)); in DMPlexCreateExodus()
1748 …for (v = numCells; v < numCells + numVertices; ++v) PetscCall(DMPlexSetCellType(*dm, v, DM_POLYTOP… in DMPlexCreateExodus()
1759 … (v_loc = 0; v_loc < num_vertex_per_cell; ++v_loc, ++v) cone[v_loc] = cs_connect[v] + numCells - 1; in DMPlexCreateExodus()
[all …]
/petsc/src/dm/field/impls/da/
H A Ddmfieldda.c308 PetscInt cStart, cEnd, numCells; in DMFieldEvaluateFV_DA() local
332 PetscCall(ISGetLocalSize(cellIS, &numCells)); in DMFieldEvaluateFV_DA()
333 PetscCall(DMGetWorkArray(dm, dim * numCells, MPIU_SCALAR, &points)); in DMFieldEvaluateFV_DA()
337 for (c = 0; c < numCells; c++) { in DMFieldEvaluateFV_DA()
350 …MultilinearEvaluate(dim, dafield->coordRange, nc, dafield->cornerCoeffs, dafield->work, numCells, … in DMFieldEvaluateFV_DA()
351 PetscCall(DMRestoreWorkArray(dm, dim * numCells, MPIU_SCALAR, &points)); in DMFieldEvaluateFV_DA()

12