2651   const PetscInt   NcI       = TI->Nc;  in PetscFEUpdateElementMat_Internal()  local
2652   const PetscReal *basisI    = &TI->T[0][(r * NqI + q) * NbI * NcI];  in PetscFEUpdateElementMat_Internal()
2653   const PetscReal *basisDerI = &TI->T[1][(r * NqI + q) * NbI * NcI * cdim];  in PetscFEUpdateElementMat_Internal()
2661     for (PetscInt fc = 0; fc < NcI; ++fc) {  in PetscFEUpdateElementMat_Internal()
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()
2690         const PetscScalar bI = tmpBasisI[f * NcI + fc]; /* Test function basis value */  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()
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()
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()
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()
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()
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()
2805   const PetscInt   NcI       = TI->Nc;  in PetscFEUpdateElementMat_Hybrid_Internal()  local
2806   const PetscReal *basisI    = &TI->T[0][(r * NqI + q) * NbI * NcI];  in PetscFEUpdateElementMat_Hybrid_Internal()
2807   const PetscReal *basisDerI = &TI->T[1][(r * NqI + q) * NbI * NcI * dE];  in PetscFEUpdateElementMat_Hybrid_Internal()
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()
2903   const PetscInt NcI     = TI->Nc;
2916     for (PetscInt i = 0; i < NcI * NcJ; i++)
2919     for (PetscInt i = 0; i < NcI * NcJ * dE; i++)
2922     for (PetscInt i = 0; i < NcI * NcJ * dE; i++)
2925     for (PetscInt i = 0; i < NcI * NcJ * dE * dE; i++)
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;
2943     for (PetscInt fc = 0; fc < NcI; ++fc) {
2944       const PetscInt fidx = f * NcI + fc; /* Test function basis index */