458   PetscInt         Nq, Nc, dimC, qNc, N;  in DMFieldEvaluateFV_DS()  local
488   PetscCall(PetscQuadratureGetData(quad, NULL, &qNc, &Nq, NULL, &weights));  in DMFieldEvaluateFV_DS()
490   N = numPoints * Nq * Nc;  in DMFieldEvaluateFV_DS()
506         for (k = 0; k < Nq; k++) {  in DMFieldEvaluateFV_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()
520         for (k = 0; k < Nq; k++) {  in DMFieldEvaluateFV_DS()
521           vol += geom->detJ[i * Nq + k] * weights[k];  in DMFieldEvaluateFV_DS()
522           for (j = 0; j < Nc; j++) rB[i * Nc + j] += weights[k] * rqB[(i * Nq + k) * Nc + j];  in DMFieldEvaluateFV_DS()
539         for (k = 0; k < Nq; k++) {  in DMFieldEvaluateFV_DS()
540           vol += geom->detJ[i * Nq + k] * weights[k];  in DMFieldEvaluateFV_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()
547 …for (l = 0; l < dimC; l++) sD[(i * Nc + j) * dimC + l] += geom->detJ[i * Nq + k] * weights[k] * pD…  in DMFieldEvaluateFV_DS()
560         for (k = 0; k < Nq; k++) {  in DMFieldEvaluateFV_DS()
561           vol += geom->detJ[i * Nq + k] * weights[k];  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()
568 …for (l = 0; l < dimC; l++) rD[(i * Nc + j) * dimC + l] += geom->detJ[i * Nq + k] * weights[k] * pD…  in DMFieldEvaluateFV_DS()
586         for (k = 0; k < Nq; k++) {  in DMFieldEvaluateFV_DS()
587           const PetscReal *invJ = &geom->invJ[(i * Nq + k) * dimC * dimC];  in DMFieldEvaluateFV_DS()
589           vol += geom->detJ[i * Nq + k] * weights[k];  in DMFieldEvaluateFV_DS()
596             const PetscScalar *spH = &sqH[((i * Nq + k) * Nc + j) * dimC * dimC];  in DMFieldEvaluateFV_DS()
606 …dimC; m++) sH[(i * Nc + j) * dimC * dimC + l * dimC + m] += geom->detJ[i * Nq + k] * weights[k] * …  in DMFieldEvaluateFV_DS()
620         for (k = 0; k < Nq; k++) {  in DMFieldEvaluateFV_DS()
621           const PetscReal *invJ = &geom->invJ[(i * Nq + k) * dimC * dimC];  in DMFieldEvaluateFV_DS()
623           vol += geom->detJ[i * Nq + k] * weights[k];  in DMFieldEvaluateFV_DS()
630             const PetscReal *rpH = &rqH[((i * Nq + k) * Nc + j) * dimC * dimC];  in DMFieldEvaluateFV_DS()
640 …dimC; m++) rH[(i * Nc + j) * dimC * dimC + l * dimC + m] += geom->detJ[i * Nq + k] * weights[k] * …  in DMFieldEvaluateFV_DS()
789   PetscInt        p, dim, dE, numFaces, Nq;  in DMFieldComputeFaceData_DS()  local
803   Nq       = geom->numPoints;  in DMFieldComputeFaceData_DS()
892         for (q = 0; q < Nq * dE * dE; q++) {  in DMFieldComputeFaceData_DS()
893           geom->suppJ[s][p * Nq * dE * dE + q]    = cellGeom->J[offset * Nq * dE * dE + q];  in DMFieldComputeFaceData_DS()
894           geom->suppInvJ[s][p * Nq * dE * dE + q] = cellGeom->invJ[offset * Nq * dE * dE + q];  in DMFieldComputeFaceData_DS()
896         for (q = 0; q < Nq; q++) geom->suppDetJ[s][p * Nq + q] = cellGeom->detJ[offset * Nq + q];  in DMFieldComputeFaceData_DS()
935     PetscCall(PetscMalloc1(dE * Nq, &cellPoints));  in DMFieldComputeFaceData_DS()
936     PetscCall(PetscMalloc1(Nq, &unusedWeights));  in DMFieldComputeFaceData_DS()
938     PetscCall(PetscQuadratureSetData(cellQuad, dE, 1, Nq, cellPoints, unusedWeights));  in DMFieldComputeFaceData_DS()
1001         PetscCall(PetscMalloc1(Nq * dim, &orientPoints[o]));  in DMFieldComputeFaceData_DS()
1008             for (q = 0; q < Nq; q++) orientPoints[o][q] = -geom->xi[q];  in DMFieldComputeFaceData_DS()
1010             for (q = 0; q < Nq; q++) orientPoints[o][q] = geom->xi[q];  in DMFieldComputeFaceData_DS()
1014           for (q = 0; q < Nq; q++) {  in DMFieldComputeFaceData_DS()
1033           for (q = 0; q < Nq; q++) {  in DMFieldComputeFaceData_DS()
1083         for (q = 0; q < Nq; q++) {  in DMFieldComputeFaceData_DS()
1104               for (q = 0; q < Nq * dE * dE; q++) {  in DMFieldComputeFaceData_DS()
1105                 geom->suppJ[s][p * Nq * dE * dE + q]    = cellGeom->J[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()
1108 …          for (q = 0; q < Nq; q++) geom->suppDetJ[s][p * Nq + q] = cellGeom->detJ[offset * Nq + q];  in DMFieldComputeFaceData_DS()