88   PetscInt dE;  in PetscFEGeomGetChunk()  local
95   dE                        = geom->dimEmbed;  in PetscFEGeomGetChunk()
102   (*chunkGeom)->v           = PetscSafePointerPlusOffset(geom->v, Nq * dE * cStart);  in PetscFEGeomGetChunk()
103   (*chunkGeom)->J           = PetscSafePointerPlusOffset(geom->J, Nq * dE * dE * cStart);  in PetscFEGeomGetChunk()
104   (*chunkGeom)->invJ        = PetscSafePointerPlusOffset(geom->invJ, Nq * dE * dE * cStart);  in PetscFEGeomGetChunk()
106   (*chunkGeom)->n           = PetscSafePointerPlusOffset(geom->n, Nq * dE * cStart);  in PetscFEGeomGetChunk()
108   (*chunkGeom)->suppJ[0]    = PetscSafePointerPlusOffset(geom->suppJ[0], Nq * dE * dE * cStart);  in PetscFEGeomGetChunk()
109   (*chunkGeom)->suppJ[1]    = PetscSafePointerPlusOffset(geom->suppJ[1], Nq * dE * dE * cStart);  in PetscFEGeomGetChunk()
110   (*chunkGeom)->suppInvJ[0] = PetscSafePointerPlusOffset(geom->suppInvJ[0], Nq * dE * dE * cStart);  in PetscFEGeomGetChunk()
111   (*chunkGeom)->suppInvJ[1] = PetscSafePointerPlusOffset(geom->suppInvJ[1], Nq * dE * dE * cStart);  in PetscFEGeomGetChunk()
163   const PetscInt dE  = geom->dimEmbed;  in PetscFEGeomGetPoint()  local
169   pgeom->dimEmbed = dE;  in PetscFEGeomGetPoint()
174       pgeom->J    = &geom->J[c * Np * dE * dE];  in PetscFEGeomGetPoint()
175       pgeom->invJ = &geom->invJ[c * Np * dE * dE];  in PetscFEGeomGetPoint()
177       pgeom->n    = PetscSafePointerPlusOffset(geom->n, c * Np * dE);  in PetscFEGeomGetPoint()
179 …if (pcoords) CoordinatesRefToReal(dE, dim, pgeom->xi, &geom->v[c * Np * dE], pgeom->J, pcoords, pg…  in PetscFEGeomGetPoint()
181     pgeom->v    = &geom->v[(c * Np + p) * dE];  in PetscFEGeomGetPoint()
182     pgeom->J    = &geom->J[(c * Np + p) * dE * dE];  in PetscFEGeomGetPoint()
183     pgeom->invJ = &geom->invJ[(c * Np + p) * dE * dE];  in PetscFEGeomGetPoint()
185     pgeom->n    = PetscSafePointerPlusOffset(geom->n, (c * Np + p) * dE);  in PetscFEGeomGetPoint()
216   const PetscInt  dE  = geom->dimEmbed;  in PetscFEGeomGetCellPoint()  local
222   pgeom->dimEmbed = dE;  in PetscFEGeomGetCellPoint()
227         pgeom->J    = &geom->suppJ[0][c * Np * dE * dE];  in PetscFEGeomGetCellPoint()
228         pgeom->invJ = &geom->suppInvJ[0][c * Np * dE * dE];  in PetscFEGeomGetCellPoint()
231         pgeom->J    = &geom->J[c * Np * dE * dE];  in PetscFEGeomGetCellPoint()
232         pgeom->invJ = &geom->invJ[c * Np * dE * dE];  in PetscFEGeomGetCellPoint()
238       pgeom->J    = &geom->suppJ[0][(c * Np + p) * dE * dE];  in PetscFEGeomGetCellPoint()
239       pgeom->invJ = &geom->suppInvJ[0][(c * Np + p) * dE * dE];  in PetscFEGeomGetCellPoint()
242       pgeom->J    = &geom->J[(c * Np + p) * dE * dE];  in PetscFEGeomGetCellPoint()
243       pgeom->invJ = &geom->invJ[(c * Np + p) * dE * dE];  in PetscFEGeomGetCellPoint()
262   PetscInt i, j, N, dE;  in PetscFEGeomComplete()  local
266   dE = geom->dimEmbed;  in PetscFEGeomComplete()
267   switch (dE) {  in PetscFEGeomComplete()
270       DMPlex_Det3D_Internal(&geom->detJ[i], &geom->J[dE * dE * i]);  in PetscFEGeomComplete()
271 …if (geom->invJ) DMPlex_Invert3D_Internal(&geom->invJ[dE * dE * i], &geom->J[dE * dE * i], geom->de…  in PetscFEGeomComplete()
276       DMPlex_Det2D_Internal(&geom->detJ[i], &geom->J[dE * dE * i]);  in PetscFEGeomComplete()
277 …if (geom->invJ) DMPlex_Invert2D_Internal(&geom->invJ[dE * dE * i], &geom->J[dE * dE * i], geom->de…  in PetscFEGeomComplete()
289 …for (j = 0; j < dE; j++) geom->n[dE * i + j] = geom->J[dE * dE * i + dE * j + dE - 1] * ((dE == 2)…  in PetscFEGeomComplete()