Lines Matching refs:NcJ
2656 const PetscInt NcJ = TJ->Nc; in PetscFEUpdateElementMat_Internal() local
2657 const PetscReal *basisJ = &TJ->T[0][(r * NqJ + q) * NbJ * NcJ]; in PetscFEUpdateElementMat_Internal()
2658 const PetscReal *basisDerJ = &TJ->T[1][(r * NqJ + q) * NbJ * NcJ * cdim]; in PetscFEUpdateElementMat_Internal()
2672 for (PetscInt gc = 0; gc < NcJ; ++gc) { in PetscFEUpdateElementMat_Internal()
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()
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()
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()
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()
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 …0; df < dE; ++df) elemMat[fOff] += tmpBasisDerI[fidx * dE + df] * g2[(fc * NcJ + gc) * dE + df] * … 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()
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 …; dg < dE; ++dg) elemMat[fOff] += tmpBasisDerI[fidx * dE + df] * g3[((fc * NcJ + gc) * dE + df) * … in PetscFEUpdateElementMat_Internal()
2810 const PetscInt NcJ = TJ->Nc; in PetscFEUpdateElementMat_Hybrid_Internal() local
2811 const PetscReal *basisJ = &TJ->T[0][(r * NqJ + q) * NbJ * NcJ]; in PetscFEUpdateElementMat_Hybrid_Internal()
2812 const PetscReal *basisDerJ = &TJ->T[1][(r * NqJ + q) * NbJ * NcJ * dE]; in PetscFEUpdateElementMat_Hybrid_Internal()
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 …; dg < dE; ++dg) elemMat[fOff] += tmpBasisDerI[fidx * dE + df] * g3[((fc * NcJ + gc) * dE + df) * … in PetscFEUpdateElementMat_Hybrid_Internal()
2905 const PetscInt NcJ = TJ->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;
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;