52 #define NB  2  macro
290   for (kk = 0; kk < NB; kk++) {  in DeformationGradient()
291     for (jj = 0; jj < NB; jj++) {  in DeformationGradient()
292       for (ii = 0; ii < NB; ii++) {  in DeformationGradient()
293         PetscInt    idx  = ii + jj * NB + kk * NB * NB;  in DeformationGradient()
315   PetscInt    idx  = ii + jj * NB + kk * NB * NB;  in DeformationGradientJacobian()
377   for (k = 0; k < NB; k++) {  in FormElements()
378     for (j = 0; j < NB; j++) {  in FormElements()
379       for (i = 0; i < NB; i++) {  in FormElements()
384 …           PetscInt idx = ii + NQ * jj + NQ * NQ * kk + NEQ * i + NEQ * NB * j + NEQ * NB * NB * k;  in FormElements()
421   for (kk = 0; kk < NB; kk++) {  in GatherElementData()
422     for (jj = 0; jj < NB; jj++) {  in GatherElementData()
423       for (ii = 0; ii < NB; ii++) {  in GatherElementData()
424         PetscInt idx = ii + jj * NB + kk * NB * NB;  in GatherElementData()
442   for (kk = 0; kk < NB; kk++) {  in QuadraturePointGeometricJacobian()
443     for (jj = 0; jj < NB; jj++) {  in QuadraturePointGeometricJacobian()
444       for (ii = 0; ii < NB; ii++) {  in QuadraturePointGeometricJacobian()
445         PetscInt idx  = ii + jj * NB + kk * NB * NB;  in QuadraturePointGeometricJacobian()
496           for (kk = 0; kk < NB; kk++) {  in FormElementJacobian()
497             for (jj = 0; jj < NB; jj++) {  in FormElementJacobian()
498               for (ii = 0; ii < NB; ii++) {  in FormElementJacobian()
499                 PetscInt    idx  = ii + jj * NB + kk * NB * NB;  in FormElementJacobian()
511           for (kk = 0; kk < NB; kk++) {  in FormElementJacobian()
512             for (jj = 0; jj < NB; jj++) {  in FormElementJacobian()
513               for (ii = 0; ii < NB; ii++) {  in FormElementJacobian()
514                 PetscInt idx  = ii + jj * NB + kk * NB * NB;  in FormElementJacobian()
525           for (k = 0; k < NB; k++) {  in FormElementJacobian()
526             for (j = 0; j < NB; j++) {  in FormElementJacobian()
527               for (i = 0; i < NB; i++) {  in FormElementJacobian()
529                   PetscInt tridx = l + 3 * (i + j * NB + k * NB * NB);  in FormElementJacobian()
536                   for (kk = 0; kk < NB; kk++) {  in FormElementJacobian()
537                     for (jj = 0; jj < NB; jj++) {  in FormElementJacobian()
538                       for (ii = 0; ii < NB; ii++) {  in FormElementJacobian()
539                         PetscInt    idx  = ii + jj * NB + kk * NB * NB;  in FormElementJacobian()
544                           PetscInt teidx = ll + 3 * (ii + jj * NB + kk * NB * NB);  in FormElementJacobian()
563   for (kk = 0; kk < NB; kk++) {  in ApplyBCsElement()
564     for (jj = 0; jj < NB; jj++) {  in ApplyBCsElement()
565       for (ii = 0; ii < NB; ii++) {  in ApplyBCsElement()
567           PetscInt tridx = ll + 3 * (ii + jj * NB + kk * NB * NB);  in ApplyBCsElement()
568           for (ek = 0; ek < NB; ek++) {  in ApplyBCsElement()
569             for (ej = 0; ej < NB; ej++) {  in ApplyBCsElement()
570               for (ei = 0; ei < NB; ei++) {  in ApplyBCsElement()
573                     PetscInt teidx = el + 3 * (ei + ej * NB + ek * NB * NB);  in ApplyBCsElement()
637         for (kk = 0; kk < NB; kk++) {  in FormJacobianLocal()
638           for (jj = 0; jj < NB; jj++) {  in FormJacobianLocal()
639             for (ii = 0; ii < NB; ii++) {  in FormJacobianLocal()
752         for (kk = 0; kk < NB; kk++) {  in FormFunctionLocal()
753           for (jj = 0; jj < NB; jj++) {  in FormFunctionLocal()
754             for (ii = 0; ii < NB; ii++) {  in FormFunctionLocal()
755               PetscInt idx = ii + jj * NB + kk * NB * NB;  in FormFunctionLocal()
840         for (kk = 0; kk < NB; kk++) {  in TangentLoad()
841           for (jj = 0; jj < NB; jj++) {  in TangentLoad()
842             for (ii = 0; ii < NB; ii++) {  in TangentLoad()
843               PetscInt idx = ii + jj * NB + kk * NB * NB;  in TangentLoad()