246   PetscReal         *v, *J, *invJ, *detJ;  in DMFieldEvaluate_DS()  local
269 …herMax * dimR, &v, gatherMax * dimR * dimR, &J, gatherMax * dimR * dimR, &invJ, gatherMax, &detJ));  in DMFieldEvaluate_DS()
300       PetscCall(DMPlexComputeCellGeometryFEM(field->dm, c, quad, v, J, invJ, detJ));  in DMFieldEvaluate_DS()
325 …for (e = 0; e < dimR; e++) vs[d] += invJ[dimR * dimR * p + e * dimR + d] * cD[(nc * p + g) * dimR …  in DMFieldEvaluate_DS()
338 …for (e = 0; e < dimR; e++) v[d] += invJ[dimR * dimR * p + e * dimR + d] * cD[(nc * p + g) * dimR +…  in DMFieldEvaluate_DS()
356 …for (e = 0; e < dimR; e++) vs[d] += invJ[dimR * dimR * p + e * dimR + d] * cH[(nc * dimR * p + g) …  in DMFieldEvaluate_DS()
366 …for (e = 0; e < dimR; e++) vs[d] += invJ[dimR * dimR * p + e * dimR + d] * cH[((nc * p + g) * dimR…  in DMFieldEvaluate_DS()
380 …for (e = 0; e < dimR; e++) v[d] += invJ[dimR * dimR * p + e * dimR + d] * cH[(nc * dimR * p + g) *…  in DMFieldEvaluate_DS()
388 …for (e = 0; e < dimR; e++) v[d] += invJ[dimR * dimR * p + e * dimR + d] * cH[((nc * p + g) * dimR …  in DMFieldEvaluate_DS()
440   PetscCall(PetscFree4(v, J, invJ, detJ));  in DMFieldEvaluate_DS()
545 …for (m = 0; m < dim; m++) pD[l] += geom->invJ[((i * Nq + k) * dimC + m) * dimC + l] * sqD[((i * Nq…  in DMFieldEvaluateFV_DS()
566 …for (m = 0; m < dim; m++) pD[l] += geom->invJ[((i * Nq + k) * dimC + m) * dimC + l] * rqD[((i * Nq…  in DMFieldEvaluateFV_DS()
587           const PetscReal *invJ = &geom->invJ[(i * Nq + k) * dimC * dimC];  in DMFieldEvaluateFV_DS()  local
601 …  for (r = 0; r < dim; r++) pH[l][m] += invJ[q * dimC + l] * invJ[r * dimC + m] * spH[q * dim + r];  in DMFieldEvaluateFV_DS()
621           const PetscReal *invJ = &geom->invJ[(i * Nq + k) * dimC * dimC];  in DMFieldEvaluateFV_DS()  local
635 …  for (r = 0; r < dim; r++) pH[l][m] += invJ[q * dimC + l] * invJ[r * dimC + m] * rpH[q * dim + r];  in DMFieldEvaluateFV_DS()
894           geom->suppInvJ[s][p * Nq * dE * dE + q] = cellGeom->invJ[offset * Nq * dE * dE + q];  in DMFieldComputeFaceData_DS()
1106                 geom->suppInvJ[s][p * Nq * dE * dE + q] = cellGeom->invJ[offset * Nq * dE * dE + q];  in DMFieldComputeFaceData_DS()