2673         const PetscInt gidx = g * NcJ + gc; /* Trial function basis index */  in PetscFEUpdateElementMat_Internal()  local
2675         tmpBasisJ[gidx] = basisJ[gidx];  in PetscFEUpdateElementMat_Internal()
2676 … for (PetscInt dg = 0; dg < cdim; ++dg) tmpBasisDerJ[gidx * dE + dg] = basisDerJ[gidx * cdim + dg];  in PetscFEUpdateElementMat_Internal()
2722             const PetscInt gidx = g * NcJ + gc; /* Trial function basis index */  in PetscFEUpdateElementMat_Internal()
2724 … df < dE; ++df) elemMat[fOff] += bI * g1[(fc * NcJ + gc) * dE + df] * tmpBasisDerJ[gidx * dE + df];  in PetscFEUpdateElementMat_Internal()
2782             const PetscInt gidx = g * NcJ + gc; /* Trial function basis index */  in PetscFEUpdateElementMat_Internal()
2785 …sisDerI[fidx * dE + df] * g3[((fc * NcJ + gc) * dE + df) * dE + dg] * tmpBasisDerJ[gidx * dE + dg];  in PetscFEUpdateElementMat_Internal()
2829       const PetscInt gidx = g * NcJ + gc; /* Trial function basis index */  in PetscFEUpdateElementMat_Hybrid_Internal()  local
2831       tmpBasisJ[gidx] = basisJ[gidx];  in PetscFEUpdateElementMat_Hybrid_Internal()
2832       for (dg = 0; dg < dE; ++dg) tmpBasisDerJ[gidx * dE + dg] = basisDerJ[gidx * dE + dg];  in PetscFEUpdateElementMat_Hybrid_Internal()
2844           const PetscInt gidx = g * NcJ + gc;           /* Trial function basis index */  in PetscFEUpdateElementMat_Hybrid_Internal()  local
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 …sisDerI[fidx * dE + df] * g3[((fc * NcJ + gc) * dE + df) * dE + dg] * tmpBasisDerJ[gidx * dE + dg];  in PetscFEUpdateElementMat_Hybrid_Internal()
2950           const PetscInt gidx = g * NcJ + gc; /* Trial function basis index */
2956               g0_idxs[4 * g0_i + 2] = gidx;
2967                 g1_idxs[4 * g1_i + 2] = gidx * dE + df;
2976                 g2_idxs[4 * g2_i + 2] = gidx;
2986                   g3_idxs[4 * g3_i + 2] = gidx * dE + dg;