Home
last modified time | relevance | path

Searched refs:geom (Results 1 – 21 of 21) sorted by relevance

/petsc/src/dm/dt/fe/interface/
H A Dfegeom.c19 …tscQuadrature quad, PetscInt numCells, PetscInt dimEmbed, PetscFEGeomMode mode, PetscFEGeom **geom) in PetscFEGeomCreate() argument
41 *geom = g; in PetscFEGeomCreate()
55 PetscErrorCode PetscFEGeomDestroy(PetscFEGeom **geom) in PetscFEGeomDestroy() argument
58 if (!*geom) PetscFunctionReturn(PETSC_SUCCESS); in PetscFEGeomDestroy()
59 PetscCall(PetscFree3((*geom)->v, (*geom)->J, (*geom)->detJ)); in PetscFEGeomDestroy()
60 PetscCall(PetscFree((*geom)->invJ)); in PetscFEGeomDestroy()
61 PetscCall(PetscFree2((*geom)->face, (*geom)->n)); in PetscFEGeomDestroy()
62 …etscCall(PetscFree6((*geom)->suppJ[0], (*geom)->suppJ[1], (*geom)->suppInvJ[0], (*geom)->suppInvJ[… in PetscFEGeomDestroy()
63 PetscCall(PetscFree(*geom)); in PetscFEGeomDestroy()
85 PetscErrorCode PetscFEGeomGetChunk(PetscFEGeom *geom, PetscInt cStart, PetscInt cEnd, PetscFEGeom *… in PetscFEGeomGetChunk() argument
[all …]
H A Dfe.c1435 … PetscInt, const PetscScalar[], PetscScalar[]), PetscInt Ne, PetscFEGeom *geom, const PetscScalar … in PetscFEIntegrateBd() argument
1442 …if (fe->ops->integratebd) PetscCall((*fe->ops->integratebd)(prob, field, obj_func, Ne, geom, coeff… in PetscFEIntegrateBd()
/petsc/src/ksp/ksp/tutorials/amrex/
H A DinitEB.cxx25 EB2::Build(gshop, geom.back(), max_level, max_level + max_coarsening_level); in initializeEB()
31 EB2::Build(gshop, geom.back(), max_level, max_level + max_coarsening_level); in initializeEB()
38 EB2::Build(gshop, geom.back(), max_level, max_level + max_coarsening_level); in initializeEB()
42 EB2::Build(gshop, geom.back(), max_level, max_level + max_coarsening_level); in initializeEB()
44 EB2::Build(geom.back(), max_level, max_level + max_coarsening_level); in initializeEB()
H A DMyTest.cxx38 if (geom[0].isPeriodic(idim)) { in solve()
50 MLEBABecLap mleb(geom, grids, dmap, info, amrex::GetVecOfConstPtrs(factory)); in solve()
91 …_name, max_level + 1, amrex::GetVecOfConstPtrs(plotmf), {"phi", "vfrac"}, geom, 0.0, Vector<int>(m… in writePlotfile()
135 geom.resize(nlevels); in initGrids()
144 geom[ilev].define(domain); in initGrids()
171 const EB2::Level &eb_level = eb_is.getLevel(geom[ilev]); in initData()
172 …factory[ilev].reset(new EBFArrayBoxFactory(eb_level, geom[ilev], grids[ilev], dmap[ilev], {2, 2, 2… in initData()
188 const auto dx = geom[ilev].CellSizeArray(); in initData()
H A DMyTest.H46 amrex::Vector<amrex::Geometry> geom; variable
/petsc/src/dm/field/impls/shell/
H A Ddmfieldshell.c35 PetscFEGeom *geom; in DMFieldShellEvaluateFEDefault() local
43 PetscCall(DMFieldCreateFEGeom(coordField, pointIS, quad, PETSC_FEGEOM_BASIC, &geom)); in DMFieldShellEvaluateFEDefault()
48 for (p = 0; p < numPoints * dimC * Nq; p++) pfArray[p] = geom->v[p]; in DMFieldShellEvaluateFEDefault()
66 …for (j = 0; j < dimC; j++) sD[(p * Nc + q) * dimC + i] += geom->J[(p * dimC + j) * dimC + i] * d[j… in DMFieldShellEvaluateFEDefault()
83 …for (j = 0; j < dimC; j++) rD[(p * Nc + q) * dimC + i] += geom->J[(p * dimC + j) * dimC + i] * d[j… in DMFieldShellEvaluateFEDefault()
107 … dimC; l++) sH[((p * Nc + q) * dimC + i) * dimC + j] += geom->J[(p * dimC + k) * dimC + i] * geom-… in DMFieldShellEvaluateFEDefault()
130 … dimC; l++) rH[((p * Nc + q) * dimC + i) * dimC + j] += geom->J[(p * dimC + k) * dimC + i] * geom-… in DMFieldShellEvaluateFEDefault()
140 PetscCall(PetscFEGeomDestroy(&geom)); in DMFieldShellEvaluateFEDefault()
148 PetscFEGeom *geom; in DMFieldShellEvaluateFVDefault() local
164 PetscCall(DMFieldCreateFEGeom(coordField, pointIS, quad, PETSC_FEGEOM_BASIC, &geom)); in DMFieldShellEvaluateFVDefault()
[all …]
/petsc/src/dm/impls/plex/tests/
H A Dex36.c75 const PetscFVCellGeom *geom; in dm_view_geometry() local
77 geom = (PetscFVCellGeom *)&cell_array[offset]; in dm_view_geometry()
78 …%g,%g vol %g\n", rank, c, (double)geom->centroid[0], (double)geom->centroid[1], (double)geom->cent… in dm_view_geometry()
/petsc/src/dm/field/impls/ds/
H A Ddmfieldds.c457 PetscFEGeom *geom; in DMFieldEvaluateFV_DS() local
487 PetscCall(DMFieldCreateFEGeom(coordField, pointIS, quad, PETSC_FEGEOM_BASIC, &geom)); in DMFieldEvaluateFV_DS()
507 vol += geom->detJ[i * Nq + k] * weights[k]; in DMFieldEvaluateFV_DS()
508 …for (j = 0; j < Nc; j++) sB[i * Nc + j] += geom->detJ[i * Nq + k] * weights[k] * sqB[(i * Nq + k) … in DMFieldEvaluateFV_DS()
521 vol += geom->detJ[i * Nq + k] * weights[k]; in DMFieldEvaluateFV_DS()
540 vol += geom->detJ[i * Nq + k] * weights[k]; in DMFieldEvaluateFV_DS()
545 …for (m = 0; m < dim; m++) pD[l] += geom->invJ[((i * Nq + k) * dimC + m) * dimC + l] * sqD[((i * Nq… in DMFieldEvaluateFV_DS()
547 …for (l = 0; l < dimC; l++) sD[(i * Nc + j) * dimC + l] += geom->detJ[i * Nq + k] * weights[k] * pD… in DMFieldEvaluateFV_DS()
561 vol += geom->detJ[i * Nq + k] * weights[k]; in DMFieldEvaluateFV_DS()
566 …for (m = 0; m < dim; m++) pD[l] += geom->invJ[((i * Nq + k) * dimC + m) * dimC + l] * rqD[((i * Nq… in DMFieldEvaluateFV_DS()
[all …]
/petsc/src/dm/interface/
H A Ddmceed.c140 CeedQFunctionUser geom = NULL; in DMCeedCreateGeometry() local
165 geom = Geometry2D; in DMCeedCreateGeometry()
169 geom = Geometry3D; in DMCeedCreateGeometry()
173 PetscCallCEED(CeedQFunctionCreateInterior(ceed, 1, geom, geomName, &sd->qf)); in DMCeedCreateGeometry()
193 …if (coarse->dmceed) PetscCall(DMCeedCreate(fine, coarse->dmceed->geom ? PETSC_TRUE : PETSC_FALSE, … in DMRefineHook_Ceed()
221 PetscCall(DMCeedCreateGeometry(cdm, cellIS, &Nqdata, &sd->erq, &sd->qd, &sd->geom)); in DMCeedCreate_Internal()
224 if (sd->geom) { in DMCeedCreate_Internal()
228 PetscCallCEED(CeedBasisGetNumQuadraturePoints(sd->geom->basis, &Nqx)); in DMCeedCreate_Internal()
346 PetscCall(DMCeedCreateGeometryFVM(cdm, faceIS, &Nqdata, &sd->erq, &sd->qd, &sd->geom)); in DMCeedCreateFVM_Internal()
418 PetscCall(DMCeedDestroy(&p->geom)); in DMCeedDestroy()
[all …]
/petsc/src/dm/dt/fe/tests/
H A Dex2.c94 …Geom(IS cellIS, DMField coordField, PetscQuadrature quad, PetscFEGeomMode mode, PetscFEGeom **geom) in CellRangeGetFEGeom() argument
105 PetscCall(PetscContainerGetPointer(container, geom)); in CellRangeGetFEGeom()
107 PetscCall(DMFieldCreateFEGeom(coordField, cellIS, quad, mode, geom)); in CellRangeGetFEGeom()
108 …PetscCall(PetscObjectContainerCompose((PetscObject)cellIS, composeStr, *geom, PetscFEGeomDestroy_V… in CellRangeGetFEGeom()
113 …FEGeom(IS cellIS, DMField coordField, PetscQuadrature quad, PetscBool faceData, PetscFEGeom **geom) in CellRangeRestoreFEGeom() argument
116 *geom = NULL; in CellRangeRestoreFEGeom()
267 PetscFEGeom *geom = affineGeom ? affineGeom : geoms[f]; in TestIntegration() local
274 PetscCall(PetscFEGeomGetChunk(geom, cS, cE, &chunkGeom)); in TestIntegration()
H A Dex1.c103 …Geom(IS cellIS, DMField coordField, PetscQuadrature quad, PetscFEGeomMode mode, PetscFEGeom **geom) in CellRangeGetFEGeom() argument
114 PetscCall(PetscContainerGetPointer(container, geom)); in CellRangeGetFEGeom()
116 PetscCall(DMFieldCreateFEGeom(coordField, cellIS, quad, mode, geom)); in CellRangeGetFEGeom()
117 …PetscCall(PetscObjectContainerCompose((PetscObject)cellIS, composeStr, *geom, PetscFEGeomDestroy_V… in CellRangeGetFEGeom()
122 …FEGeom(IS cellIS, DMField coordField, PetscQuadrature quad, PetscBool faceData, PetscFEGeom **geom) in CellRangeRestoreFEGeom() argument
125 *geom = NULL; in CellRangeRestoreFEGeom()
219 PetscFEGeom *geom = affineGeom ? affineGeom : geoms[f]; in TestIntegration() local
226 PetscCall(PetscFEGeomGetChunk(geom, cS, cE, &chunkGeom)); in TestIntegration()
/petsc/src/dm/impls/plex/
H A Dplexegads.c485 ego geom, body, obj; in DMSnapToGeomModel_EGADS_Internal() local
501 PetscCall(EGlite_getTopology(model, &geom, &oclass, &mtype, NULL, &Nb, &bodies, &senses)); in DMSnapToGeomModel_EGADS_Internal()
503 PetscCall(EG_getTopology(model, &geom, &oclass, &mtype, NULL, &Nb, &bodies, &senses)); in DMSnapToGeomModel_EGADS_Internal()
617 ego geom, *bodies, *nobjs, *mobjs, *lobjs, *shobjs, *fobjs, *eobjs; in DMPlexGeomPrintModel_Internal() local
624 PetscCall(EGlite_getTopology(model, &geom, &oclass, &mtype, NULL, &Nb, &bodies, &senses)); in DMPlexGeomPrintModel_Internal()
626 PetscCall(EG_getTopology(model, &geom, &oclass, &mtype, NULL, &Nb, &bodies, &senses)); in DMPlexGeomPrintModel_Internal()
694 PetscCall(EGlite_getTopology(body, &geom, &oclass, &mtype, NULL, &Nsh, &shobjs, &shsenses)); in DMPlexGeomPrintModel_Internal()
696 PetscCall(EG_getTopology(body, &geom, &oclass, &mtype, NULL, &Nsh, &shobjs, &shsenses)); in DMPlexGeomPrintModel_Internal()
712 PetscCall(EGlite_getTopology(shell, &geom, &oclass, &mtype, NULL, &Nf, &fobjs, &fsenses)); in DMPlexGeomPrintModel_Internal()
714 PetscCall(EG_getTopology(shell, &geom, &oclass, &mtype, NULL, &Nf, &fobjs, &fsenses)); in DMPlexGeomPrintModel_Internal()
[all …]
H A Dplexfem.c52 PetscFEGeom *geom = *(PetscFEGeom **)ctx; in PetscContainerCtxDestroy_PetscFEGeom() local
55 PetscCall(PetscFEGeomDestroy(&geom)); in PetscContainerCtxDestroy_PetscFEGeom()
59 …eom(DMField coordField, IS pointIS, PetscQuadrature quad, PetscFEGeomMode mode, PetscFEGeom **geom) in DMPlexGetFEGeom() argument
70 PetscCall(PetscContainerGetPointer(container, geom)); in DMPlexGetFEGeom()
72 PetscCall(DMFieldCreateFEGeom(coordField, pointIS, quad, mode, geom)); in DMPlexGetFEGeom()
74 PetscCall(PetscContainerSetPointer(container, (void *)*geom)); in DMPlexGetFEGeom()
82 …eom(DMField coordField, IS pointIS, PetscQuadrature quad, PetscFEGeomMode mode, PetscFEGeom **geom) in DMPlexRestoreFEGeom() argument
85 *geom = NULL; in DMPlexRestoreFEGeom()
4359 …eom(DMField coordField, IS pointIS, PetscQuadrature quad, PetscFEGeomMode mode, PetscFEGeom **geom) in DMSNESGetFEGeom() argument
4370 PetscCall(PetscContainerGetPointer(container, geom)); in DMSNESGetFEGeom()
[all …]
/petsc/src/ts/tutorials/
H A Dex11.h230 const CeedScalar *xL = in[0], *xR = in[1], *geom = in[2]; in PhysicsRiemann_SW_Rusanov_CEED() local
243 const CeedScalar n[2] = {geom[i + Q * 0], geom[i + Q * 1]}; in PhysicsRiemann_SW_Rusanov_CEED()
256 cL[i + Q * j] = -flux[j] / geom[i + Q * 2]; in PhysicsRiemann_SW_Rusanov_CEED()
257 cR[i + Q * j] = flux[j] / geom[i + Q * 3]; in PhysicsRiemann_SW_Rusanov_CEED()
261 … j < DIM + 1; ++j) PetscPrintf(PETSC_COMM_SELF, " | %g | (%g)\n", cL[i + Q * j], geom[i + Q * 2]); in PhysicsRiemann_SW_Rusanov_CEED()
263 … j < DIM + 1; ++j) PetscPrintf(PETSC_COMM_SELF, " | %g | (%g)\n", cR[i + Q * j], geom[i + Q * 3]); in PhysicsRiemann_SW_Rusanov_CEED()
842 const CeedScalar *xL = in[0], *xR = in[1], *geom = in[2]; in PhysicsRiemann_Euler_Godunov_CEED() local
852 const CeedScalar n[DIM] = {geom[i + Q * 0], geom[i + Q * 1]}; in PhysicsRiemann_Euler_Godunov_CEED()
865 cL[i + Q * j] = -flux[j] / geom[i + Q * 2]; in PhysicsRiemann_Euler_Godunov_CEED()
866 cR[i + Q * j] = flux[j] / geom[i + Q * 3]; in PhysicsRiemann_Euler_Godunov_CEED()
[all …]
/petsc/lib/petsc/bin/
H A Dpetsc_gen_xdmf.py248 geom = h5['viz']['geometry']
251 geom = h5['geometry']
264 vertices = geom['vertices']
/petsc/include/
H A Dpetscdmceed.h40 DMCeed geom; // Operator computing geometric data at quadrature points member
/petsc/src/dm/impls/plex/generators/tetgen/
H A Dtetgenerate.cxx230 ego model, geom; in DMPlexGenerate_Tetgen() local
239 PetscCall(EG_getTopology(model, &geom, &oclass, &mtype, nullptr, &Nb, &bodies, &senses)); in DMPlexGenerate_Tetgen()
244 … PetscCall(EGlite_getTopology(model, &geom, &oclass, &mtype, nullptr, &Nb, &bodies, &senses)); in DMPlexGenerate_Tetgen()
509 ego model, geom; in DMPlexRefine_Tetgen() local
518 PetscCall(EG_getTopology(model, &geom, &oclass, &mtype, nullptr, &Nb, &bodies, &senses)); in DMPlexRefine_Tetgen()
523 … PetscCall(EGlite_getTopology(model, &geom, &oclass, &mtype, nullptr, &Nb, &bodies, &senses)); in DMPlexRefine_Tetgen()
/petsc/src/dm/impls/forest/p4est/
H A Dpforest.h172 p4est_geometry_t *geom; member
206 static void GeometryMapping_pforest(p4est_geometry_t *geom, p4est_topidx_t which_tree, const double… in GeometryMapping_pforest() argument
208 DM_Forest_geometry_pforest *geom_pforest = (DM_Forest_geometry_pforest *)geom->user; in GeometryMapping_pforest()
224 static void GeometryDestroy_pforest(p4est_geometry_t *geom) in GeometryDestroy_pforest() argument
226 DM_Forest_geometry_pforest *geom_pforest = (DM_Forest_geometry_pforest *)geom->user; in GeometryDestroy_pforest()
230 ierr = PetscFree(geom->user); in GeometryDestroy_pforest()
232 ierr = PetscFree(geom); in GeometryDestroy_pforest()
245 if ((*topo)->geom) PetscCallP4est(p4est_geometry_destroy, ((*topo)->geom)); in DMFTopologyDestroy_pforest()
278 (*topo)->geom = NULL; in DMFTopologyCreateBrick_pforest()
330 (*topo)->geom = NULL; in DMFTopologyCreate_pforest()
[all …]
/petsc/src/dm/field/interface/
H A Ddmfield.c439 …teFEGeom(DMField field, IS pointIS, PetscQuadrature quad, PetscFEGeomMode mode, PetscFEGeom **geom) in DMFieldCreateFEGeom() argument
526 *geom = g; in DMFieldCreateFEGeom()
/petsc/src/dm/impls/plex/generators/ctetgen/
H A Dctetgenerate.c214 ego model, geom; in DMPlexGenerate_CTetgen() local
221 PetscCall(EG_getTopology(model, &geom, &oclass, &mtype, NULL, &Nb, &bodies, &senses)); in DMPlexGenerate_CTetgen()
489 ego model, geom; in DMPlexRefine_CTetgen() local
496 PetscCall(EG_getTopology(model, &geom, &oclass, &mtype, NULL, &Nb, &bodies, &senses)); in DMPlexRefine_CTetgen()
/petsc/share/petsc/datafiles/meshes/
H A Dtestcase3D.cas132 (virtualboundary/geom ())
510 (spark/user-geom? (#f))
511 (spark/time-dep-geom? (#t))
512 (spark/fixed-geom? (#f))
589 (domains (((1 geom-domain mixture) (children) (material . air))))
713 (heatxc/geom ())
1907 (dpm/volume-injection/show-bound-geom? #f)
5255 (species/edc-geom-mean-tmix? #f)
6186 (udf/spark-geom-fcn ("none"))
6643 (dynamesh/auxgeom_prefix "geom::")
[all …]