151   PetscInt         nq, nc, dim, meshDim, numCells, feDim, i, K = H ? 2 : (D ? 1 : (B ? 0 : -1));  in DMFieldEvaluateFE_DS()  local
163   PetscCall(PetscQuadratureGetData(quad, &dim, NULL, &nq, &qpoints, NULL));  in DMFieldEvaluateFE_DS()
177   PetscCall(PetscFECreateTabulation(fe, 1, nq, qpoints, K, &T));  in DMFieldEvaluateFE_DS()
189         PetscScalar *cB = &((PetscScalar *)B)[nc * nq * i];  in DMFieldEvaluateFE_DS()
191         DMFieldDSdot(cB, T->T[0], elem, nq, feDim, nc, (PetscScalar));  in DMFieldEvaluateFE_DS()
193         PetscReal *cB = &((PetscReal *)B)[nc * nq * i];  in DMFieldEvaluateFE_DS()
195         DMFieldDSdot(cB, T->T[0], elem, nq, feDim, nc, PetscRealPart);  in DMFieldEvaluateFE_DS()
200         PetscScalar *cD = &((PetscScalar *)D)[nc * nq * dim * i];  in DMFieldEvaluateFE_DS()
202         DMFieldDSdot(cD, T->T[1], elem, nq, feDim, (nc * dim), (PetscScalar));  in DMFieldEvaluateFE_DS()
204         PetscReal *cD = &((PetscReal *)D)[nc * nq * dim * i];  in DMFieldEvaluateFE_DS()
206         DMFieldDSdot(cD, T->T[1], elem, nq, feDim, (nc * dim), PetscRealPart);  in DMFieldEvaluateFE_DS()
211         PetscScalar *cH = &((PetscScalar *)H)[nc * nq * dim * dim * i];  in DMFieldEvaluateFE_DS()
213         DMFieldDSdot(cH, T->T[2], elem, nq, feDim, (nc * dim * dim), (PetscScalar));  in DMFieldEvaluateFE_DS()
215         PetscReal *cH = &((PetscReal *)H)[nc * nq * dim * dim * i];  in DMFieldEvaluateFE_DS()
217         DMFieldDSdot(cH, T->T[2], elem, nq, feDim, (nc * dim * dim), PetscRealPart);  in DMFieldEvaluateFE_DS()
281     PetscInt nq = cellDegrees[c], p;  in DMFieldEvaluate_DS()  local
283     if (nq) {  in DMFieldEvaluate_DS()
293       for (p = 0; p < dim * nq; p++) coordsReal[p] = PetscRealPart(cellPoints[dim * offset + p]);  in DMFieldEvaluate_DS()
294       PetscCall(DMPlexCoordinatesToReference(field->dm, c, nq, coordsReal, coordsRef));  in DMFieldEvaluate_DS()
295       PetscCall(PetscFECreateTabulation(cellFE, 1, nq, coordsRef, K, &T));  in DMFieldEvaluate_DS()
297       PetscCall(PetscMalloc1(dimR * nq, &quadPoints));  in DMFieldEvaluate_DS()
298       for (p = 0; p < dimR * nq; p++) quadPoints[p] = coordsRef[p];  in DMFieldEvaluate_DS()
299       PetscCall(PetscQuadratureSetData(quad, dimR, 0, nq, quadPoints, NULL));  in DMFieldEvaluate_DS()
307           DMFieldDSdot(cB, T->T[0], elem, nq, feDim, nc, (PetscScalar));  in DMFieldEvaluate_DS()
311           DMFieldDSdot(cB, T->T[0], elem, nq, feDim, nc, PetscRealPart);  in DMFieldEvaluate_DS()
318           DMFieldDSdot(cD, T->T[1], elem, nq, feDim, (nc * dim), (PetscScalar));  in DMFieldEvaluate_DS()
319           for (p = 0; p < nq; p++) {  in DMFieldEvaluate_DS()
333           DMFieldDSdot(cD, T->T[1], elem, nq, feDim, (nc * dim), PetscRealPart);  in DMFieldEvaluate_DS()
334           for (p = 0; p < nq; p++) {  in DMFieldEvaluate_DS()
349           DMFieldDSdot(cH, T->T[2], elem, nq, feDim, (nc * dim * dim), (PetscScalar));  in DMFieldEvaluate_DS()
350           for (p = 0; p < nq; p++) {  in DMFieldEvaluate_DS()
375           DMFieldDSdot(cH, T->T[2], elem, nq, feDim, (nc * dim * dim), PetscRealPart);  in DMFieldEvaluate_DS()
376           for (p = 0; p < nq; p++) {  in DMFieldEvaluate_DS()
399     offset += nq;  in DMFieldEvaluate_DS()