2589 #define petsc_elemmat_kernel_g1(_NbI, _NcI, _NbJ, _NcJ, _dE) \  argument
2592       for (PetscInt gc = 0; gc < (_NcJ); ++gc) { \
2593         const PetscScalar *G = g1 + (fc * (_NcJ) + gc) * _dE; \
2597             const PetscScalar *tBDJ = tmpBasisDerJ + (g * (_NcJ) + gc) * (_dE); \
2607 #define petsc_elemmat_kernel_g2(_NbI, _NcI, _NbJ, _NcJ, _dE) \  argument
2609     for (PetscInt gc = 0; gc < (_NcJ); ++gc) { \
2611         const PetscScalar *G = g2 + (fc * (_NcJ) + gc) * _dE; \
2613           const PetscScalar tBJv = tmpBasisJ[g * (_NcJ) + gc]; \
2625 #define petsc_elemmat_kernel_g3(_NbI, _NcI, _NbJ, _NcJ, _dE) \  argument
2628       for (PetscInt gc = 0; gc < (_NcJ); ++gc) { \
2629         const PetscScalar *G = g3 + (fc * (_NcJ) + gc) * (_dE) * (_dE); \
2634             const PetscScalar *tBDJ = tmpBasisDerJ + (g * (_NcJ) + gc) * (_dE); \