37   PetscInt     dimC, dim, numPoints, Nq, p, Nc;  in DMFieldShellEvaluateFEDefault()  local
45   PetscCall(PetscQuadratureGetData(quad, &dim, NULL, &Nq, NULL, NULL));  in DMFieldShellEvaluateFEDefault()
47   PetscCall(PetscMalloc1(dimC * Nq * numPoints, &pfArray));  in DMFieldShellEvaluateFEDefault()
48   for (p = 0; p < numPoints * dimC * Nq; p++) pfArray[p] = geom->v[p];  in DMFieldShellEvaluateFEDefault()
49 …PetscCall(VecCreateMPIWithArray(PetscObjectComm((PetscObject)pointIS), dimC, dimC * Nq * numPoints…  in DMFieldShellEvaluateFEDefault()
57       for (p = 0; p < numPoints * Nq; p++) {  in DMFieldShellEvaluateFEDefault()
74       for (p = 0; p < numPoints * Nq; p++) {  in DMFieldShellEvaluateFEDefault()
94       for (p = 0; p < numPoints * Nq; p++) {  in DMFieldShellEvaluateFEDefault()
117       for (p = 0; p < numPoints * Nq; p++) {  in DMFieldShellEvaluateFEDefault()
150   PetscInt        dimC, dim, numPoints, Nq, p;  in DMFieldShellEvaluateFVDefault()  local
162   PetscCall(PetscQuadratureGetData(quad, NULL, NULL, &Nq, NULL, NULL));  in DMFieldShellEvaluateFVDefault()
163   PetscCheck(Nq == 1, comm, PETSC_ERR_ARG_WRONGSTATE, "quadrature must have only one point");  in DMFieldShellEvaluateFVDefault()