2589 #define petsc_elemmat_kernel_g1(_NbI, _NcI, _NbJ, _NcJ, _dE) \  argument
2593         const PetscScalar *G = g1 + (fc * (_NcJ) + gc) * _dE; \
2597             const PetscScalar *tBDJ = tmpBasisDerJ + (g * (_NcJ) + gc) * (_dE); \
2599             for (PetscInt df = 0; df < _dE; ++df) s += G[df] * tBDJ[df]; \
2607 #define petsc_elemmat_kernel_g2(_NbI, _NcI, _NbJ, _NcJ, _dE) \  argument
2611         const PetscScalar *G = g2 + (fc * (_NcJ) + gc) * _dE; \
2615             const PetscScalar *tBDI = tmpBasisDerI + (f * (_NcI) + fc) * (_dE); \
2617             for (PetscInt df = 0; df < _dE; ++df) s += tBDI[df] * G[df]; \
2625 #define petsc_elemmat_kernel_g3(_NbI, _NcI, _NbJ, _NcJ, _dE) \  argument
2629         const PetscScalar *G = g3 + (fc * (_NcJ) + gc) * (_dE) * (_dE); \
2631           const PetscScalar *tBDI = tmpBasisDerI + (f * (_NcI) + fc) * (_dE); \
2634             const PetscScalar *tBDJ = tmpBasisDerJ + (g * (_NcJ) + gc) * (_dE); \
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]; \