2655   const PetscInt   NbJ       = TJ->Nb;  in PetscFEUpdateElementMat_Internal()  local
2657   const PetscReal *basisJ    = &TJ->T[0][(r * NqJ + q) * NbJ * NcJ];  in PetscFEUpdateElementMat_Internal()
2658   const PetscReal *basisDerJ = &TJ->T[1][(r * NqJ + q) * NbJ * NcJ * cdim];  in PetscFEUpdateElementMat_Internal()
2671     for (PetscInt g = 0; g < NbJ; ++g) {  in PetscFEUpdateElementMat_Internal()
2679     PetscCall(PetscFEPushforward(feJ, fegeom, NbJ, tmpBasisJ));  in PetscFEUpdateElementMat_Internal()
2680     PetscCall(PetscFEPushforwardGradient(feJ, fegeom, NbJ, tmpBasisDerJ));  in PetscFEUpdateElementMat_Internal()
2692         for (PetscInt g = 0; g < NbJ; ++g) {  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()
2717         for (PetscInt g = 0; g < NbJ; ++g) {  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()
2741     for (PetscInt g = 0; g < NbJ; ++g) {  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()
2777         for (PetscInt g = 0; g < NbJ; ++g) {  in PetscFEUpdateElementMat_Internal()
2809   const PetscInt   NbJ       = TJ->Nb;  in PetscFEUpdateElementMat_Hybrid_Internal()  local
2811   const PetscReal *basisJ    = &TJ->T[0][(r * NqJ + q) * NbJ * NcJ];  in PetscFEUpdateElementMat_Hybrid_Internal()
2812   const PetscReal *basisDerJ = &TJ->T[1][(r * NqJ + q) * NbJ * NcJ * dE];  in PetscFEUpdateElementMat_Hybrid_Internal()
2827   for (g = 0; g < NbJ; ++g) {  in PetscFEUpdateElementMat_Hybrid_Internal()
2835   PetscCall(PetscFEPushforward(feJ, fegeom, NbJ, tmpBasisJ));  in PetscFEUpdateElementMat_Hybrid_Internal()
2842       for (g = 0; g < NbJ; ++g) {  in PetscFEUpdateElementMat_Hybrid_Internal()
2845           const PetscInt j    = offsetJ + NbJ * to + g; /* Element matrix column */  in PetscFEUpdateElementMat_Hybrid_Internal()
2904   const PetscInt NbJ     = TJ->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;
2946       for (PetscInt g = 0; g < NbJ; ++g) {