982   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
2494     for (b = 0; b < Nb; ++b) {  in PetscFEEvaluateFaceFields_Internal()
2495       for (c = 0; c < Nc; ++c) u[c] += coefficients[b] * faceBasis[(faceLoc * Nb + b) * Nc + c];  in PetscFEEvaluateFaceFields_Internal()
2507   const PetscInt   Nb       = T->Nb;  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()
2514     for (b = 0; b < Nb; ++b) {  in PetscFEUpdateElementVec_Internal()
2518         tmpBasis[bcidx] = basis[q * Nb * Nc + bcidx];  in PetscFEUpdateElementVec_Internal()
2519 …for (d = 0; d < dEt; ++d) tmpBasisDer[bcidx * dE + d] = basisDer[q * Nb * Nc * dEt + bcidx * dEt +…  in PetscFEUpdateElementVec_Internal()
2524     PetscCall(PetscFEPushforward(fe, &pgeom, Nb, tmpBasis));  in PetscFEUpdateElementVec_Internal()
2525     PetscCall(PetscFEPushforwardGradient(fe, &pgeom, Nb, tmpBasisDer));  in PetscFEUpdateElementVec_Internal()
2526     for (b = 0; b < Nb; ++b) {  in PetscFEUpdateElementVec_Internal()
2543   const PetscInt   Nb       = T->Nb;  in PetscFEUpdateElementVec_Hybrid_Internal()  local
2545   const PetscReal *basis    = &T->T[0][r * Nq * Nb * Nc];  in PetscFEUpdateElementVec_Hybrid_Internal()
2546   const PetscReal *basisDer = &T->T[1][r * Nq * Nb * Nc * dE];  in PetscFEUpdateElementVec_Hybrid_Internal()
2549     for (PetscInt b = 0; b < Nb; ++b) {  in PetscFEUpdateElementVec_Hybrid_Internal()
2553         tmpBasis[bcidx] = basis[q * Nb * Nc + bcidx];  in PetscFEUpdateElementVec_Hybrid_Internal()
2554 …for (PetscInt d = 0; d < dE; ++d) tmpBasisDer[bcidx * dE + d] = basisDer[q * Nb * Nc * dE + bcidx …  in PetscFEUpdateElementVec_Hybrid_Internal()
2557     PetscCall(PetscFEPushforward(fe, fegeom, Nb, tmpBasis));  in PetscFEUpdateElementVec_Hybrid_Internal()
2563         for (PetscInt b = 0; b < Nb; ++b) {  in PetscFEUpdateElementVec_Hybrid_Internal()
2568             elemVec[Nb * s + b] += tmpBasis[bcidx] * f0[qcidx];  in PetscFEUpdateElementVec_Hybrid_Internal()
2569 …for (PetscInt d = 0; d < dE; ++d) elemVec[Nb * s + b] += tmpBasisDer[bcidx * dE + d] * f1[qcidx * …  in PetscFEUpdateElementVec_Hybrid_Internal()
2575       for (PetscInt b = 0; b < Nb; ++b) {  in PetscFEUpdateElementVec_Hybrid_Internal()
2580           elemVec[Nb * side + b] += tmpBasis[bcidx] * f0[qcidx];  in PetscFEUpdateElementVec_Hybrid_Internal()
2581 …for (PetscInt d = 0; d < dE; ++d) elemVec[Nb * side + b] += tmpBasisDer[bcidx * dE + d] * f1[qcidx…  in PetscFEUpdateElementVec_Hybrid_Internal()
2650   const PetscInt   NbI       = TI->Nb;  in PetscFEUpdateElementMat_Internal()
2655   const PetscInt   NbJ       = TJ->Nb;  in PetscFEUpdateElementMat_Internal()
2804   const PetscInt   NbI       = TI->Nb;  in PetscFEUpdateElementMat_Hybrid_Internal()
2809   const PetscInt   NbJ       = TJ->Nb;  in PetscFEUpdateElementMat_Hybrid_Internal()
2902   const PetscInt NbI     = TI->Nb;
2904   const PetscInt NbJ     = TJ->Nb;