2650   const PetscInt   NbI       = TI->Nb;  in PetscFEUpdateElementMat_Internal()  local
2652   const PetscReal *basisI    = &TI->T[0][(r * NqI + q) * NbI * NcI];  in PetscFEUpdateElementMat_Internal()
2653   const PetscReal *basisDerI = &TI->T[1][(r * NqI + q) * NbI * NcI * cdim];  in PetscFEUpdateElementMat_Internal()
2660   for (PetscInt f = 0; f < NbI; ++f) {  in PetscFEUpdateElementMat_Internal()
2668   PetscCall(PetscFEPushforward(feI, fegeom, NbI, tmpBasisI));  in PetscFEUpdateElementMat_Internal()
2669   PetscCall(PetscFEPushforwardGradient(feI, fegeom, NbI, tmpBasisDerI));  in PetscFEUpdateElementMat_Internal()
2686     for (PetscInt f = 0; f < NbI; ++f) {  in PetscFEUpdateElementMat_Internal()
2704       petsc_elemmat_kernel_g1(NbI, NcI, NbJ, NcJ, 2);  in PetscFEUpdateElementMat_Internal()
2706       petsc_elemmat_kernel_g1(NbI, NcI, NbJ, NcJ, 3);  in PetscFEUpdateElementMat_Internal()
2708       petsc_elemmat_kernel_g1(NbI, NcI, NbJ, NcJ, dE);  in PetscFEUpdateElementMat_Internal()
2711     for (PetscInt f = 0; f < NbI; ++f) {  in PetscFEUpdateElementMat_Internal()
2734       petsc_elemmat_kernel_g2(NbI, NcI, NbJ, NcJ, 2);  in PetscFEUpdateElementMat_Internal()
2736       petsc_elemmat_kernel_g2(NbI, NcI, NbJ, NcJ, 3);  in PetscFEUpdateElementMat_Internal()
2738       petsc_elemmat_kernel_g2(NbI, NcI, NbJ, NcJ, dE);  in PetscFEUpdateElementMat_Internal()
2747         for (PetscInt f = 0; f < NbI; ++f) {  in PetscFEUpdateElementMat_Internal()
2764       petsc_elemmat_kernel_g3(NbI, NcI, NbJ, NcJ, 2);  in PetscFEUpdateElementMat_Internal()
2766       petsc_elemmat_kernel_g3(NbI, NcI, NbJ, NcJ, 3);  in PetscFEUpdateElementMat_Internal()
2768       petsc_elemmat_kernel_g3(NbI, NcI, NbJ, NcJ, dE);  in PetscFEUpdateElementMat_Internal()
2771     for (PetscInt f = 0; f < NbI; ++f) {  in PetscFEUpdateElementMat_Internal()
2804   const PetscInt   NbI       = TI->Nb;  in PetscFEUpdateElementMat_Hybrid_Internal()  local
2806   const PetscReal *basisI    = &TI->T[0][(r * NqI + q) * NbI * NcI];  in PetscFEUpdateElementMat_Hybrid_Internal()
2807   const PetscReal *basisDerI = &TI->T[1][(r * NqI + q) * NbI * NcI * dE];  in PetscFEUpdateElementMat_Hybrid_Internal()
2817   for (f = 0; f < NbI; ++f) {  in PetscFEUpdateElementMat_Hybrid_Internal()
2825   PetscCall(PetscFEPushforward(feI, fegeom, NbI, tmpBasisI));  in PetscFEUpdateElementMat_Hybrid_Internal()
2826   PetscCall(PetscFEPushforwardGradient(feI, fegeom, NbI, tmpBasisDerI));  in PetscFEUpdateElementMat_Hybrid_Internal()
2838   for (f = 0; f < NbI; ++f) {  in PetscFEUpdateElementMat_Hybrid_Internal()
2841       const PetscInt i    = offsetI + NbI * so + f; /* Element matrix row */  in PetscFEUpdateElementMat_Hybrid_Internal()
2902   const PetscInt NbI     = TI->Nb;
2917       if (g0[i]) g0_i += NbI * NbJ;
2920       if (g1[i]) g1_i += NbI * NbJ;
2923       if (g2[i]) g2_i += NbI * NbJ;
2926       if (g3[i]) g3_i += NbI * NbJ;
2927   if (g0_i == NbI * NbJ * NcI * NcJ) g0_i = 0;
2928   if (g1_i == NbI * NbJ * NcI * NcJ * dE) g1_i = 0;
2929   if (g2_i == NbI * NbJ * NcI * NcJ * dE) g2_i = 0;
2930   if (g3_i == NbI * NbJ * NcI * NcJ * dE * dE) g3_i = 0;
2941   for (PetscInt f = 0; f < NbI; ++f) {