624 … totDim, totDimIn, totDimAux = 0, cOffset = 0, cOffsetIn = 0, cOffsetAux = 0, fOffset, e, NcI, NcS;  in PetscFEIntegrateHybridResidual_Basic()  local
665   NcS = NcI;  in PetscFEIntegrateHybridResidual_Basic()
668     NcS *= 2;  in PetscFEIntegrateHybridResidual_Basic()
683     PetscCall(PetscArrayzero(f0, Nq * NcS));  in PetscFEIntegrateHybridResidual_Basic()
684     PetscCall(PetscArrayzero(f1, Nq * NcS * dE));  in PetscFEIntegrateHybridResidual_Basic()
710 …_t, u_x, aOff, aOff_x, a, NULL, a_x, t, fegeom.v, fegeom.n, numConstants, constants, &f0[q * NcS]);  in PetscFEIntegrateHybridResidual_Basic()
711       for (c = 0; c < NcS; ++c) f0[q * NcS + c] *= w;  in PetscFEIntegrateHybridResidual_Basic()
712 …_x, aOff, aOff_x, a, NULL, a_x, t, fegeom.v, fegeom.n, numConstants, constants, &f1[q * NcS * dE]);  in PetscFEIntegrateHybridResidual_Basic()
713       for (c = 0; c < NcS; ++c)  in PetscFEIntegrateHybridResidual_Basic()
714         for (d = 0; d < dE; ++d) f1[(q * NcS + c) * dE + d] *= w;  in PetscFEIntegrateHybridResidual_Basic()
722         for (c = 0; c < NcS; ++c) {  in PetscFEIntegrateHybridResidual_Basic()
723 …rintf(PETSC_COMM_SELF, "  f0[%" PetscInt_FMT "] %g\n", c, (double)PetscRealPart(f0[q * NcS + c])));  in PetscFEIntegrateHybridResidual_Basic()
725 …" PetscInt_FMT ",%" PetscInt_FMT "] %g", c, d, (double)PetscRealPart(f1[(q * NcS + c) * dE + d])));  in PetscFEIntegrateHybridResidual_Basic()
1078   PetscInt            NcI = 0, NcJ = 0, NcS, NcT;  in PetscFEIntegrateHybridJacobian_Basic()  local
1134   NcS = isCohesiveFieldI ? NcI : 2 * NcI;  in PetscFEIntegrateHybridJacobian_Basic()
1138     NcS *= 2;  in PetscFEIntegrateHybridJacobian_Basic()
1146   PetscCall(PetscArrayzero(g0, NcS * NcT));  in PetscFEIntegrateHybridJacobian_Basic()
1147   PetscCall(PetscArrayzero(g1, NcS * NcT * dim));  in PetscFEIntegrateHybridJacobian_Basic()
1148   PetscCall(PetscArrayzero(g2, NcS * NcT * dim));  in PetscFEIntegrateHybridJacobian_Basic()
1149   PetscCall(PetscArrayzero(g3, NcS * NcT * dim * dim));  in PetscFEIntegrateHybridJacobian_Basic()
1181         PetscCall(PetscArrayzero(g0, NcS * NcT));  in PetscFEIntegrateHybridJacobian_Basic()
1183         for (c = 0; c < NcS * NcT; ++c) g0[c] *= w;  in PetscFEIntegrateHybridJacobian_Basic()
1186         PetscCall(PetscArrayzero(g1, NcS * NcT * dim));  in PetscFEIntegrateHybridJacobian_Basic()
1188         for (c = 0; c < NcS * NcT * dim; ++c) g1[c] *= w;  in PetscFEIntegrateHybridJacobian_Basic()
1191         PetscCall(PetscArrayzero(g2, NcS * NcT * dim));  in PetscFEIntegrateHybridJacobian_Basic()
1193         for (c = 0; c < NcS * NcT * dim; ++c) g2[c] *= w;  in PetscFEIntegrateHybridJacobian_Basic()
1196         PetscCall(PetscArrayzero(g3, NcS * NcT * dim * dim));  in PetscFEIntegrateHybridJacobian_Basic()
1198         for (c = 0; c < NcS * NcT * dim * dim; ++c) g3[c] *= w;  in PetscFEIntegrateHybridJacobian_Basic()