Lines Matching refs:Nq
89 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()
324 …, grid, ctx->batch_sz, ctx->num_grids, ctx->mat_offset), totDim = loc_Nf * Nq, elemMatSize = totDi… in LandauFormJacobian_Internal()
326 const PetscReal *invJe = &invJ_a[(ip_offset[grid] + loc_elem * Nq) * dim * dim]; in LandauFormJacobian_Internal()
334 for (PetscInt qj = 0; qj < Nq; ++qj) { in LandauFormJacobian_Internal()
335 const PetscInt jpidx_glb = ip_offset[grid] + qj + loc_elem * Nq; in LandauFormJacobian_Internal()
350 PetscInt nip_loc_r = numCells[grid_r] * Nq, Nfloc_r = Nf[grid_r]; in LandauFormJacobian_Internal()
352 for (PetscInt qi = 0; qi < Nq; qi++, ipidx++) { in LandauFormJacobian_Internal()
367 … const PetscInt idx = b_id * IPf_sz_glb + ipf_offset[grid_r] + f * nip_loc_r + ei_r * Nq + qi; in LandauFormJacobian_Internal()
913 PetscInt Nq, Nb, *Nb2, cStart, cEnd, c, dim, qj, k; in adaptToleranceFEM() local
927 PetscCall(PetscQuadratureGetData(quad, NULL, NULL, &Nq, NULL, NULL)); in adaptToleranceFEM()
928 …PetscCheck(Nq <= LANDAU_MAX_NQND, ctx->comm, PETSC_ERR_ARG_WRONG, "Order too high. Nq = %" PetscIn… in adaptToleranceFEM()
944 for (qj = 0; qj < Nq; ++qj) { in adaptToleranceFEM()
951 …ell %" PetscInt_FMT ", qp %" PetscInt_FMT "/%" PetscInt_FMT "\n", grid, (double)r, c, qj + 1, Nq)); in adaptToleranceFEM()
957 … %" PetscInt_FMT "/%" PetscInt_FMT ", d=%e\n", grid, (double)r, c, qj + 1, Nq, (double)(r - minRad… in adaptToleranceFEM()
1007 …" cells, %" PetscInt_FMT " total quadrature points\n", grid, cEnd - cStart, Nq * (cEnd - cStart))); in adaptToleranceFEM()
1409 …PetscInt numCells[LANDAU_MAX_GRIDS], Nq, Nb, Nf[LANDAU_MAX_GRIDS], ncellsTot = 0, MAP_BF_S… in CreateStaticData() local
1462 PetscCall(PetscQuadratureGetData(quad, NULL, NULL, &Nq, NULL, &quadWeights)); in CreateStaticData()
1463 …PetscCheck(Nq <= LANDAU_MAX_NQND, ctx->comm, PETSC_ERR_ARG_WRONG, "Order too high. Nq = %" PetscIn… in CreateStaticData()
1509 maps[grid].num_face = (PetscInt)(pow(Nq, 1. / ((double)dim)) + .001); // Q in CreateStaticData()
1704 for (PetscInt grid = 0; grid < ctx->num_grids; grid++) nip_glb += Nq * numCells[grid]; in CreateStaticData()
1712 …" N=%" PetscInt_FMT " nnz= %" PetscInt_FMT "\n", 0, "FormLandau", nip_glb, ncells, Nb, Nq, dim, Nb, in CreateStaticData()
1756 invJe = invJ_a + outer_ej * Nq * dim * dim; in CreateStaticData()
1760 for (PetscInt qj = 0; qj < Nq; qj++, outer_ipidx++) { in CreateStaticData()
1816 …PetscCall(LandauKokkosStaticDataSet(ctx->plex[0], Nq, Nb, ctx->batch_sz, ctx->num_grids, numCells,… in CreateStaticData()