2359   PetscInt dOffset = 0, fOffset = 0, f, g;  in PetscFEEvaluateFieldJets_Internal()  local
2377     for (c = 0; c < Ncf; ++c) u[fOffset + c] = 0.0;  in PetscFEEvaluateFieldJets_Internal()
2378     for (d = 0; d < dE * Ncf; ++d) u_x[fOffset * dE + d] = 0.0;  in PetscFEEvaluateFieldJets_Internal()
2383         u[fOffset + c] += Bq[cidx] * coefficients[dOffset + b];  in PetscFEEvaluateFieldJets_Internal()
2384 …for (d = 0; d < cdim; ++d) u_x[(fOffset + c) * dE + d] += Dq[cidx * cdim + d] * coefficients[dOffs…  in PetscFEEvaluateFieldJets_Internal()
2389       for (d = 0; d < dE * dE * Ncf; ++d) u_x[hOffset + fOffset * dE * dE + d] = 0.0;  in PetscFEEvaluateFieldJets_Internal()
2394 …for (d = 0; d < cdim * cdim; ++d) u_x[hOffset + (fOffset + c) * dE * dE + d] += Hq[cidx * cdim * c…  in PetscFEEvaluateFieldJets_Internal()
2397       PetscCall(PetscFEPushforwardHessian(fe, fegeom, 1, &u_x[hOffset + fOffset * dE * dE]));  in PetscFEEvaluateFieldJets_Internal()
2399     PetscCall(PetscFEPushforward(fe, fegeom, 1, &u[fOffset]));  in PetscFEEvaluateFieldJets_Internal()
2400     PetscCall(PetscFEPushforwardGradient(fe, fegeom, 1, &u_x[fOffset * dE]));  in PetscFEEvaluateFieldJets_Internal()
2402       for (c = 0; c < Ncf; ++c) u_t[fOffset + c] = 0.0;  in PetscFEEvaluateFieldJets_Internal()
2407           u_t[fOffset + c] += Bq[cidx] * coefficients_t[dOffset + b];  in PetscFEEvaluateFieldJets_Internal()
2410       PetscCall(PetscFEPushforward(fe, fegeom, 1, &u_t[fOffset]));  in PetscFEEvaluateFieldJets_Internal()
2412     fOffset += Ncf;  in PetscFEEvaluateFieldJets_Internal()
2420   PetscInt dOffset = 0, fOffset = 0, f;  in PetscFEEvaluateFieldJets_Hybrid_Internal()  local
2448         for (c = 0; c < Ncf; ++c) u[fOffset + c] = 0.0;  in PetscFEEvaluateFieldJets_Hybrid_Internal()
2449         for (d = 0; d < dE * Ncf; ++d) u_x[fOffset * dE + d] = 0.0;  in PetscFEEvaluateFieldJets_Hybrid_Internal()
2454             u[fOffset + c] += Bq[cidx] * coefficients[dOffset + b];  in PetscFEEvaluateFieldJets_Hybrid_Internal()
2455 …for (d = 0; d < dEt; ++d) u_x[(fOffset + c) * dE + d] += Dq[cidx * dEt + d] * coefficients[dOffset…  in PetscFEEvaluateFieldJets_Hybrid_Internal()
2458         PetscCall(PetscFEPushforward(fe, isCohesive ? fegeom : &fegeomNbr[s], 1, &u[fOffset]));  in PetscFEEvaluateFieldJets_Hybrid_Internal()
2459 …scCall(PetscFEPushforwardGradient(fe, isCohesive ? fegeom : &fegeomNbr[s], 1, &u_x[fOffset * dE]));  in PetscFEEvaluateFieldJets_Hybrid_Internal()
2461           for (c = 0; c < Ncf; ++c) u_t[fOffset + c] = 0.0;  in PetscFEEvaluateFieldJets_Hybrid_Internal()
2466               u_t[fOffset + c] += Bq[cidx] * coefficients_t[dOffset + b];  in PetscFEEvaluateFieldJets_Hybrid_Internal()
2469           PetscCall(PetscFEPushforward(fe, fegeom, 1, &u_t[fOffset]));  in PetscFEEvaluateFieldJets_Hybrid_Internal()
2471         fOffset += Ncf;  in PetscFEEvaluateFieldJets_Hybrid_Internal()