2599             for (PetscInt df = 0; df < _dE; ++df) s += G[df] * tBDJ[df]; \
2617             for (PetscInt df = 0; df < _dE; ++df) s += tBDI[df] * G[df]; \
2635             for (PetscInt df = 0; df < (_dE); ++df) { \
2636 …           for (PetscInt dg = 0; dg < (_dE); ++dg) s += tBDI[df] * G[df * (_dE) + dg] * tBDJ[dg]; \
2665 … for (PetscInt df = 0; df < cdim; ++df) tmpBasisDerI[fidx * dE + df] = basisDerI[fidx * cdim + df];  in PetscFEUpdateElementMat_Internal()  local
2724 …for (PetscInt df = 0; df < dE; ++df) elemMat[fOff] += bI * g1[(fc * NcJ + gc) * dE + df] * tmpBasi…  in PetscFEUpdateElementMat_Internal()
2754 …for (PetscInt df = 0; df < dE; ++df) elemMat[fOff] += tmpBasisDerI[fidx * dE + df] * g2[(fc * NcJ …  in PetscFEUpdateElementMat_Internal()
2784             for (PetscInt df = 0; df < dE; ++df) {  in PetscFEUpdateElementMat_Internal()
2785 …; dg < dE; ++dg) elemMat[fOff] += tmpBasisDerI[fidx * dE + df] * g3[((fc * NcJ + gc) * dE + df) * …  in PetscFEUpdateElementMat_Internal()
2815   PetscInt         f, fc, g, gc, df, dg;  in PetscFEUpdateElementMat_Hybrid_Internal()  local
2822       for (df = 0; df < dE; ++df) tmpBasisDerI[fidx * dE + df] = basisDerI[fidx * dE + df];  in PetscFEUpdateElementMat_Hybrid_Internal()
2849           for (df = 0; df < dE; ++df) {  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 …; dg < dE; ++dg) elemMat[fOff] += tmpBasisDerI[fidx * dE + df] * g3[((fc * NcJ + gc) * dE + df) * …  in PetscFEUpdateElementMat_Hybrid_Internal()
2962           for (PetscInt df = 0; df < dE; ++df) {
2964               if (g1[(fc * NcJ + gc) * dE + df]) {
2966                 g1_idxs[4 * g1_i + 1] = (fc * NcJ + gc) * dE + df;
2967                 g1_idxs[4 * g1_i + 2] = gidx * dE + df;
2973               if (g2[(fc * NcJ + gc) * dE + df]) {
2974                 g2_idxs[4 * g2_i + 0] = fidx * dE + df;
2975                 g2_idxs[4 * g2_i + 1] = (fc * NcJ + gc) * dE + df;
2983                 if (g3[((fc * NcJ + gc) * dE + df) * dE + dg]) {
2984                   g3_idxs[4 * g3_i + 0] = fidx * dE + df;
2985                   g3_idxs[4 * g3_i + 1] = ((fc * NcJ + gc) * dE + df) * dE + dg;