7   PetscInt   depth, Nf, f, pStart, pEnd;  in DMPlexCreateSectionFields()  local
15   for (f = 0; f < Nf; ++f) {  in DMPlexCreateSectionFields()
19     PetscCall(DMGetField(dm, f, NULL, &obj));  in DMPlexCreateSectionFields()
22       isFE[f] = PETSC_TRUE;  in DMPlexCreateSectionFields()
24       isFE[f] = PETSC_FALSE;  in DMPlexCreateSectionFields()
32       for (f = 0; f < Nf; ++f) {  in DMPlexCreateSectionFields()
33         PetscCall(PetscSectionSetFieldComponents(*section, f, numComp[f]));  in DMPlexCreateSectionFields()
34         if (isFE[f]) {  in DMPlexCreateSectionFields()
41           PetscCall(DMGetField(dm, f, NULL, (PetscObject *)&fe));  in DMPlexCreateSectionFields()
75             PetscCall(PetscSectionSetFieldSym(*section, f, sym));  in DMPlexCreateSectionFields()
94   PetscInt       Nf, f, Nds, n, dim, d, dep, p;  in DMPlexCreateSectionDof()  local
115   for (f = 0; f < Nf; ++f) {  in DMPlexCreateSectionDof()
119     PetscCall(DMGetField(dm, f, NULL, &obj));  in DMPlexCreateSectionDof()
122     isFE[f] = id == PETSCFE_CLASSID ? PETSC_TRUE : PETSC_FALSE;  in DMPlexCreateSectionDof()
126   for (f = 0; f < Nf; ++f) {  in DMPlexCreateSectionDof()
129     PetscCall(DMGetFieldAvoidTensor(dm, f, &avoidTensor));  in DMPlexCreateSectionDof()
131     if (label && label[f]) {  in DMPlexCreateSectionDof()
136       PetscCall(DMLabelGetStratumIS(label[f], 1, &pointIS));  in DMPlexCreateSectionDof()
157         dof = d < 0 ? 0 : numDof[f * (dim + 1) + d];  in DMPlexCreateSectionDof()
158         PetscCall(PetscSectionSetFieldDof(section, point, f, dof));  in DMPlexCreateSectionDof()
169           const PetscInt dof = numDof[f * (dim + 1) + d];  in DMPlexCreateSectionDof()
177             if (avoidTensor && isFE[f]) continue;  in DMPlexCreateSectionDof()
181           PetscCall(PetscSectionSetFieldDof(section, p, f, dof));  in DMPlexCreateSectionDof()
250       PetscInt dof, f;  in DMPlexCreateSectionBCDof()  local
257         for (f = 0; f < Nf; f++) {  in DMPlexCreateSectionBCDof()
258           PetscCall(PetscSectionGetFieldDof(section, a, f, &dof));  in DMPlexCreateSectionBCDof()
259           PetscCall(PetscSectionSetFieldConstraintDof(section, a, f, dof));  in DMPlexCreateSectionBCDof()
274   PetscInt     Nf, cdof, maxDof = 0, pStart, pEnd, p, bc, f, d;  in DMPlexCreateSectionBCIndicesField()  local
288 … for (f = 0; f < Nf; ++f) PetscCall(PetscSectionSetFieldConstraintIndices(section, p, f, indices));  in DMPlexCreateSectionBCIndicesField()
342       PetscInt dof, f;  in DMPlexCreateSectionBCIndicesField()  local
347 … for (f = 0; f < Nf; f++) PetscCall(PetscSectionSetFieldConstraintIndices(section, a, f, indices));  in DMPlexCreateSectionBCIndicesField()
359   PetscInt  Nf, maxDof, pStart, pEnd, p, f, d;  in DMPlexCreateSectionBCIndices()  local
375         for (f = 0; f < Nf; ++f) {  in DMPlexCreateSectionBCIndices()
379           PetscCall(PetscSectionGetFieldDof(section, p, f, &fdof));  in DMPlexCreateSectionBCIndices()
380           PetscCall(PetscSectionGetFieldConstraintDof(section, p, f, &fcdof));  in DMPlexCreateSectionBCIndices()
382           PetscCall(PetscSectionGetFieldConstraintIndices(section, p, f, &find));  in DMPlexCreateSectionBCIndices()
455   PetscInt     depth, dim, numBC = 0, Nf, Nds, s, bc = 0, f;  in DMCreateLocalSection_Plex()  local
464   for (f = 0; f < Nf; ++f) {  in DMCreateLocalSection_Plex()
468     PetscCall(DMGetField(dm, f, NULL, &obj));  in DMCreateLocalSection_Plex()
471       isFE[f] = PETSC_TRUE;  in DMCreateLocalSection_Plex()
473       isFE[f] = PETSC_FALSE;  in DMCreateLocalSection_Plex()
474 …((PetscObject)dm), PETSC_ERR_ARG_WRONG, "Unknown discretization type for field %" PetscInt_FMT, f);  in DMCreateLocalSection_Plex()
496   for (f = 0; f < Nf; ++f)  in DMCreateLocalSection_Plex()
497     if (!isFE[f] && cEndInterior >= 0) ++numBC;  in DMCreateLocalSection_Plex()
500   for (f = 0; f < Nf; ++f) {  in DMCreateLocalSection_Plex()
503     if (isFE[f] || cEndInterior < 0) continue;  in DMCreateLocalSection_Plex()
506     bcFields[bc] = f;  in DMCreateLocalSection_Plex()
588   for (f = 0; f < Nf; ++f) {  in DMCreateLocalSection_Plex()
589     labels[f] = dm->fields[f].label;  in DMCreateLocalSection_Plex()
590     if (isFE[f]) {  in DMCreateLocalSection_Plex()
591       PetscFE         fe = (PetscFE)dm->fields[f].disc;  in DMCreateLocalSection_Plex()
595       PetscCall(PetscFEGetNumComponents(fe, &numComp[f]));  in DMCreateLocalSection_Plex()
598       for (d = 0; d < PetscMin(dim, fedim) + 1; ++d) numDof[f * (dim + 1) + d] = numFieldDof[d];  in DMCreateLocalSection_Plex()
600       PetscFV fv = (PetscFV)dm->fields[f].disc;  in DMCreateLocalSection_Plex()
602       PetscCall(PetscFVGetNumComponents(fv, &numComp[f]));  in DMCreateLocalSection_Plex()
603       numDof[f * (dim + 1) + dim] = numComp[f];  in DMCreateLocalSection_Plex()
607   for (f = 0; f < Nf; ++f) {  in DMCreateLocalSection_Plex()
609 …for (d = 1; d < dim; ++d) PetscCheck(numDof[f * (dim + 1) + d] <= 0 || depth >= dim, PetscObjectCo…  in DMCreateLocalSection_Plex()
615   for (f = 0; f < Nf; ++f) {  in DMCreateLocalSection_Plex()
619     PetscCall(DMGetField(dm, f, NULL, (PetscObject *)&fe));  in DMCreateLocalSection_Plex()
622     PetscCall(PetscSectionSetFieldName(section, f, name));  in DMCreateLocalSection_Plex()