52 FaceListType& faces =
pp;
61 for (
const label celli : static_cast<const labelHashSet&>(
set))
63 const cell& cFaces =
mesh.cells()[celli];
65 for (
const label facei : cFaces)
67 if (
mesh.isInternalFace(facei))
69 label otherCelli =
mesh.faceOwner()[facei];
71 if (otherCelli == celli)
73 otherCelli =
mesh.faceNeighbour()[facei];
76 if (!
set.
found(otherCelli))
78 cellFaces.
insert(facei, celli);
83 cellFaces.
insert(facei, celli);
89 faces.addressing() = cellFaces.
sortedToc();
95 writer.open(file, parallel);
110 const labelList& faceIds = faces.addressing();
113 const label cellIdOffset =
115 writer.parallel() ? globalIndex::calcOffset(
mesh.nCells()) : 0
120 faceValues[facei] = cellFaces[faceIds[facei]] + cellIdOffset;
123 writer.write(
"faceID", faceValues);
Write concrete PrimitivePatch faces/points (optionally with fields) as a vtp file or a legacy vtk fil...
A class for handling file names.
vtk::lineWriter writer(edgeCentres, edgeList::null(), fileName(aMesh.time().globalPath()/"finiteArea-edgesCentres"))
bool writeCellSetFaces(const polyMesh &mesh, const cellSet &set, const vtk::outputOptions opts, const fileName &file, bool parallel=UPstream::parRun())
Write perimeter faces of cellset to vtk polydata file.
#define forAll(list, i)
Loop across all elements in list.
Encapsulated combinations of output format options. This is primarily useful when defining the output...
bool insert(const Key &key, const T &obj)
Copy insert a new entry, not overwriting existing entries.
A list of faces which address into the list of points.
word name(const expressions::valueTypeCode typeCode)
A word representation of a valueTypeCode. Empty for expressions::valueTypeCode::INVALID.
Write topoSet in VTK format.
A cell is defined as a list of faces with extra functionality.
A collection of cell labels.
List< Key > sortedToc() const
The table of contents (the keys) in sorted order.
Mesh consisting of general polyhedral cells.
List< label > labelList
A List of labels.
A List with indirect addressing.
uindirectPrimitivePatch pp(UIndirectList< face >(mesh.faces(), faceLabels), mesh.points())