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()
440   PetscCall(PetscFree4(v, J, invJ, detJ));  in DMFieldEvaluate_DS()
507           vol += geom->detJ[i * Nq + k] * weights[k];  in DMFieldEvaluateFV_DS()
508 …for (j = 0; j < Nc; j++) sB[i * Nc + j] += geom->detJ[i * Nq + k] * weights[k] * sqB[(i * Nq + k) …  in DMFieldEvaluateFV_DS()
521           vol += geom->detJ[i * Nq + k] * weights[k];  in DMFieldEvaluateFV_DS()
540           vol += geom->detJ[i * Nq + k] * weights[k];  in DMFieldEvaluateFV_DS()
547 …for (l = 0; l < dimC; l++) sD[(i * Nc + j) * dimC + l] += geom->detJ[i * Nq + k] * weights[k] * pD…  in DMFieldEvaluateFV_DS()
561           vol += geom->detJ[i * Nq + k] * weights[k];  in DMFieldEvaluateFV_DS()
568 …for (l = 0; l < dimC; l++) rD[(i * Nc + j) * dimC + l] += geom->detJ[i * Nq + k] * weights[k] * pD…  in DMFieldEvaluateFV_DS()
589           vol += geom->detJ[i * Nq + k] * weights[k];  in DMFieldEvaluateFV_DS()
606 … 0; m < dimC; m++) sH[(i * Nc + j) * dimC * dimC + l * dimC + m] += geom->detJ[i * Nq + k] * weigh…  in DMFieldEvaluateFV_DS()
623           vol += geom->detJ[i * Nq + k] * weights[k];  in DMFieldEvaluateFV_DS()
640 … 0; m < dimC; m++) rH[(i * Nc + j) * dimC * dimC + l * dimC + m] += geom->detJ[i * Nq + k] * weigh…  in DMFieldEvaluateFV_DS()
896         for (q = 0; q < Nq; q++) geom->suppDetJ[s][p * Nq + q] = cellGeom->detJ[offset * Nq + q];  in DMFieldComputeFaceData_DS()
1071       PetscReal J[9], detJ;  in DMFieldComputeFaceData_DS()  local
1076       PetscCall(DMPlexComputeCellGeometryFEM(K, face, NULL, v0, J, NULL, &detJ));  in DMFieldComputeFaceData_DS()
1108 …          for (q = 0; q < Nq; q++) geom->suppDetJ[s][p * Nq + q] = cellGeom->detJ[offset * Nq + q];  in DMFieldComputeFaceData_DS()