2592       for (PetscInt gc = 0; gc < (_NcJ); ++gc) { \
2593         const PetscScalar *G = g1 + (fc * (_NcJ) + gc) * _dE; \
2597             const PetscScalar *tBDJ = tmpBasisDerJ + (g * (_NcJ) + gc) * (_dE); \
2609     for (PetscInt gc = 0; gc < (_NcJ); ++gc) { \
2611         const PetscScalar *G = g2 + (fc * (_NcJ) + gc) * _dE; \
2613           const PetscScalar tBJv = tmpBasisJ[g * (_NcJ) + gc]; \
2628       for (PetscInt gc = 0; gc < (_NcJ); ++gc) { \
2629         const PetscScalar *G = g3 + (fc * (_NcJ) + gc) * (_dE) * (_dE); \
2634             const PetscScalar *tBDJ = tmpBasisDerJ + (g * (_NcJ) + gc) * (_dE); \
2672       for (PetscInt gc = 0; gc < NcJ; ++gc) {  in PetscFEUpdateElementMat_Internal()  local
2673         const PetscInt gidx = g * NcJ + gc; /* Trial function basis index */  in PetscFEUpdateElementMat_Internal()
2696 …for (PetscInt gc = 0; gc < NcJ; ++gc) elemMat[fOff] += bI * g0[fc * NcJ + gc] * tmpBasisJ[g * NcJ …  in PetscFEUpdateElementMat_Internal()  local
2721           for (PetscInt gc = 0; gc < NcJ; ++gc) {  in PetscFEUpdateElementMat_Internal()
2722             const PetscInt gidx = g * NcJ + gc; /* Trial function basis index */  in PetscFEUpdateElementMat_Internal()
2724 …for (PetscInt df = 0; df < dE; ++df) elemMat[fOff] += bI * g1[(fc * NcJ + gc) * dE + df] * tmpBasi…  in PetscFEUpdateElementMat_Internal()
2744       for (PetscInt gc = 0; gc < NcJ; ++gc) {  in PetscFEUpdateElementMat_Internal()
2745         const PetscScalar bJ = tmpBasisJ[g * NcJ + gc]; /* Trial function basis value */  in PetscFEUpdateElementMat_Internal()
2754 … df < dE; ++df) elemMat[fOff] += tmpBasisDerI[fidx * dE + df] * g2[(fc * NcJ + gc) * dE + df] * bJ;  in PetscFEUpdateElementMat_Internal()
2781           for (PetscInt gc = 0; gc < NcJ; ++gc) {  in PetscFEUpdateElementMat_Internal()
2782             const PetscInt gidx = g * NcJ + gc; /* Trial function basis index */  in PetscFEUpdateElementMat_Internal()
2785 … dE; ++dg) elemMat[fOff] += tmpBasisDerI[fidx * dE + df] * g3[((fc * NcJ + gc) * dE + df) * dE + d…  in PetscFEUpdateElementMat_Internal()
2815   PetscInt         f, fc, g, gc, df, dg;  in PetscFEUpdateElementMat_Hybrid_Internal()  local
2828     for (gc = 0; gc < NcJ; ++gc) {  in PetscFEUpdateElementMat_Hybrid_Internal()
2829       const PetscInt gidx = g * NcJ + gc; /* Trial function basis index */  in PetscFEUpdateElementMat_Hybrid_Internal()
2843         for (gc = 0; gc < NcJ; ++gc) {  in PetscFEUpdateElementMat_Hybrid_Internal()
2844           const PetscInt gidx = g * NcJ + gc;           /* Trial function basis index */  in PetscFEUpdateElementMat_Hybrid_Internal()
2848           elemMat[fOff] += tmpBasisI[fidx] * g0[fc * NcJ + gc] * tmpBasisJ[gidx];  in PetscFEUpdateElementMat_Hybrid_Internal()
2850 …   elemMat[fOff] += tmpBasisI[fidx] * g1[(fc * NcJ + gc) * dE + df] * tmpBasisDerJ[gidx * dE + df];  in PetscFEUpdateElementMat_Hybrid_Internal()
2851 …   elemMat[fOff] += tmpBasisDerI[fidx * dE + df] * g2[(fc * NcJ + gc) * dE + df] * tmpBasisJ[gidx];  in PetscFEUpdateElementMat_Hybrid_Internal()
2852 … dE; ++dg) elemMat[fOff] += tmpBasisDerI[fidx * dE + df] * g3[((fc * NcJ + gc) * dE + df) * dE + d…  in PetscFEUpdateElementMat_Hybrid_Internal()
2949         for (PetscInt gc = 0; gc < NcJ; ++gc) {
2950           const PetscInt gidx = g * NcJ + gc; /* Trial function basis index */
2953             if (g0[fc * NcJ + gc]) {
2955               g0_idxs[4 * g0_i + 1] = fc * NcJ + gc;
2964               if (g1[(fc * NcJ + gc) * dE + df]) {
2966                 g1_idxs[4 * g1_i + 1] = (fc * NcJ + gc) * dE + df;
2973               if (g2[(fc * NcJ + gc) * dE + df]) {
2975                 g2_idxs[4 * g2_i + 1] = (fc * NcJ + gc) * dE + df;
2983                 if (g3[((fc * NcJ + gc) * dE + df) * dE + dg]) {
2985                   g3_idxs[4 * g3_i + 1] = ((fc * NcJ + gc) * dE + df) * dE + dg;