Home
last modified time | relevance | path

Searched refs:Nb (Results 1 – 25 of 46) sorted by relevance

12

/petsc/include/petsc/private/
H A Dpetscfeimpl.h304 const PetscInt Nb = T->Nb; in PetscFEInterpolate_Static() local
308 for (PetscInt f = 0; f < Nb; ++f) interpolant[fc] += x[f] * basis[(q * Nb + f) * Nc + fc]; in PetscFEInterpolate_Static()
322 const PetscInt Nb = T->Nb; in PetscFEInterpolateAtPoints_Static() local
326 for (f = 0; f < Nb; ++f) interpolant[fc] += x[f] * basis[(q * Nb + f) * Nc + fc]; in PetscFEInterpolateAtPoints_Static()
343 const PetscInt Nb = T->Nb; in PetscFEInterpolateGradient_Static() local
350 for (f = 0; f < Nb; ++f) { in PetscFEInterpolateGradient_Static()
351 …for (d = 0; d < cdim; ++d) interpolant[fc * cdim + d] += x[f] * basisDer[((q * Nb + f) * Nc + fc) … in PetscFEInterpolateGradient_Static()
359 for (f = 0; f < Nb; ++f) { in PetscFEInterpolateGradient_Static()
360 …dim * cdim; ++d) interpolant[off + fc * cdim + d] += x[f] * basisHes[((q * Nb + f) * Nc + fc) * cd… in PetscFEInterpolateGradient_Static()
373 PetscInt Nb, Nc, fc, f, d, g; in PetscFEFreeInterpolateGradient_Static() local
[all …]
/petsc/src/dm/dt/space/impls/poly/
H A Dspacepoly.c140 PetscInt Nb, jet, Njet; in PetscSpaceEvaluate_Polynomial() local
150 PetscCall(PetscDTBinomialInt(dim + sp->degree, dim, &Nb)); in PetscSpaceEvaluate_Polynomial()
159 PetscCall(DMGetWorkArray(dm, Nb * Njet * npoints, MPIU_REAL, &pScalar)); in PetscSpaceEvaluate_Polynomial()
169 PetscInt p_strl = Nb; in PetscSpaceEvaluate_Polynomial()
175 PetscCall(PetscArrayzero(B, npoints * Nb)); in PetscSpaceEvaluate_Polynomial()
176 for (PetscInt b = 0; b < Nb; b++) { in PetscSpaceEvaluate_Polynomial()
181 PetscInt p_strl = dim * Nb; in PetscSpaceEvaluate_Polynomial()
189 PetscCall(PetscArrayzero(D, npoints * Nb * dim)); in PetscSpaceEvaluate_Polynomial()
191 for (PetscInt b = 0; b < Nb; b++) { in PetscSpaceEvaluate_Polynomial()
197 PetscInt p_strl = dim * dim * Nb; in PetscSpaceEvaluate_Polynomial()
[all …]
/petsc/src/dm/dt/space/impls/ptrimmed/tests/
H A Dex1.c9 PetscInt Nf, Nb; in test() local
31 PetscCall(PetscSpaceGetDimension(sp, &Nb)); in test()
32 …PetscCheck(Nb == Nbexp, comm, PETSC_ERR_PLIB, "Space dimension mismatch, %" PetscInt_FMT " != %" P… in test()
42 Bsize = npoints * Nb * Nf * nCopies; in test()
/petsc/src/dm/dt/tests/
H A Dex15.c14 PetscInt Nb, Nb_min; in testQuadrature() local
21 PetscCall(PetscDTBinomialInt(dim + p_degree, dim, &Nb)); in testQuadrature()
23 PetscCall(PetscMalloc1(num_points * Nb, &eval)); in testQuadrature()
26 for (PetscInt j = i; j < Nb; j++) { in testQuadrature()
H A Dex16.c14 PetscInt Nb, vNb, vNb_s, vNb_d, Nc, vNc, cdim; in PetscFEVectorTest() local
30 PetscCall(PetscFEGetDimension(orig_fe, &Nb)); in PetscFEVectorTest()
43 for (PetscInt b = 0; b < Nb; b++) { in PetscFEVectorTest()
53 const PetscReal *s_Tp = &s_Tk[(p * Nb + b) * Nc * dblock]; in PetscFEVectorTest()
/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
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()
343 … const PetscScalar Aij, const PetscInt f, const PetscInt g, const PetscInt Nb, PetscInt moffset, c… in landau_mat_assemble() argument
381 PetscErrorCode LandauKokkosJacobian(DM plex[], const PetscInt Nq, const PetscInt Nb, const PetscInt… in LandauKokkosJacobian() argument
480 …const int totDim_max = Nf_max * Nb, elem_mat_size_max = … in LandauKokkosJacobian()
497 closure_sz += Nb * nfloc * a_numCells[grid]; in LandauKokkosJacobian()
511 const int scr_bytes_fdf = real2_scr_t::shmem_size(Nf_max, Nb); in LandauKokkosJacobian()
[all …]
/petsc/src/dm/dt/fe/tests/
H A Dex2.c177 PetscInt dim, Nb, b, Nc, c, Np, p; in TestEvaluation() local
205 PetscCall(PetscSpaceGetDimension(sp, &Nb)); in TestEvaluation()
207 PetscCall(DMGetWorkArray(dm, Np * Nb * Nc, MPIU_REAL, &B)); in TestEvaluation()
208 PetscCall(DMGetWorkArray(dm, Np * Nb * Nc * dim, MPIU_REAL, &D)); in TestEvaluation()
209 PetscCall(DMGetWorkArray(dm, Np * Nb * Nc * dim * dim, MPIU_REAL, &H)); in TestEvaluation()
213 for (b = 0; b < Nb; ++b) { in TestEvaluation()
215 …for (c = 0; c < Nc; ++c) PetscCall(PetscPrintf(PETSC_COMM_SELF, " %g", (double)B[(p * Nb + b) * Nc… in TestEvaluation()
220 …for (d = 0; d < dim; ++d) PetscCall(PetscPrintf(PETSC_COMM_SELF, " %g", (double) B[((p*Nb+b)*Nc+c)… in TestEvaluation()
226 PetscCall(DMRestoreWorkArray(dm, Np * Nb, MPIU_REAL, &B)); in TestEvaluation()
227 PetscCall(DMRestoreWorkArray(dm, Np * Nb * dim, MPIU_REAL, &D)); in TestEvaluation()
[all …]
/petsc/src/dm/dt/space/impls/ptrimmed/
H A Dspaceptrimmed.c143 PetscInt Nb; in PetscSpaceEvaluate_Ptrimmed() local
164 PetscCall(PetscDTPTrimmedSize(dim, degree, f, &Nb)); in PetscSpaceEvaluate_Ptrimmed()
165 PetscCall(DMGetWorkArray(dm, Nb * Nf * Njet * npoints, MPIU_REAL, &eval)); in PetscSpaceEvaluate_Ptrimmed()
168 PetscInt p_strl = Nf * Nb; in PetscSpaceEvaluate_Ptrimmed()
177 for (PetscInt b = 0; b < Nb; b++) { in PetscSpaceEvaluate_Ptrimmed()
183 PetscInt p_strl = dim * Nf * Nb; in PetscSpaceEvaluate_Ptrimmed()
195 for (PetscInt b = 0; b < Nb; b++) { in PetscSpaceEvaluate_Ptrimmed()
202 PetscInt p_strl = dim * dim * Nf * Nb; in PetscSpaceEvaluate_Ptrimmed()
224 for (PetscInt b = 0; b < Nb; b++) { in PetscSpaceEvaluate_Ptrimmed()
232 PetscCall(DMRestoreWorkArray(dm, Nb * Nf * Njet * npoints, MPIU_REAL, &eval)); in PetscSpaceEvaluate_Ptrimmed()
/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…
H A Dex30_conserve.out6 0) FormLandau 144 IPs, 12 cells total, Nb=10, Nq=12, dim=2, Tab: Nb=10 Nf=1 Np=10 cdim=2 N=268 nnz=…
H A Dex30_3d.out6 0) FormLandau 432 IPs, 16 cells total, Nb=27, Nq=27, dim=3, Tab: Nb=27 Nf=2 Np=27 cdim=3 N=250 nnz=…
H A Dex30_0.out16 0) FormLandau 384 IPs, 24 cells total, Nb=16, Nq=16, dim=2, Tab: Nb=16 Nf=3 Np=16 cdim=2 N=1092 nnz…
/petsc/src/dm/dt/fe/interface/
H A Dfe.c982 PetscInt Nb; /* Dimension of FE space P */ in PetscFECreateTabulation() local
998 PetscCall(PetscDualSpaceGetDimension(Q, &Nb)); in PetscFECreateTabulation()
1004 (*T)->Nb = Nb; in PetscFECreateTabulation()
1008 …for (k = 0; k <= (*T)->K; ++k) PetscCall(PetscCalloc1(nrepl * npoints * Nb * Nc * PetscPowInt(cdim… in PetscFECreateTabulation()
1049 PetscInt Nb; /* Dimension of FE space P */ in PetscFEComputeTabulation() local
1056 PetscCall(PetscDualSpaceGetDimension(Q, &Nb)); in PetscFEComputeTabulation()
1059 …tscCheck(T->Nb == Nb, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Tabulation Nb %" PetscInt_FMT " must ma… in PetscFEComputeTabulation()
2369 const PetscInt Nbf = T[f]->Nb; in PetscFEEvaluateFieldJets_Internal()
2436 const PetscInt Nbf = T->Nb; in PetscFEEvaluateFieldJets_Hybrid_Internal()
2490 const PetscInt Nb = Tc->Nb; in PetscFEEvaluateFaceFields_Internal() local
[all …]
/petsc/src/dm/dt/interface/ftn-custom/
H A Dzdtdsf90.c20 PetscInt Nb; member
41 size = tab[*i - 1]->Nr * tab[*i - 1]->Np * tab[*i - 1]->Nb * tab[*i - 1]->Nc; in petscdsgettabulationsetpointers_()
/petsc/src/ts/utils/dmplexlandau/tutorials/output/
H A Dex2_single.out6 0) FormLandau 512 IPs, 32 cells total, Nb=16, Nq=16, dim=2, Tab: Nb=16 Nf=1 Np=16 cdim=2 N=325 nnz=…
H A Dex2_sphere_3d.out6 0) FormLandau 448 IPs, 7 cells total, Nb=64, Nq=64, dim=3, Tab: Nb=64 Nf=1 Np=64 cdim=3 N=232 nnz= …
H A Dex1_sphere.out6 0) FormLandau 216 IPs, 24 cells total, Nb=9, Nq=9, dim=2, Tab: Nb=9 Nf=2 Np=9 cdim=2 N=110 nnz= 1522
H A Dex1_0.out6 0) FormLandau 736 IPs, 46 cells total, Nb=16, Nq=16, dim=2, Tab: Nb=16 Nf=2 Np=16 cdim=2 N=443 nnz=…
H A Dex1f90_0.out16 0) FormLandau 352 IPs, 22 cells total, Nb=16, Nq=16, dim=2, Tab: Nb=16 Nf=3 Np=16 cdim=2 N=324 nnz=…
H A Dex2_simplex.out6 0) FormLandau 576 IPs, 96 cells total, Nb=6, Nq=6, dim=2, Tab: Nb=6 Nf=2 Np=6 cdim=2 N=226 nnz=2338
/petsc/src/ts/utils/dmplexlandau/
H A Dplexland.c89 PetscInt numCells[LANDAU_MAX_GRIDS], Nq, Nb; in LandauFormJacobian_Internal() local
123 PetscCall(PetscFEGetDimension(ctx->fe[0], &Nb)); in LandauFormJacobian_Internal()
125 …PetscCheck(Nb <= LANDAU_MAX_NQND, ctx->comm, PETSC_ERR_ARG_WRONG, "Order too high. Nb = %" PetscIn… in LandauFormJacobian_Internal()
155 …for (PetscInt grid = 0; grid < ctx->num_grids; grid++) cellClosure_sz += Nb * Nf[grid] * numCells[… in LandauFormJacobian_Internal()
173 …PetscCall(PetscMemcpy(cellClosure_it, coef, Nb * Nf[grid] * sizeof(*cellClosure_it))); /* change i… in LandauFormJacobian_Internal()
175 cellClosure_it += Nb * Nf[grid]; in LandauFormJacobian_Internal()
204 …PetscCall(LandauKokkosJacobian(ctx->plex, Nq, Nb, ctx->batch_sz, ctx->num_grids, numCells, Eq_m, c… in LandauFormJacobian_Internal()
255 …coef = &cellClosure[b_id * IPf_sz_glb + ipf_offset[grid] + loc_elem * Nb * loc_Nf]; // this is con… in LandauFormJacobian_Internal()
260 for (b = 0; b < Nb; ++b) { in LandauFormJacobian_Internal()
263 coef[f * Nb + b] = xdata[idx + moffset]; in LandauFormJacobian_Internal()
[all …]
/petsc/src/dm/impls/plex/generators/tetgen/
H A Dtetgenerate.cxx231 int Nb, oclass, mtype, *senses; 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()
274 for (b = 0; b < Nb; ++b) { in DMPlexGenerate_Tetgen()
281 if (b < Nb) { in DMPlexGenerate_Tetgen()
510 int Nb, oclass, mtype, *senses; 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()
553 for (b = 0; b < Nb; ++b) { in DMPlexRefine_Tetgen()
560 if (b < Nb) { in DMPlexRefine_Tetgen()
/petsc/src/dm/impls/plex/generators/ctetgen/
H A Dctetgenerate.c215 int Nb, oclass, mtype, *senses; in DMPlexGenerate_CTetgen() local
221 PetscCall(EG_getTopology(model, &geom, &oclass, &mtype, NULL, &Nb, &bodies, &senses)); in DMPlexGenerate_CTetgen()
249 for (b = 0; b < Nb; ++b) { in DMPlexGenerate_CTetgen()
252 if (b < Nb) { in DMPlexGenerate_CTetgen()
490 int Nb, oclass, mtype, *senses; in DMPlexRefine_CTetgen() local
496 PetscCall(EG_getTopology(model, &geom, &oclass, &mtype, NULL, &Nb, &bodies, &senses)); in DMPlexRefine_CTetgen()
524 for (b = 0; b < Nb; ++b) { in DMPlexRefine_CTetgen()
527 if (b < Nb) { in DMPlexRefine_CTetgen()
/petsc/src/dm/impls/plex/
H A Dplexfem.c1519 PetscInt Nb, Nc, q, fc; in DMPlexComputeL2DiffLocal() local
1525 PetscCall(PetscFEGetDimension((PetscFE)obj, &Nb)); in DMPlexComputeL2DiffLocal()
1528 Nb = 1; in DMPlexComputeL2DiffLocal()
1533 PetscCall(DMPrintCellVector(c, title, Nb, &x[fieldOffset])); in DMPlexComputeL2DiffLocal()
1569 fieldOffset += Nb; in DMPlexComputeL2DiffLocal()
1634 PetscInt Nb, Nc, q, fc; in DMComputeL2GradientDiff_Plex() local
1637 PetscCall(PetscFEGetDimension(fe, &Nb)); in DMComputeL2GradientDiff_Plex()
1642 PetscCall(DMPrintCellVector(c, title, Nb, &x[fieldOffset])); in DMComputeL2GradientDiff_Plex()
1683 fieldOffset += Nb; in DMComputeL2GradientDiff_Plex()
1774 PetscInt Nb, Nc, q, fc; in DMComputeL2FieldDiff_Plex() local
[all …]

12