16 static PetscErrorCode CreateHalfCellsLabel(DM dm, PetscBool lower, DMLabel *label)  in CreateHalfCellsLabel()  argument
22   PetscCall(DMGetLabel(dm, "cells", label));  in CreateHalfCellsLabel()
23   PetscCall(DMLabelClearStratum(*label, 1));  in CreateHalfCellsLabel()
32   for (PetscInt c = cStartSub; c < cEndSub; ++c) PetscCall(DMLabelSetValue(*label, c, 1));  in CreateHalfCellsLabel()
33   PetscCall(DMPlexLabelComplete(dm, *label));  in CreateHalfCellsLabel()
38 …atic PetscErrorCode CreateHalfDomainLabel(DM dm, PetscBool lower, PetscReal height, DMLabel *label)  in CreateHalfDomainLabel()  argument
46   PetscCall(DMGetLabel(dm, "cells", label));  in CreateHalfDomainLabel()
47   PetscCall(DMLabelClearStratum(*label, 1));  in CreateHalfDomainLabel()
54       if (centroid[0] < 0.5) PetscCall(DMLabelSetValue(*label, c, 1));  in CreateHalfDomainLabel()
56       if (centroid[0] > 0.5) PetscCall(DMLabelSetValue(*label, c, 1));  in CreateHalfDomainLabel()
59   PetscCall(DMPlexLabelComplete(dm, *label));  in CreateHalfDomainLabel()
64 static PetscErrorCode CreateLineLabel(DM dm, PetscReal x, DMLabel *label)  in CreateLineLabel()  argument
72   PetscCall(DMGetLabel(dm, "faces", label));  in CreateLineLabel()
76     if (PetscAbsReal(centroid[0] - x) < PETSC_SMALL) PetscCall(DMLabelSetValue(*label, f, 1));  in CreateLineLabel()
78   PetscCall(DMPlexLabelComplete(dm, *label));  in CreateLineLabel()
84   DMLabel label, map;  in CreateVolumeSubmesh()  local
87   if (domain) PetscCall(CreateHalfDomainLabel(dm, lower, height, &label));  in CreateVolumeSubmesh()
88   else PetscCall(CreateHalfCellsLabel(dm, lower, &label));  in CreateVolumeSubmesh()
89 …PetscCall(DMPlexFilter(dm, label, 1, PETSC_FALSE, PETSC_FALSE, PetscObjectComm((PetscObject)dm), N…  in CreateVolumeSubmesh()
101   DMLabel  label, map;  in TestBoundaryField()  local
107   PetscCall(CreateLineLabel(dm, 0.5, &label));  in TestBoundaryField()
108 …PetscCall(DMPlexFilter(dm, label, 1, PETSC_FALSE, PETSC_FALSE, PetscObjectComm((PetscObject)dm), N…  in TestBoundaryField()