Home
last modified time | relevance | path

Searched refs:Nq (Results 1 – 25 of 62) sorted by relevance

123

/petsc/src/dm/dt/fe/interface/
H A Dfegeom.c22 PetscInt dim, Nq, N; in PetscFEGeomCreate() local
26 PetscCall(PetscQuadratureGetData(quad, &dim, NULL, &Nq, &p, NULL)); in PetscFEGeomCreate()
31 g->numPoints = Nq; in PetscFEGeomCreate()
34 N = numCells * Nq; in PetscFEGeomCreate()
87 PetscInt Nq; in PetscFEGeomGetChunk() local
94 Nq = geom->numPoints; in PetscFEGeomGetChunk()
102 (*chunkGeom)->v = PetscSafePointerPlusOffset(geom->v, Nq * dE * cStart); in PetscFEGeomGetChunk()
103 (*chunkGeom)->J = PetscSafePointerPlusOffset(geom->J, Nq * dE * dE * cStart); in PetscFEGeomGetChunk()
104 (*chunkGeom)->invJ = PetscSafePointerPlusOffset(geom->invJ, Nq * dE * dE * cStart); in PetscFEGeomGetChunk()
105 (*chunkGeom)->detJ = PetscSafePointerPlusOffset(geom->detJ, Nq * cStart); in PetscFEGeomGetChunk()
[all …]
H A Dfe.c2368 const PetscInt Nq = T[f]->Np; in PetscFEEvaluateFieldJets_Internal() local
2371 const PetscReal *Bq = &T[f]->T[0][(r * Nq + q) * Nbf * Ncf]; in PetscFEEvaluateFieldJets_Internal()
2372 const PetscReal *Dq = &T[f]->T[1][(r * Nq + q) * Nbf * Ncf * cdim]; in PetscFEEvaluateFieldJets_Internal()
2373 const PetscReal *Hq = k > 1 ? &T[f]->T[2][(r * Nq + q) * Nbf * Ncf * cdim * cdim] : NULL; in PetscFEEvaluateFieldJets_Internal()
2435 const PetscInt Nq = T->Np; in PetscFEEvaluateFieldJets_Hybrid_Internal() local
2442 const PetscReal *Bq = &T->T[0][(r * Nq + q) * Nbf * Ncf]; in PetscFEEvaluateFieldJets_Hybrid_Internal()
2443 const PetscReal *Dq = &T->T[1][(r * Nq + q) * Nbf * Ncf * dEt]; in PetscFEEvaluateFieldJets_Hybrid_Internal()
2506 const PetscInt Nq = T->Np; in PetscFEUpdateElementVec_Internal() local
2509 const PetscReal *basis = &T->T[0][r * Nq * Nb * Nc]; in PetscFEUpdateElementVec_Internal()
2510 const PetscReal *basisDer = &T->T[1][r * Nq * Nb * Nc * dEt]; in PetscFEUpdateElementVec_Internal()
[all …]
/petsc/src/dm/field/impls/ds/
H A Ddmfieldds.c458 PetscInt Nq, Nc, dimC, qNc, N; in DMFieldEvaluateFV_DS() local
488 PetscCall(PetscQuadratureGetData(quad, NULL, &qNc, &Nq, NULL, &weights)); in DMFieldEvaluateFV_DS()
490 N = numPoints * Nq * Nc; in DMFieldEvaluateFV_DS()
506 for (k = 0; k < Nq; k++) { 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()
520 for (k = 0; k < Nq; k++) { in DMFieldEvaluateFV_DS()
521 vol += geom->detJ[i * Nq + k] * weights[k]; in DMFieldEvaluateFV_DS()
522 for (j = 0; j < Nc; j++) rB[i * Nc + j] += weights[k] * rqB[(i * Nq + k) * Nc + j]; in DMFieldEvaluateFV_DS()
539 for (k = 0; k < Nq; k++) { in DMFieldEvaluateFV_DS()
[all …]
/petsc/src/dm/field/impls/shell/
H A Ddmfieldshell.c37 PetscInt dimC, dim, numPoints, Nq, p, Nc; in DMFieldShellEvaluateFEDefault() local
45 PetscCall(PetscQuadratureGetData(quad, &dim, NULL, &Nq, NULL, NULL)); in DMFieldShellEvaluateFEDefault()
47 PetscCall(PetscMalloc1(dimC * Nq * numPoints, &pfArray)); in DMFieldShellEvaluateFEDefault()
48 for (p = 0; p < numPoints * dimC * Nq; p++) pfArray[p] = geom->v[p]; in DMFieldShellEvaluateFEDefault()
49 …PetscCall(VecCreateMPIWithArray(PetscObjectComm((PetscObject)pointIS), dimC, dimC * Nq * numPoints… in DMFieldShellEvaluateFEDefault()
57 for (p = 0; p < numPoints * Nq; p++) { in DMFieldShellEvaluateFEDefault()
74 for (p = 0; p < numPoints * Nq; p++) { in DMFieldShellEvaluateFEDefault()
94 for (p = 0; p < numPoints * Nq; p++) { in DMFieldShellEvaluateFEDefault()
117 for (p = 0; p < numPoints * Nq; p++) { in DMFieldShellEvaluateFEDefault()
150 PetscInt dimC, dim, numPoints, Nq, p; in DMFieldShellEvaluateFVDefault() local
[all …]
/petsc/src/dm/interface/
H A Ddmceed.c144 CeedInt Nq; in DMCeedCreateGeometry() local
157 PetscCall(CeedBasisGetNumQuadraturePoints(sd->basis, &Nq)); in DMCeedCreateGeometry()
161 …PetscCallCEED(CeedElemRestrictionCreateStrided(ceed, Ncell, Nq, *Nqdata, Ncell * Nq * (*Nqdata), C… in DMCeedCreateGeometry()
204 CeedInt Nq; in DMCeedCreate_Internal() local
214 PetscCall(CeedBasisGetNumQuadraturePoints(sd->basis, &Nq)); in DMCeedCreate_Internal()
229 …Nq, PetscObjectComm((PetscObject)dm), PETSC_ERR_ARG_INCOMP, "Number of qpoints for solution %" Cee… in DMCeedCreate_Internal()
280 PetscInt dim, cdim, fStart, fEnd, Nface, Nq = 1; in DMCeedCreateGeometryFVM() local
292 strides[1] = Nq; in DMCeedCreateGeometryFVM()
293 strides[2] = Nq * (*Nqdata); in DMCeedCreateGeometryFVM()
294 …PetscCallCEED(CeedElemRestrictionCreateStrided(ceed, Nface, Nq, *Nqdata, Nface * Nq * (*Nqdata), s… in DMCeedCreateGeometryFVM()
[all …]
/petsc/src/ts/utils/dmplexlandau/kokkos/
H A Dlandau.kokkos.cxx137 PetscErrorCode LandauKokkosStaticDataSet(DM plex, const PetscInt Nq, const PetscInt Nb, const Petsc… in LandauKokkosStaticDataSet() argument
158 nip += a_numCells[grid] * Nq; in LandauKokkosStaticDataSet()
160 IPf_sz += Nq * nfloc * a_numCells[grid]; in LandauKokkosStaticDataSet()
178 … Kokkos::LayoutLeft, Kokkos::HostSpace, Kokkos::MemoryTraits<Kokkos::Unmanaged>> h_BB(BB, Nq * Nb); in LandauKokkosStaticDataSet()
179 … B = new Kokkos::View<PetscReal *, Kokkos::LayoutLeft>("B", Nq * Nb); in LandauKokkosStaticDataSet()
181 …s::LayoutLeft, Kokkos::HostSpace, Kokkos::MemoryTraits<Kokkos::Unmanaged>> h_DD(DD, Nq * Nb * dim); in LandauKokkosStaticDataSet()
182 … D = new Kokkos::View<PetscReal *, Kokkos::LayoutLeft>("D", Nq * Nb * dim); in LandauKokkosStaticDataSet()
381 PetscErrorCode LandauKokkosJacobian(DM plex[], const PetscInt Nq, const PetscInt Nb, const PetscInt… in LandauKokkosJacobian() argument
388 int nfaces = 0, vector_size = 256 / Nq; in LandauKokkosJacobian()
394 …efaultExecutionSpace().concurrency(), openmp = !!(conc < 1000), team_size = (openmp == 0) ? Nq : 1; in LandauKokkosJacobian()
[all …]
/petsc/src/dm/dt/dualspace/impls/simple/
H A Ddspacesimple.c73 PetscInt Nc, c, Nq, p; in PetscDualSpaceSimpleSetFunctional_Simple() local
79 …PetscCall(PetscQuadratureGetData(sp->functional[f], NULL, &Nc, &Nq, NULL, (const PetscReal **)&wei… in PetscDualSpaceSimpleSetFunctional_Simple()
83 for (p = 0; p < Nq; ++p) vol += weights[p * Nc + c]; in PetscDualSpaceSimpleSetFunctional_Simple()
84 for (p = 0; p < Nq; ++p) weights[p * Nc + c] /= (vol == 0.0 ? 1.0 : vol); in PetscDualSpaceSimpleSetFunctional_Simple()
/petsc/src/dm/impls/swarm/tests/
H A Dex2.c274 PetscInt cStart, cEnd, c, Nq, q, dim, Nfc; in CreateParticles_Shape() local
283 PetscCall(PetscQuadratureGetData(quad, NULL, NULL, &Nq, &qpoints, &qweights)); in CreateParticles_Shape()
295 PetscCall(DMSwarmSetLocalSizes(*sw, (cEnd - cStart) * Nq, 0)); in CreateParticles_Shape()
303 for (q = 0; q < Nq; ++q) { in CreateParticles_Shape()
305 swarm_cellid[c * Nq + q] = c; in CreateParticles_Shape()
306 CoordinatesRefToReal(dim, dim, xi0, v0, J, &qpoints[q * dim], &coords[(c * Nq + q) * dim]); in CreateParticles_Shape()
307 PetscCall(linear(dim, 0.0, &coords[(c * Nq + q) * dim], 1, &vals[c * Nq + q], (void *)user)); in CreateParticles_Shape()
308 vals[c * Nq + q] *= qweights[q] * detJ; in CreateParticles_Shape()
450 …PetscInt dim, coordDim, numFields, numComponents = 0, qNc, Nq, cStart, cEnd, vStart, vEnd,… in InterpolateGradient() local
480 PetscCall(PetscQuadratureGetData(quad, NULL, &qNc, &Nq, &quadPoints, &quadWeights)); in InterpolateGradient()
[all …]
/petsc/src/dm/dt/fe/impls/basic/
H A Dfebasic.c60 PetscInt Nc, Nq, q, k, c; in PetscFESetUp_Basic() local
63 PetscCall(PetscQuadratureGetData(f, NULL, &Nc, &Nq, &points, &weights)); in PetscFESetUp_Basic()
64 PetscCall(PetscMalloc1(Nc * Nq * pdim, &Bf)); in PetscFESetUp_Basic()
65 PetscCall(PetscSpaceEvaluate(fem->basisSpace, Nq, points, Bf, NULL, NULL)); in PetscFESetUp_Basic()
70 for (q = 0; q < Nq; ++q) { in PetscFESetUp_Basic()
176 PetscInt qNc, Nq, q; in PetscFEIntegrate_Basic() local
203 PetscCall(PetscQuadratureGetData(quad, NULL, &qNc, &Nq, &quadPoints, &quadWeights)); in PetscFEIntegrate_Basic()
221 for (q = 0; q < Nq; ++q) { in PetscFEIntegrate_Basic()
267 PetscInt qNc, Nq, q, Np, dE; in PetscFEIntegrateBd_Basic() local
297 PetscCall(PetscQuadratureGetData(quad, NULL, &qNc, &Nq, &quadPoints, &quadWeights)); in PetscFEIntegrateBd_Basic()
[all …]
/petsc/src/dm/impls/plex/tests/output/
H A Dex65_ext_coord_space.out3 Field Q1 FEM 1 component (implicit) (Nq 2 Nqc 1) 1-jet
19 Field Q1 FEM 2 components (implicit) (Nq 4 Nqc 1) 1-jet
H A Dex5_quad_0.out25 Field displacement FEM 2 components (implicit) (Nq 2 Nqc 1) 1-jet
50 Field fault traction FEM 2 components (implicit) (Nq 2 Nqc 1) 1-jet
H A Dex5_hex_0.out61 Field displacement FEM 3 components (implicit) (Nq 4 Nqc 1) 1-jet
86 Field fault traction FEM 3 components (implicit) (Nq 4 Nqc 1) 1-jet
H A Dex5_hex_t1_0.out113 Field displacement FEM 3 components (implicit) (Nq 4 Nqc 1) 1-jet
138 Field fault traction FEM 3 components (implicit) (Nq 4 Nqc 1) 1-jet
H A Dex5_hex_t2_0.out86 Field displacement FEM 3 components (implicit) (Nq 4 Nqc 1) 1-jet
111 Field fault traction FEM 3 components (implicit) (Nq 4 Nqc 1) 1-jet
H A Dex5_tri_0.out84 Field displacement FEM 2 components (implicit) (Nq 2 Nqc 1) 1-jet
105 Field fault traction FEM 2 components (implicit) (Nq 2 Nqc 1) 1-jet
H A Dex5_tri_0_perm.out84 Field displacement FEM 2 components (implicit) (Nq 2 Nqc 1) 1-jet
105 Field fault traction FEM 2 components (implicit) (Nq 2 Nqc 1) 1-jet
H A Dex5_tri_t1_0.out129 Field displacement FEM 2 components (implicit) (Nq 2 Nqc 1) 1-jet
150 Field fault traction FEM 2 components (implicit) (Nq 2 Nqc 1) 1-jet
H A Dex5_tri_t5_0.out120 Field displacement FEM 2 components (implicit) (Nq 2 Nqc 1) 1-jet
141 Field fault traction FEM 2 components (implicit) (Nq 2 Nqc 1) 1-jet
H A Dex5_tri_t1_0_perm.out129 Field displacement FEM 2 components (implicit) (Nq 2 Nqc 1) 1-jet
150 Field fault traction FEM 2 components (implicit) (Nq 2 Nqc 1) 1-jet
/petsc/src/dm/impls/plex/
H A Dplexfem.c1469 …PetscInt dim, coordDim, numFields, numComponents = 0, qNc, Nq, cellHeight, cStart, cEnd, c… in DMPlexComputeL2DiffLocal() local
1502 PetscCall(PetscQuadratureGetData(quad, NULL, &qNc, &Nq, NULL, &quadWeights)); in DMPlexComputeL2DiffLocal()
1504 …, &interpolant, coordDim * (Nq + 1), &coords, Nq, &fegeom.detJ, coordDim * coordDim * Nq, &fegeom.… in DMPlexComputeL2DiffLocal()
1535 for (q = 0; q < Nq; ++q) { in DMPlexComputeL2DiffLocal()
1545 gcoords = &coords[coordDim * Nq]; in DMPlexComputeL2DiffLocal()
1594 …PetscInt dim, coordDim, qNc = 0, Nq = 0, numFields, numComponents = 0, cStart, cEnd, c, fi… in DMComputeL2GradientDiff_Plex() local
1618 PetscCall(PetscQuadratureGetData(quad, NULL, &qNc, &Nq, NULL, &quadWeights)); in DMComputeL2GradientDiff_Plex()
1621 …dDim * (Nq + 1), &coords, coordDim * coordDim * Nq, &fegeom.J, coordDim * coordDim * Nq, &fegeom.i… in DMComputeL2GradientDiff_Plex()
1644 for (q = 0; q < Nq; ++q) { in DMComputeL2GradientDiff_Plex()
1654 gcoords = &coords[coordDim * Nq]; in DMComputeL2GradientDiff_Plex()
[all …]
/petsc/src/dm/dt/fe/impls/composite/
H A Dfecomposite.c75 PetscInt Nc, Nq, q, k; in PetscFESetUp_Composite() local
78 PetscCall(PetscQuadratureGetData(f, NULL, &Nc, &Nq, &points, &weights)); in PetscFESetUp_Composite()
80 PetscCall(PetscSpaceEvaluate(fem->basisSpace, Nq, points, Bf, NULL, NULL)); in PetscFESetUp_Composite()
84 …for (q = 0; q < Nq; ++q) invVscalar[(s * spdim + j) * spdim + k] += Bf[q * spdim + k] * weights[q]; in PetscFESetUp_Composite()
/petsc/src/ts/tutorials/hamiltonian/output/
H A Dex3_sphere_3d.out6 0) FormLandau 1512 IPs, 56 cells total, Nb=27, Nq=27, dim=3, Tab: Nb=27 Nf=1 Np=27 cdim=3 N=517 nnz…
/petsc/src/ts/tests/output/
H A Dex30_cgns.out6 0) FormLandau 128 IPs, 8 cells total, Nb=16, Nq=16, dim=2, Tab: Nb=16 Nf=1 Np=16 cdim=2 N=364 nnz=7…
/petsc/src/dm/dt/dualspace/interface/
H A Ddualspace.c1144 PetscInt dim, dE, qNc, c, Nq, q; in PetscDualSpaceApplyDefault() local
1152 PetscCall(PetscQuadratureGetData(n, &dim, &qNc, &Nq, &points, &weights)); in PetscDualSpaceApplyDefault()
1159 for (q = 0; q < Nq; ++q) { in PetscDualSpaceApplyDefault()
1570 PetscInt dimEmbed, qNc, c, Nq, q; in PetscDualSpaceApplyFVM() local
1578 PetscCall(PetscQuadratureGetData(n, NULL, &qNc, &Nq, &points, &weights)); in PetscDualSpaceApplyFVM()
1582 for (q = 0; q < Nq; ++q) { in PetscDualSpaceApplyFVM()
2142 PetscErrorCode PetscDualSpacePullback(PetscDualSpace dsp, PetscFEGeom *fegeom, PetscInt Nq, PetscIn… in PetscDualSpacePullback() argument
2168 PetscCall(PetscDualSpaceTransform(dsp, trans, PETSC_TRUE, fegeom, Nq, Nc, pointEval)); in PetscDualSpacePullback()
2194 PetscErrorCode PetscDualSpacePushforward(PetscDualSpace dsp, PetscFEGeom *fegeom, PetscInt Nq, Pets… in PetscDualSpacePushforward() argument
2220 PetscCall(PetscDualSpaceTransform(dsp, trans, PETSC_FALSE, fegeom, Nq, Nc, pointEval)); in PetscDualSpacePushforward()
[all …]
/petsc/src/ts/utils/dmplexlandau/
H A Dplexland.c89 PetscInt numCells[LANDAU_MAX_GRIDS], Nq, Nb; in LandauFormJacobian_Internal() local
122 PetscCall(PetscQuadratureGetData(quad, NULL, NULL, &Nq, NULL, NULL)); in LandauFormJacobian_Internal()
124 …PetscCheck(Nq <= LANDAU_MAX_NQND, ctx->comm, PETSC_ERR_ARG_WRONG, "Order too high. Nq = %" PetscIn… in LandauFormJacobian_Internal()
204 …PetscCall(LandauKokkosJacobian(ctx->plex, Nq, Nb, ctx->batch_sz, ctx->num_grids, numCells, Eq_m, c… in LandauFormJacobian_Internal()
228 ip_offset[grid + 1] = ip_offset[grid] + numCells[grid] * Nq; in LandauFormJacobian_Internal()
229 ipf_offset[grid + 1] = ipf_offset[grid] + Nq * nfloc * numCells[grid]; in LandauFormJacobian_Internal()
249 …const PetscInt loc_nip = numCells[grid] * Nq, loc_Nf = ctx->species_offset[grid + 1] - ctx->specie… in LandauFormJacobian_Internal()
252 …PetscReal *invJe = &invJ_a[(ip_offset[grid] + loc_elem * Nq) * dim * dim]; // ingJ is static d… in LandauFormJacobian_Internal()
277 for (PetscInt qi = 0; qi < Nq; qi++) { in LandauFormJacobian_Internal()
284 … const PetscInt idx = b_id * IPf_sz_glb + ipf_offset[grid] + f * loc_nip + loc_elem * Nq + qi; in LandauFormJacobian_Internal()
[all …]

123