576     PetscInt i, j, k, l, m, q, r;  in DMFieldEvaluateFV_DS()  local
600                 for (q = 0; q < dim; q++) {  in DMFieldEvaluateFV_DS()
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()
634                 for (q = 0; q < dim; q++) {  in DMFieldEvaluateFV_DS()
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()
829       PetscInt Np = geom->numPoints, q, dE = geom->dimEmbed, d;  in DMFieldComputeFaceData_DS()  local
831       for (q = 0; q < Np; ++q)  in DMFieldComputeFaceData_DS()
832         for (d = 0; d < dE; ++d) geom->n[(p * Np + q) * dE + d] = -geom->n[(p * Np + q) * dE + d];  in DMFieldComputeFaceData_DS()
886       PetscInt        numSupp, s, q;  in DMFieldComputeFaceData_DS()  local
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()
913     PetscInt        f, i, o, q, s;  in DMFieldComputeFaceData_DS()  local
999         PetscInt q;  in DMFieldComputeFaceData_DS()  local
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()
1019             lambda[0] = -(geom->xi[2 * q] + geom->xi[2 * q + 1]) / 2.;  in DMFieldComputeFaceData_DS()
1020             lambda[1] = (geom->xi[2 * q] + 1.) / 2.;  in DMFieldComputeFaceData_DS()
1021             lambda[2] = (geom->xi[2 * q + 1] + 1.) / 2.;  in DMFieldComputeFaceData_DS()
1028             orientPoints[o][2 * q + 0] = -(lambdao[0] + lambdao[2]) + lambdao[1];  in DMFieldComputeFaceData_DS()
1029             orientPoints[o][2 * q + 1] = -(lambdao[0] + lambdao[1]) + lambdao[2];  in DMFieldComputeFaceData_DS()
1033           for (q = 0; q < Nq; q++) {  in DMFieldComputeFaceData_DS()
1037             xi[0] = geom->xi[2 * q];  in DMFieldComputeFaceData_DS()
1038             xi[1] = geom->xi[2 * q + 1];  in DMFieldComputeFaceData_DS()
1059             orientPoints[o][2 * q + 0] = xio[0];  in DMFieldComputeFaceData_DS()
1060             orientPoints[o][2 * q + 1] = xio[1];  in DMFieldComputeFaceData_DS()
1083         for (q = 0; q < Nq; q++) {  in DMFieldComputeFaceData_DS()
1086 … CoordinatesRefToReal(dE, dim, xi0, v0, J, &orientPoints[o][dim * q + 0], &cellPoints[dE * q + 0]);  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()