86 …tClosure_Internal(DMField field, PetscInt cell, PetscBool *isDG, PetscInt *Nc, const PetscScalar *…  in DMFieldGetClosure_Internal()  argument
96   *Nc     = 0;  in DMFieldGetClosure_Internal()
106   PetscCall(PetscSectionGetDof(s, cell, Nc));  in DMFieldGetClosure_Internal()
107   if (!*Nc) goto cg;  in DMFieldGetClosure_Internal()
113   PetscCall(DMGetWorkArray(fdm, *Nc, MPIU_SCALAR, values));  in DMFieldGetClosure_Internal()
114   PetscCall(PetscArraycpy(*values, cvalues, *Nc));  in DMFieldGetClosure_Internal()
123   PetscCall(DMPlexVecGetClosure(fdm, s, dsfield->vec, cell, Nc, values));  in DMFieldGetClosure_Internal()
127 …eClosure_Internal(DMField field, PetscInt cell, PetscBool *isDG, PetscInt *Nc, const PetscScalar *…  in DMFieldRestoreClosure_Internal()  argument
135     PetscCall(DMRestoreWorkArray(dsfield->dmDG, *Nc, MPIU_SCALAR, values));  in DMFieldRestoreClosure_Internal()
139     PetscCall(DMPlexVecRestoreClosure(fdm, s, dsfield->vec, cell, Nc, values));  in DMFieldRestoreClosure_Internal()
458   PetscInt         Nq, Nc, dimC, qNc, N;  in DMFieldEvaluateFV_DS()  local
467   Nc = field->numComponents;  in DMFieldEvaluateFV_DS()
490   N = numPoints * Nq * Nc;  in DMFieldEvaluateFV_DS()
505         for (j = 0; j < Nc; j++) sB[i * Nc + j] = 0.;  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()
510         for (k = 0; k < Nc; k++) sB[i * Nc + k] /= vol;  in DMFieldEvaluateFV_DS()
519         for (j = 0; j < Nc; j++) rB[i * Nc + j] = 0.;  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()
524         for (k = 0; k < Nc; k++) rB[i * Nc + k] /= vol;  in DMFieldEvaluateFV_DS()
538         for (j = 0; j < Nc * dimC; j++) sD[i * Nc * dimC + j] = 0.;  in DMFieldEvaluateFV_DS()
541           for (j = 0; j < Nc; j++) {  in DMFieldEvaluateFV_DS()
545 …pD[l] += geom->invJ[((i * Nq + k) * dimC + m) * dimC + l] * sqD[((i * Nq + k) * Nc + j) * dim + m];  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()
550         for (k = 0; k < Nc * dimC; k++) sD[i * Nc * dimC + k] /= vol;  in DMFieldEvaluateFV_DS()
559         for (j = 0; j < Nc * dimC; j++) rD[i * Nc * dimC + j] = 0.;  in DMFieldEvaluateFV_DS()
562           for (j = 0; j < Nc; j++) {  in DMFieldEvaluateFV_DS()
566 …pD[l] += geom->invJ[((i * Nq + k) * dimC + m) * dimC + l] * rqD[((i * Nq + k) * Nc + j) * dim + m];  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()
571         for (k = 0; k < Nc * dimC; k++) rD[i * Nc * dimC + k] /= vol;  in DMFieldEvaluateFV_DS()
585         for (j = 0; j < Nc * dimC * dimC; j++) sH[i * Nc * dimC * dimC + j] = 0.;  in DMFieldEvaluateFV_DS()
590           for (j = 0; j < Nc; j++) {  in DMFieldEvaluateFV_DS()
596             const PetscScalar *spH = &sqH[((i * Nq + k) * Nc + j) * dimC * dimC];  in DMFieldEvaluateFV_DS()
606 …for (m = 0; m < dimC; m++) sH[(i * Nc + j) * dimC * dimC + l * dimC + m] += geom->detJ[i * Nq + k]…  in DMFieldEvaluateFV_DS()
610         for (k = 0; k < Nc * dimC * dimC; k++) sH[i * Nc * dimC * dimC + k] /= vol;  in DMFieldEvaluateFV_DS()
619         for (j = 0; j < Nc * dimC * dimC; j++) rH[i * Nc * dimC * dimC + j] = 0.;  in DMFieldEvaluateFV_DS()
624           for (j = 0; j < Nc; j++) {  in DMFieldEvaluateFV_DS()
630             const PetscReal *rpH = &rqH[((i * Nq + k) * Nc + j) * dimC * dimC];  in DMFieldEvaluateFV_DS()
640 …for (m = 0; m < dimC; m++) rH[(i * Nc + j) * dimC * dimC + l * dimC + m] += geom->detJ[i * Nq + k]…  in DMFieldEvaluateFV_DS()
644         for (k = 0; k < Nc * dimC * dimC; k++) rH[i * Nc * dimC * dimC + k] /= vol;  in DMFieldEvaluateFV_DS()
844       PetscInt   tdim, Nc, Np;  in DMFieldComputeFaceData_DS()  local
851         PetscCall(PetscQuadratureGetData(quad, NULL, &Nc, &Np, NULL, NULL));  in DMFieldComputeFaceData_DS()
853         PetscCall(PetscCalloc1(Nc * Np, &weights));  in DMFieldComputeFaceData_DS()
854         PetscCall(PetscQuadratureSetData(cellQuad, dim + 1, Nc, Np, points, weights));  in DMFieldComputeFaceData_DS()