801 …PetscCheck(!fem->T || k <= fem->T->K || (!fem->T->cdim && !fem->T->K), PetscObjectComm((PetscObjec…  in PetscFEGetCellTabulation()
984   PetscInt       cdim; /* Reference coordinate dimension */  in PetscFECreateTabulation()  local
997   PetscCall(DMGetDimension(dm, &cdim));  in PetscFECreateTabulation()
1001   (*T)->K    = !cdim ? 0 : K;  in PetscFECreateTabulation()
1006   (*T)->cdim = cdim;  in PetscFECreateTabulation()
1008 …T)->K; ++k) PetscCall(PetscCalloc1(nrepl * npoints * Nb * Nc * PetscPowInt(cdim, k), &(*T)->T[k]));  in PetscFECreateTabulation()
1051     PetscInt       cdim; /* Reference coordinate dimension */  in PetscFEComputeTabulation()  local
1055     PetscCall(DMGetDimension(dm, &cdim));  in PetscFEComputeTabulation()
1058 …->K == (!cdim ? 0 : K), PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Tabulation K %" PetscInt_FMT " must m…  in PetscFEComputeTabulation()
1061 …ck(T->cdim == cdim, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Tabulation cdim %" PetscInt_FMT " must ma…  in PetscFEComputeTabulation()
2366     const PetscInt   cdim    = T[f]->cdim;  in PetscFEEvaluateFieldJets_Internal()  local
2372     const PetscReal *Dq      = &T[f]->T[1][(r * Nq + q) * Nbf * Ncf * cdim];  in PetscFEEvaluateFieldJets_Internal()
2373     const PetscReal *Hq      = k > 1 ? &T[f]->T[2][(r * Nq + q) * Nbf * Ncf * cdim * cdim] : NULL;  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()
2394 …for (d = 0; d < cdim * cdim; ++d) u_x[hOffset + (fOffset + c) * dE * dE + d] += Hq[cidx * cdim * c…  in PetscFEEvaluateFieldJets_Internal()
2433       const PetscInt  dEt = T->cdim;  in PetscFEEvaluateFieldJets_Hybrid_Internal()
2504   const PetscInt   dEt      = T->cdim;  in PetscFEUpdateElementVec_Internal()
2541   const PetscInt   dE       = T->cdim;  in PetscFEUpdateElementVec_Hybrid_Internal()
2647   const PetscInt   cdim      = TI->cdim;  in PetscFEUpdateElementMat_Internal()  local
2653   const PetscReal *basisDerI = &TI->T[1][(r * NqI + q) * NbI * NcI * cdim];  in PetscFEUpdateElementMat_Internal()
2658   const PetscReal *basisDerJ = &TJ->T[1][(r * NqJ + q) * NbJ * NcJ * cdim];  in PetscFEUpdateElementMat_Internal()
2665 … for (PetscInt df = 0; df < cdim; ++df) tmpBasisDerI[fidx * dE + df] = basisDerI[fidx * cdim + df];  in PetscFEUpdateElementMat_Internal()
2676 … for (PetscInt dg = 0; dg < cdim; ++dg) tmpBasisDerJ[gidx * dE + dg] = basisDerJ[gidx * cdim + dg];  in PetscFEUpdateElementMat_Internal()
2802   const PetscInt   dE        = TI->cdim;  in PetscFEUpdateElementMat_Hybrid_Internal()
2866   PetscInt        dim, cdim, Nq;  in PetscFECreateCellGeometry()  local
2872   PetscCall(DMGetCoordinateDim(dm, &cdim));  in PetscFECreateCellGeometry()
2876   PetscCall(PetscMalloc1(Nq * cdim, &cgeom->v));  in PetscFECreateCellGeometry()
2877   PetscCall(PetscMalloc1(Nq * cdim * cdim, &cgeom->J));  in PetscFECreateCellGeometry()
2878   PetscCall(PetscMalloc1(Nq * cdim * cdim, &cgeom->invJ));  in PetscFECreateCellGeometry()
2881   cgeom->dimEmbed  = cdim;  in PetscFECreateCellGeometry()