2591     for (PetscInt fc = 0; fc < (_NcI); ++fc) { \
2593         const PetscScalar *G = g1 + (fc * (_NcJ) + gc) * _dE; \
2595           const PetscScalar tBIv = tmpBasisI[f * (_NcI) + fc]; \
2610       for (PetscInt fc = 0; fc < (_NcI); ++fc) { \
2611         const PetscScalar *G = g2 + (fc * (_NcJ) + gc) * _dE; \
2615             const PetscScalar *tBDI = tmpBasisDerI + (f * (_NcI) + fc) * (_dE); \
2627     for (PetscInt fc = 0; fc < (_NcI); ++fc) { \
2629         const PetscScalar *G = g3 + (fc * (_NcJ) + gc) * (_dE) * (_dE); \
2631           const PetscScalar *tBDI = tmpBasisDerI + (f * (_NcI) + fc) * (_dE); \
2661     for (PetscInt fc = 0; fc < NcI; ++fc) {  in PetscFEUpdateElementMat_Internal()  local
2662       const PetscInt fidx = f * NcI + fc; /* Test function basis index */  in PetscFEUpdateElementMat_Internal()
2689       for (PetscInt fc = 0; fc < NcI; ++fc) {  in PetscFEUpdateElementMat_Internal()  local
2690         const PetscScalar bI = tmpBasisI[f * NcI + fc]; /* Test function basis value */  in PetscFEUpdateElementMat_Internal()
2696 …for (PetscInt gc = 0; gc < NcJ; ++gc) elemMat[fOff] += bI * g0[fc * NcJ + gc] * tmpBasisJ[g * NcJ …  in PetscFEUpdateElementMat_Internal()
2714       for (PetscInt fc = 0; fc < NcI; ++fc) {  in PetscFEUpdateElementMat_Internal()
2715         const PetscScalar bI = tmpBasisI[f * NcI + fc]; /* Test function basis value */  in PetscFEUpdateElementMat_Internal()
2724 …for (PetscInt df = 0; df < dE; ++df) elemMat[fOff] += bI * g1[(fc * NcJ + gc) * dE + df] * tmpBasi…  in PetscFEUpdateElementMat_Internal()
2751           for (PetscInt fc = 0; fc < NcI; ++fc) {  in PetscFEUpdateElementMat_Internal()
2752             const PetscInt fidx = f * NcI + fc; /* Test function basis index */  in PetscFEUpdateElementMat_Internal()
2754 …df = 0; df < dE; ++df) elemMat[fOff] += tmpBasisDerI[fidx * dE + df] * g2[(fc * NcJ + gc) * dE + d…  in PetscFEUpdateElementMat_Internal()
2774       for (PetscInt fc = 0; fc < NcI; ++fc) {  in PetscFEUpdateElementMat_Internal()
2775         const PetscInt fidx = f * NcI + fc; /* Test function basis index */  in PetscFEUpdateElementMat_Internal()
2785 …g = 0; dg < dE; ++dg) elemMat[fOff] += tmpBasisDerI[fidx * dE + df] * g3[((fc * NcJ + gc) * dE + d…  in PetscFEUpdateElementMat_Internal()
2815   PetscInt         f, fc, g, gc, df, dg;  in PetscFEUpdateElementMat_Hybrid_Internal()  local
2818     for (fc = 0; fc < NcI; ++fc) {  in PetscFEUpdateElementMat_Hybrid_Internal()
2819       const PetscInt fidx = f * NcI + fc; /* Test function basis index */  in PetscFEUpdateElementMat_Hybrid_Internal()
2839     for (fc = 0; fc < NcI; ++fc) {  in PetscFEUpdateElementMat_Hybrid_Internal()
2840       const PetscInt fidx = f * NcI + fc;           /* Test 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 …for (dg = 0; dg < dE; ++dg) elemMat[fOff] += tmpBasisDerI[fidx * dE + df] * g3[((fc * NcJ + gc) * …  in PetscFEUpdateElementMat_Hybrid_Internal()
2943     for (PetscInt fc = 0; fc < NcI; ++fc) {
2944       const PetscInt fidx = f * NcI + fc; /* Test 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;