| /petsc/src/dm/dt/tests/ |
| H A D | ex11.c | 8 PetscInt dim, numFaces; in main() local 20 numFaces = 9; in main() 21 PetscCall(PetscMalloc2(dim * numFaces, &dx, dim * numFaces, &grad)); in main() 24 PetscCall(PetscFVLeastSquaresSetMaxFaces(fvm, numFaces)); in main() 27 PetscCall(PetscFVComputeGradient(fvm, numFaces, dx, grad)); in main()
|
| H A D | ex4.c | 53 PetscInt point = closure[2 * i], numFaces, j; in CheckSymmetry() local 64 numFaces = DMPolytopeTypeGetNumArrangements(ct) / 2; in CheckSymmetry() 66 for (j = -numFaces; j < numFaces; j++) { in CheckSymmetry()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexinterpolate.c | 46 …ces_Internal(DM dm, DMPolytopeType ct, const PetscInt cone[], PetscInt *numFaces, const DMPolytope… in DMPlexGetRawFaces_Internal() argument 74 if (numFaces) *numFaces = 0; in DMPlexGetRawFaces_Internal() 80 if (numFaces) *numFaces = 2; in DMPlexGetRawFaces_Internal() 98 if (numFaces) *numFaces = 2; in DMPlexGetRawFaces_Internal() 116 if (numFaces) *numFaces = 3; in DMPlexGetRawFaces_Internal() 141 if (numFaces) *numFaces = 4; in DMPlexGetRawFaces_Internal() 169 if (numFaces) *numFaces = 4; in DMPlexGetRawFaces_Internal() 198 if (numFaces) *numFaces = 4; in DMPlexGetRawFaces_Internal() 241 if (numFaces) *numFaces = 6; in DMPlexGetRawFaces_Internal() 289 if (numFaces) *numFaces = 5; in DMPlexGetRawFaces_Internal() [all …]
|
| H A D | plexfluent.c | 193 PetscInt numEntries, numFaces, maxsize = 0, offset = 0; in DMPlexCreateFluent_ReadSection() local 214 numFaces = s->last - s->first + 1; in DMPlexCreateFluent_ReadSection() 217 PetscCall(PetscMalloc1(numEntries * numFaces, (PetscInt **)&s->data)); in DMPlexCreateFluent_ReadSection() 219 for (PetscInt f = 0; f < numFaces; f++) { in DMPlexCreateFluent_ReadSection() 225 maxsize = (numFaceVert + 3) * numFaces; in DMPlexCreateFluent_ReadSection() 688 PetscInt numFaces = 0; in DMPlexCreateFluent() local 712 numFaces = PETSC_DETERMINE; in DMPlexCreateFluent() 740 numFaces = (PetscInt)(s.last - s.first + 1); in DMPlexCreateFluent() 741 …ject)viewer, "CASE: Found number of faces %" PetscInt_FMT " face vertices: %d\n", numFaces, s.nd)); in DMPlexCreateFluent() 746 …PetscCheck(numFaces >= 0, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "No header section for facets in F… in DMPlexCreateFluent() [all …]
|
| H A D | plexfem.c | 1076 PetscInt numFaces, f; in DMPlexInsertBoundaryValuesRiemann() local 1080 PetscCall(ISGetLocalSize(faceIS, &numFaces)); in DMPlexInsertBoundaryValuesRiemann() 1082 for (f = 0; f < numFaces; ++f) { in DMPlexInsertBoundaryValuesRiemann() 2767 PetscInt numFaces, face, Nq, field; in DMPlexComputeBdIntegral_Internal() local 2770 PetscCall(ISGetLocalSize(pointIS, &numFaces)); in DMPlexComputeBdIntegral_Internal() 2772 PetscCall(PetscCalloc2(numFaces * totDim, &u, (locA ? (size_t)numFaces * totDimAux : 0), &a)); in DMPlexComputeBdIntegral_Internal() 2774 for (face = 0; face < numFaces; ++face) { in DMPlexComputeBdIntegral_Internal() 2815 numChunks = numFaces / chunkSize; in DMPlexComputeBdIntegral_Internal() 2816 Nr = numFaces % chunkSize; in DMPlexComputeBdIntegral_Internal() 2817 offset = numFaces - Nr; in DMPlexComputeBdIntegral_Internal() [all …]
|
| H A D | plexsubmesh.c | 936 PetscInt numFaces, f, numBdFaces = 0; in DMPlexConstructGhostCells_Internal() local 939 PetscCall(ISGetLocalSize(faceIS, &numFaces)); in DMPlexConstructGhostCells_Internal() 941 for (f = 0; f < numFaces; ++f) { in DMPlexConstructGhostCells_Internal() 972 PetscInt numFaces, f; in DMPlexConstructGhostCells_Internal() local 975 PetscCall(ISGetLocalSize(faceIS, &numFaces)); in DMPlexConstructGhostCells_Internal() 977 for (f = 0; f < numFaces; ++f) { in DMPlexConstructGhostCells_Internal() 999 PetscInt numFaces, f; in DMPlexConstructGhostCells_Internal() local 1002 PetscCall(ISGetLocalSize(faceIS, &numFaces)); in DMPlexConstructGhostCells_Internal() 1004 for (f = 0; f < numFaces; ++f) { in DMPlexConstructGhostCells_Internal() 2466 …dm, DMLabel vertexLabel, PetscInt value, DMLabel subpointMap, PetscInt *numFaces, PetscInt *nFV, D… in DMPlexMarkSubmesh_Uninterpolated() argument [all …]
|
| H A D | plexcreate.c | 705 …PetscInt numFaces = 2 * faces[0] * faces[1] + 2 * faces[1] * faces[2] + 2 * faces[0] * f… in DMPlexCreateBoxSurfaceMesh_Tensor_3D_Internal() local 740 PetscCall(DMPlexSetChart(dm, 0, numFaces + numVertices)); in DMPlexCreateBoxSurfaceMesh_Tensor_3D_Internal() 741 for (f = 0; f < numFaces; ++f) PetscCall(DMPlexSetConeSize(dm, f, 4)); in DMPlexCreateBoxSurfaceMesh_Tensor_3D_Internal() 747 voffset = numFaces + vertices[0] * vertices[1] * (vertices[2] - 1) + vy * vertices[0] + vx; in DMPlexCreateBoxSurfaceMesh_Tensor_3D_Internal() 765 voffset = numFaces + vy * (faces[0] + 1) + vx; in DMPlexCreateBoxSurfaceMesh_Tensor_3D_Internal() 783 voffset = numFaces + vz * vertices[0] * vertices[1] + vx; in DMPlexCreateBoxSurfaceMesh_Tensor_3D_Internal() 801 voffset = numFaces + vz * vertices[0] * vertices[1] + vertices[0] * (vertices[1] - 1) + vx; in DMPlexCreateBoxSurfaceMesh_Tensor_3D_Internal() 819 voffset = numFaces + vz * vertices[0] * vertices[1] + vy * vertices[0]; in DMPlexCreateBoxSurfaceMesh_Tensor_3D_Internal() 837 voffset = numFaces + vz * vertices[0] * vertices[1] + vy * vertices[0] + faces[0]; in DMPlexCreateBoxSurfaceMesh_Tensor_3D_Internal() 858 PetscCall(PetscSectionSetChart(coordSection, numFaces, numFaces + numVertices)); in DMPlexCreateBoxSurfaceMesh_Tensor_3D_Internal() [all …]
|
| H A D | plexgeometry.c | 2807 PetscInt numFaces, f, fOff = 0, p, d; in DMPlexComputeGeometryFVM_3D_Internal() local 2828 PetscCall(DMPlexGetRawFaces_Internal(dm, ct, order, &numFaces, &faceTypes, &faceSizes, &faces)); in DMPlexComputeGeometryFVM_3D_Internal() 2830 for (f = 0; f < numFaces; ++f) { in DMPlexComputeGeometryFVM_3D_Internal() 2899 …PetscCall(DMPlexRestoreRawFaces_Internal(dm, ct, order, &numFaces, &faceTypes, &faceSizes, &faces)… in DMPlexComputeGeometryFVM_3D_Internal() 3220 PetscInt numFaces, usedFaces, f, d; in BuildGradientReconstruction_Internal() local 3230 PetscCall(DMPlexGetConeSize(dm, c, &numFaces)); in BuildGradientReconstruction_Internal() 3232 …numFaces >= dim, PETSC_COMM_SELF, PETSC_ERR_ARG_INCOMP, "Cell %" PetscInt_FMT " has only %" PetscI… in BuildGradientReconstruction_Internal() 3233 for (f = 0, usedFaces = 0; f < numFaces; ++f) { in BuildGradientReconstruction_Internal() 3252 for (f = 0, usedFaces = 0; f < numFaces; ++f) { in BuildGradientReconstruction_Internal() 3337 PetscInt numFaces, f, d, off, ghost = -1; in BuildGradientReconstruction_Internal_Tree() local [all …]
|
| H A D | plexegads.c | 1481 …PetscInt dim = -1, cdim = -1, numCorners = 0, numVertices = 0, numEdges = 0, numFaces = 0, … in DMPlexCreateGeom() local 1521 if (!BIfound) PetscCall(PetscHMapISet(bodyIndexMap, b, numFaces + numEdges + numVertices)); in DMPlexCreateGeom() 1525 if (!BFfound) PetscCall(PetscHMapISet(bodyFaceMap, b, numFaces)); in DMPlexCreateGeom() 1617 numFaces += Nf; in DMPlexCreateGeom() 1627 numPoints = numVertices + numEdges + numFaces; // total number of coordinate points in DMPlexCreateGeom() 5323 …Code DMPlexGetGeomModelBodyFaces(DM dm, PetscGeom body, PetscGeom **faces, PetscInt *numFaces) PeNS in DMPlexGetGeomModelBodyFaces() argument 5338 PetscCall(EGlite_getBodyTopos(body, NULL, FACE, numFaces, faces)); in DMPlexGetGeomModelBodyFaces() 5340 PetscCall(EG_getBodyTopos(body, NULL, FACE, numFaces, faces)); in DMPlexGetGeomModelBodyFaces() 5404 …ModelShellFaces(DM dm, PetscGeom body, PetscGeom shell, PetscGeom **faces, PetscInt *numFaces) PeNS in DMPlexGetGeomModelShellFaces() argument 5419 PetscCall(EGlite_getBodyTopos(body, shell, FACE, numFaces, faces)); in DMPlexGetGeomModelShellFaces() [all …]
|
| H A D | plex.c | 9584 PetscInt numFaces, coneSize, f; in DMPlexCheckFaces() local 9598 … PetscCall(DMPlexGetRawFaces_Internal(dm, ct, closure, &numFaces, &faceTypes, &faceSizes, &faces)); in DMPlexCheckFaces() 9599 …numFaces, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Cell %" PetscInt_FMT " of type %s has %" PetscInt… in DMPlexCheckFaces() 9600 for (f = 0; f < numFaces; ++f) { in DMPlexCheckFaces() 9626 …PetscCall(DMPlexRestoreRawFaces_Internal(dm, ct, closure, &numFaces, &faceTypes, &faceSizes, &face… in DMPlexCheckFaces()
|
| /petsc/src/dm/impls/plex/ftn-custom/ |
| H A D | zplexfemf90.c | 77 PetscInt numFaces = *fEnd - *fStart, totDim; in dmplexgetfacefields_() local 85 …*ierr = F90Array1dCreate((void *)uL, MPIU_SCALAR, 1, numFaces * totDim, uLPtr PETSC_F90_2PTR_PARAM… in dmplexgetfacefields_() 87 …*ierr = F90Array1dCreate((void *)uR, MPIU_SCALAR, 1, numFaces * totDim, uRPtr PETSC_F90_2PTR_PARAM… in dmplexgetfacefields_() 111 …PetscInt numFaces = *fEnd - *fStart, structSize = sizeof(PetscFVFaceGeom) / sizeof(PetscSc… in dmplexgetfacegeometry_() local 115 …*ierr = F90Array1dCreate((void *)g, MPIU_SCALAR, 1, numFaces * structSize, gPtr PETSC_F90_2PTR_PAR… in dmplexgetfacegeometry_() 117 *ierr = F90Array1dCreate((void *)v, MPIU_REAL, 1, numFaces * 2, vPtr PETSC_F90_2PTR_PARAM(vPtrd)); in dmplexgetfacegeometry_()
|
| /petsc/src/dm/impls/plex/generators/ctetgen/ |
| H A D | ctetgenerate.c | 198 PetscInt numFaces; in DMPlexGenerate_CTetgen() local 200 PetscCall(DMPlexGetFullJoin(*dm, 3, vertices, &numFaces, &faces)); in DMPlexGenerate_CTetgen() 201 …PetscCheck(numFaces == 1, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Three vertices must cover only one fac… in DMPlexGenerate_CTetgen() 203 PetscCall(DMPlexRestoreJoin(*dm, 3, vertices, &numFaces, &faces)); in DMPlexGenerate_CTetgen() 473 PetscInt numFaces; in DMPlexRefine_CTetgen() local 475 PetscCall(DMPlexGetFullJoin(*dmRefined, 3, vertices, &numFaces, &faces)); in DMPlexRefine_CTetgen() 476 …PetscCheck(numFaces == 1, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Three vertices must cover only one fac… in DMPlexRefine_CTetgen() 478 PetscCall(DMPlexRestoreJoin(*dmRefined, 3, vertices, &numFaces, &faces)); in DMPlexRefine_CTetgen()
|
| /petsc/src/dm/impls/plex/generators/tetgen/ |
| H A D | tetgenerate.cxx | 211 PetscInt numFaces; in DMPlexGenerate_Tetgen() local 213 PetscCall(DMPlexGetFullJoin(*dm, 3, vertices, &numFaces, &faces)); in DMPlexGenerate_Tetgen() 214 …PetscCheck(numFaces == 1, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Three vertices must cover only one fac… in DMPlexGenerate_Tetgen() 216 PetscCall(DMPlexRestoreJoin(*dm, 3, vertices, &numFaces, &faces)); in DMPlexGenerate_Tetgen() 492 PetscInt numFaces; in DMPlexRefine_Tetgen() local 494 PetscCall(DMPlexGetFullJoin(*dmRefined, 3, vertices, &numFaces, &faces)); in DMPlexRefine_Tetgen() 495 …PetscCheck(numFaces == 1, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Three vertices must cover only one fac… in DMPlexRefine_Tetgen() 497 PetscCall(DMPlexRestoreJoin(*dmRefined, 3, vertices, &numFaces, &faces)); in DMPlexRefine_Tetgen()
|
| /petsc/src/dm/dt/fv/interface/ |
| H A D | fv.c | 1700 PetscErrorCode PetscFVComputeGradient(PetscFV fvm, PetscInt numFaces, PetscScalar dx[], PetscScalar… in PetscFVComputeGradient() argument 1704 PetscTryTypeMethod(fvm, computegradient, numFaces, dx, grad); in PetscFVComputeGradient() 1891 static PetscErrorCode PetscFVComputeGradient_Upwind(PetscFV fv, PetscInt numFaces, const PetscScala… in PetscFVComputeGradient_Upwind() argument 1897 for (PetscInt f = 0; f < numFaces; ++f) { in PetscFVComputeGradient_Upwind() 2112 PetscInt numFaces, f; 2115 PetscCall(DMPlexGetConeSize(dm, cell, &numFaces)); 2117 for (f = 0; f < numFaces; ++f) { 2151 static PetscErrorCode PetscFVComputeGradient_LeastSquares(PetscFV fvm, PetscInt numFaces, const Pet… in PetscFVComputeGradient_LeastSquares() argument 2159 if (numFaces > maxFaces) { in PetscFVComputeGradient_LeastSquares() 2161 …RANGE, "Number of input faces %" PetscInt_FMT " > %" PetscInt_FMT " maxfaces", numFaces, maxFaces); in PetscFVComputeGradient_LeastSquares() [all …]
|
| /petsc/src/dm/impls/plex/adaptors/mmg/ |
| H A D | mmgadapt.c | 257 PetscInt numCoveredPoints, numFaces = 0, facePoints[3]; in DMAdaptMetric_Mmg_Plex() local 264 numFaces++; in DMAdaptMetric_Mmg_Plex() 268 …heck(numFaces == 1, comm, PETSC_ERR_ARG_OUTOFRANGE, "%" PetscInt_FMT " vertices cannot define more… in DMAdaptMetric_Mmg_Plex()
|
| /petsc/src/dm/dt/dualspace/impls/lagrange/tests/ |
| H A D | ex1.c | 274 PetscInt numFaces, o; in testLagrange() local 280 numFaces = DMPolytopeTypeGetNumArrangements(ct) / 2; in testLagrange() 282 for (o = -numFaces; o < numFaces; ++o) { in testLagrange()
|
| /petsc/src/dm/field/impls/ds/ |
| H A D | dmfieldds.c | 789 PetscInt p, dim, dE, numFaces, Nq; in DMFieldComputeFaceData_DS() local 802 numFaces = geom->numCells; in DMFieldComputeFaceData_DS() 805 for (p = 0; p < numFaces; p++) { in DMFieldComputeFaceData_DS() 862 for (p = 0, numCells = 0; p < numFaces; p++) { in DMFieldComputeFaceData_DS() 873 for (p = 0, offset = 0; p < numFaces; p++) { in DMFieldComputeFaceData_DS() 884 for (p = 0, offset = 0; p < numFaces; p++) { in DMFieldComputeFaceData_DS() 934 PetscCall(PetscMalloc2(numFaces, &co, coneSize, &counts)); in DMFieldComputeFaceData_DS() 941 for (p = 0; p < numFaces; p++) { /* record the orientation of the facet wrt the support cells */ in DMFieldComputeFaceData_DS() 988 for (p = 0; p < numFaces; p++) { in DMFieldComputeFaceData_DS() 1088 for (p = 0, numCells = 0; p < numFaces; p++) { in DMFieldComputeFaceData_DS() [all …]
|
| /petsc/src/dm/impls/plex/exodusii/ |
| H A D | plexexodusii2.c | 778 …End, numCells, vStart, vEnd, numVertices, eStart, eEnd, numEdges, fStart, fEnd, numFaces, numNodes; in DMView_PlexExodusII() local 848 numFaces = fEnd - fStart; in DMView_PlexExodusII() 850 numFaces = 0; in DMView_PlexExodusII() 906 numNodes += numFaces; in DMView_PlexExodusII() 996 if (nodes[cs][2] == 0) connect[i + off] -= numFaces; in DMView_PlexExodusII() 1878 PetscInt faceSize = fs_vertex_count_list[f], numFaces; in DMPlexCreateExodus() local 1883 PetscCall(DMPlexGetFullJoin(*dm, faceSize, faceVertices, &numFaces, &faces)); in DMPlexCreateExodus() 1884 …PetscCheck(numFaces == 1, comm, PETSC_ERR_ARG_WRONG, "Invalid ExodusII side %d in set %d maps to %… in DMPlexCreateExodus() 1889 PetscCall(DMPlexRestoreJoin(*dm, faceSize, faceVertices, &numFaces, &faces)); in DMPlexCreateExodus()
|
| /petsc/src/dm/impls/plex/adaptors/parmmg/ |
| H A D | parmmgadapt.c | 422 PetscInt numCoveredPoints, numFaces = 0, facePoints[3]; in DMAdaptMetric_ParMmg_Plex() local 434 numFaces++; in DMAdaptMetric_ParMmg_Plex() 438 …heck(numFaces == 1, comm, PETSC_ERR_ARG_OUTOFRANGE, "%" PetscInt_FMT " vertices cannot define more… in DMAdaptMetric_ParMmg_Plex()
|
| /petsc/src/dm/dt/dualspace/impls/lagrange/ |
| H A D | dspacelagrange.c | 2681 PetscInt numFaces; in PetscDualSpaceGetSymmetries_Lagrange() local 2692 numFaces = DMPolytopeTypeGetNumArrangements(ct) / 2; in PetscDualSpaceGetSymmetries_Lagrange() 2708 lag->numSelfSym = 2 * numFaces; in PetscDualSpaceGetSymmetries_Lagrange() 2709 lag->selfSymOff = numFaces; in PetscDualSpaceGetSymmetries_Lagrange() 2710 PetscCall(PetscCalloc1(2 * numFaces, &cellSymperms)); in PetscDualSpaceGetSymmetries_Lagrange() 2711 PetscCall(PetscCalloc1(2 * numFaces, &cellSymflips)); in PetscDualSpaceGetSymmetries_Lagrange() 2713 symperms[0] = &cellSymperms[numFaces]; in PetscDualSpaceGetSymmetries_Lagrange() 2714 symflips[0] = &cellSymflips[numFaces]; in PetscDualSpaceGetSymmetries_Lagrange() 2717 …for (ornt = -numFaces; ornt < numFaces; ornt++) { /* for every symmetry, compute the symmetry matr… in PetscDualSpaceGetSymmetries_Lagrange() 2772 for (ornt = -numFaces; ornt < numFaces; ornt++) in PetscDualSpaceGetSymmetries_Lagrange() [all …]
|
| /petsc/src/ts/tutorials/ |
| H A D | ex11.c | 777 PetscInt numFaces; in CreateMassMatrix() local 779 PetscCall(DMPlexGetSupportSize(dmMass, v, &numFaces)); in CreateMassMatrix() 780 PetscCall(PetscSectionSetDof(sectionMass, v, numFaces * numFaces)); in CreateMassMatrix() 798 PetscInt numFaces, sides[2], f, g; in CreateMassMatrix() local 801 PetscCall(DMPlexGetSupportSize(dmMass, v, &numFaces)); in CreateMassMatrix() 803 for (f = 0; f < numFaces; ++f) { in CreateMassMatrix() 806 for (g = 0; g < numFaces; ++g) { in CreateMassMatrix() 818 m[f * numFaces + g] = Dot2Real(fgA->normal, fgB->normal) * area * 0.5; in CreateMassMatrix()
|
| /petsc/src/dm/dt/fe/interface/ |
| H A D | fe.c | 935 PetscInt dim, numFaces, f; in PetscFEGetFaceCentroidTabulation() local 940 PetscCall(DMPlexGetConeSize(dm, 0, &numFaces)); in PetscFEGetFaceCentroidTabulation() 942 PetscCall(PetscMalloc1(numFaces * dim, ¢roids)); in PetscFEGetFaceCentroidTabulation() 943 …for (f = 0; f < numFaces; ++f) PetscCall(DMPlexComputeCellGeometryFVM(dm, cone[f], NULL, ¢roid… in PetscFEGetFaceCentroidTabulation() 944 PetscCall(PetscFECreateTabulation(fem, 1, numFaces, centroids, 0, &fem->Tc)); in PetscFEGetFaceCentroidTabulation()
|