46 Foam::polyBoundaryMesh::polyBoundaryMesh
65 <<
"Specified IOobject::MUST_READ_IF_MODIFIED but class" 66 <<
" does not support automatic rereading." 74 Istream& is = readStream(typeName);
86 patchEntries[
patchi].keyword(),
97 "polyBoundaryMesh::polyBoundaryMesh" 98 "(const IOobject&, const polyMesh&)" 106 Foam::polyBoundaryMesh::polyBoundaryMesh
119 Foam::polyBoundaryMesh::polyBoundaryMesh
141 <<
"Specified IOobject::MUST_READ_IF_MODIFIED but class" 142 <<
" does not support automatic rereading." 149 Istream& is = readStream(typeName);
161 patchEntries[
patchi].keyword(),
172 "polyBoundaryMesh::polyBoundaryMesh" 173 "(const IOobject&, const polyMesh&, const polyPatchList&)" 200 operator[](
patchi).clearGeom();
207 neighbourEdgesPtr_.clear();
209 groupPatchIDsPtr_.clear();
213 operator[](
patchi).clearAddressing();
220 void Foam::polyBoundaryMesh::calcGeometry()
232 operator[](
patchi).initGeometry(pBufs);
239 operator[](
patchi).calcGeometry(pBufs);
249 forAll(patchSchedule, patchEvali)
251 const label patchi = patchSchedule[patchEvali].patch;
253 if (patchSchedule[patchEvali].init)
255 operator[](patchi).initGeometry(pBufs);
259 operator[](patchi).calcGeometry(pBufs);
272 <<
"Neighbour edge addressing not correct across parallel" 273 <<
" boundaries." <<
endl;
276 if (!neighbourEdgesPtr_.valid())
282 label nEdgePairs = 0;
313 edgei < edges.
size();
318 const edge&
e = edges[edgei];
324 pointsToEdge.
find(meshEdge);
326 if (fnd == pointsToEdge.
end())
348 neighbourEdges[edgeInfo[0]][edgeInfo[1]]
354 pointsToEdge.
erase(meshEdge);
359 if (pointsToEdge.
size())
362 <<
"Not all boundary edges of patches match up." <<
nl 363 <<
"Is the outside of your mesh multiply connected?" 377 if (edgeInfo[0] == -1 || edgeInfo[1] == -1)
383 <<
"Not all boundary edges of patches match up." <<
nl 384 <<
"Edge " << edgeI <<
" on patch " << pp.
name()
386 << pp.
localPoints()[e[1]] <<
" is not matched to an" 387 <<
" edge on any other patch." <<
nl 388 <<
"Is the outside of your mesh multiply connected?" 395 return neighbourEdgesPtr_();
401 if (!patchIDPtr_.valid())
408 - mesh_.nInternalFaces()
417 label bFacei = bm[
patchi].start() - mesh_.nInternalFaces();
420 patchID[bFacei++] =
patchi;
424 return patchIDPtr_();
431 if (!groupPatchIDsPtr_.valid())
451 if (iter != groupPatchIDs.
end())
468 <<
"Removing patchGroup '" << bm[
patchi].
name()
469 <<
"' which clashes with patch " <<
patchi 470 <<
" of the same name." 476 return groupPatchIDsPtr_();
482 const word& groupName,
486 groupPatchIDsPtr_.clear();
519 if (groups[i] != groupName)
521 groups[newI++] = groups[i];
579 const bool usePatchGroups
590 if (usePatchGroups && groupPatchIDs().size())
594 const wordList allGroupNames = groupPatchIDs().toc();
598 const word& grpName = allGroupNames[groupIDs[i]];
599 const labelList& patchIDs = groupPatchIDs()[grpName];
602 if (indexSet.
insert(patchIDs[j]))
604 indices.
append(patchIDs[j]);
618 if (key ==
operator[](i).
name())
625 if (usePatchGroups && groupPatchIDs().size())
628 groupPatchIDs().find(key);
630 if (iter != groupPatchIDs().end())
637 if (indexSet.
insert(patchIDs[j]))
639 indices.
append(patchIDs[j]);
660 if (!indices.
empty())
669 if (key ==
operator[](i).
name())
697 Pout<<
"label polyBoundaryMesh::findPatchID(const word&) const" 698 <<
"Patch named " << patchName <<
" not found. " 699 <<
"List of available patch names: " << names() <<
endl;
713 if (faceIndex <
mesh().nInternalFaces())
717 else if (faceIndex >=
mesh().nFaces())
720 <<
"given label " << faceIndex
721 <<
" greater than the number of geometric faces " <<
mesh().
nFaces()
732 faceIndex >= bp.
start()
733 && faceIndex < bp.
start() + bp.size()
742 <<
"Cannot find face " << faceIndex <<
" in any of the patches " 744 <<
"It seems your patches are not consistent with the mesh :" 746 <<
" total number of faces:" <<
mesh().
nFaces()
756 const bool warnNotFound,
757 const bool usePatchGroups
760 const wordList allPatchNames(this->names());
765 const wordRe& patchName = patchNames[i];
776 if (patchIDs.
empty())
780 const wordList allGroupNames = groupPatchIDs().toc();
787 const word&
name = allGroupNames[groupIDs[i]];
790 forAll(extraPatchIDs, extraI)
792 ids.
insert(extraPatchIDs[extraI]);
796 if (groupIDs.
empty() && warnNotFound)
799 <<
"Cannot find any patch or group names matching " 804 else if (warnNotFound)
807 <<
"Cannot find any patch names matching " << patchName
834 groupPatchIDs.
begin();
835 iter != groupPatchIDs.
end();
846 label nMatch = nonGroupPatches.
erase(groupPatchSet);
848 if (nMatch == groupPatchSet.
size())
850 matchedGroups.
append(iter.key());
852 else if (nMatch != 0)
855 nonGroupPatches.
transfer(oldNonGroupPatches);
873 bool hasError =
false;
883 if (!isA<processorPolyPatch>(bm[
patchi]))
885 if (nonProci != patchi)
892 Pout<<
" ***Problem with boundary patch " << patchi
894 <<
" of type " << bm[
patchi].type()
895 <<
". The patch seems to be preceded by processor" 896 <<
" patches. This is can give problems." 903 types[nonProci] = bm[
patchi].type();
909 types.setSize(nonProci);
923 for (
label proci = 1; proci < allNames.
size(); ++proci)
927 (allNames[proci] != allNames[0])
928 || (allTypes[proci] != allTypes[0])
935 Info<<
" ***Inconsistent patches across processors, " 936 "processor 0 has patch names:" << allNames[0]
937 <<
" patch types:" << allTypes[0]
938 <<
" processor " << proci <<
" has patch names:" 940 <<
" patch types:" << allTypes[proci]
955 bool hasError =
false;
961 if (bm[
patchi].start() != nextPatchStart && !hasError)
965 Info<<
" ****Problem with boundary patch " <<
patchi 967 <<
" of type " << bm[
patchi].type()
968 <<
". The patch should start on face no " << nextPatchStart
969 <<
" and the patch specifies " << bm[
patchi].start()
971 <<
"Possibly consecutive patches have this same problem." 972 <<
" Suppressing future warnings." <<
endl;
979 Info<<
" ****Duplicate boundary patch " <<
patchi 981 <<
" of type " << bm[
patchi].type()
983 <<
"Suppressing future warnings." <<
endl;
995 Pout<<
" ***Boundary definition is in error." <<
endl;
999 Info<<
" Boundary definition OK." <<
endl;
1019 operator[](
patchi).initMovePoints(pBufs, p);
1026 operator[](
patchi).movePoints(pBufs, p);
1036 forAll(patchSchedule, patchEvali)
1038 const label patchi = patchSchedule[patchEvali].patch;
1040 if (patchSchedule[patchEvali].init)
1042 operator[](patchi).initMovePoints(pBufs, p);
1046 operator[](patchi).movePoints(pBufs, p);
1055 neighbourEdgesPtr_.clear();
1056 patchIDPtr_.clear();
1057 groupPatchIDsPtr_.clear();
1069 operator[](
patchi).initUpdateMesh(pBufs);
1076 operator[](
patchi).updateMesh(pBufs);
1086 forAll(patchSchedule, patchEvali)
1088 const label patchi = patchSchedule[patchEvali].patch;
1090 if (patchSchedule[patchEvali].init)
1092 operator[](patchi).initUpdateMesh(pBufs);
1096 operator[](patchi).updateMesh(pBufs);
1106 const bool validBoundary
1144 os.
check(
"polyBoundaryMesh::writeData(Ostream& os) const");
1165 const word& patchName
1173 <<
"Patch named " << patchName <<
" not found." <<
nl 1174 <<
"Available patch names: " << names() <<
endl 1178 return operator[](patchi);
1184 const word& patchName
1192 <<
"Patch named " << patchName <<
" not found." <<
nl 1193 <<
"Available patch names: " << names() <<
endl 1197 return operator[](patchi);
PtrList< polyPatch > polyPatchList
container classes for polyPatch
A class for handling keywords in dictionaries.
static void scatterList(const List< commsStruct > &comms, List< T > &Values, const int tag, const label comm)
Scatter data. Reverse of gatherList.
const List< labelPairList > & neighbourEdges() const
Per patch the edges on the neighbouring patch. Is for every external.
#define forAll(list, i)
Loop across all elements in list.
bool empty() const
Return true if the UList is empty (ie, size() is zero)
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
bool writeData(Ostream &) const
writeData member function required by regIOobject
const word & name() const
Return name.
const word & name() const
Return name.
Ostream & indent(Ostream &os)
Indent stream.
bool set(const label) const
Is element set.
void finishedSends(const bool block=true)
Mark all sends as having been done. This will start receives.
An STL-conforming const_iterator.
static iteratorEnd end()
iteratorEnd set to beyond the end of any HashTable
void reorder(const labelUList &)
Reorders elements. Ordering does not have to be done in.
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
virtual bool check(const char *operation) const
Check IOstream status for given operation.
label nInternalFaces() const
const lduSchedule & patchSchedule() const
Order in which the patches should be initialised/evaluated.
bool inGroup(const word &) const
Test if in group.
An Istream is an abstract base class for all input systems (streams, files, token lists etc)...
void size(const label)
Override size to be inconsistent with allocated storage.
static int myProcNo(const label communicator=0)
Number of this process (starting from masterNo() = 0)
Ostream & endl(Ostream &os)
Add newline and flush stream.
const labelList & patchID() const
Per boundary face label the patch index.
static bool master(const label communicator=0)
Am I the master process.
label findPatchID(const word &patchName) const
Find patch index given a name.
label nInternalEdges() const
Number of internal edges.
void clearGeom()
Clear geometry at this level and at patches.
const wordList & inGroups() const
Return the optional groups patch belongs to.
Operations on lists of strings.
void transfer(HashTable< T, Key, Hash > &)
Transfer the contents of the argument table into this table.
wordList types() const
Return a list of patch types.
labelList findIndices(const ListType &, typename ListType::const_reference, const label start=0)
Find all occurrences of given element. Linear search.
bool insert(const Key &key)
Insert a new entry.
labelHashSet patchSet(const UList< wordRe > &patchNames, const bool warnNotFound=true, const bool usePatchGroups=true) const
Return the set of patch IDs corresponding to the given names.
bool good() const
Return true if next operation might succeed.
label size() const
Return number of elements in table.
const labelList & meshPoints() const
Return labelList of mesh points in patch. They are constructed.
bool erase(const iterator &)
Erase a hashedEntry specified by given iterator.
bool findStrings(const wordReListMatcher &matcher, const std::string &str)
Return true if string matches one of the regular expressions.
friend Ostream & operator(Ostream &, const UPtrList< T > &)
Write UPtrList to Ostream.
virtual bool writeObject(IOstream::streamFormat, IOstream::versionNumber, IOstream::compressionType, const bool valid) const
Write using given format, version and compression.
bool insert(const Key &, const T &newElmt)
Insert a new hashedEntry.
iterator find(const Key &)
Find and return an iterator set at the hashedEntry.
HashSet< label, Hash< label > > labelHashSet
A HashSet with label keys.
An STL-conforming iterator.
void matchGroups(const labelUList &patchIDs, wordList &groups, labelHashSet &nonGroupPatches) const
Match the patches to groups. Returns all the (fully matched) groups.
An ordered pair of two objects of type <T> with first() and second() elements.
~polyBoundaryMesh()
Destructor.
void movePoints(const pointField &)
Correct polyBoundaryMesh after moving points.
An edge is a list of two point labels. The functionality it provides supports the discretisation on a...
A class for handling words, derived from string.
bool writeObject(IOstream::streamFormat fmt, IOstream::versionNumber ver, IOstream::compressionType cmp, const bool valid) const
Write using given format, version and form uncompression.
bool checkDefinition(const bool report=false) const
Check boundary definition. Return true if in error.
wordList names() const
Return a list of patch names.
void append(const T &)
Append an element at the end of the list.
void setCapacity(const label)
Alter the size of the underlying storage.
wordList patchNames(nPatches)
streamFormat
Enumeration for the format of data in the stream.
DynamicList< T, SizeInc, SizeMult, SizeDiv > & append(const T &)
Append an element at the end of the list.
Pair< label > labelPair
Label pair.
const globalMeshData & globalData() const
Return parallel info.
const edgeList & edges() const
Return list of edges, address into LOCAL point list.
List< label > labelList
A List of labels.
wordList physicalTypes() const
Return a list of physical types.
iterator begin()
Iterator set to the beginning of the HashTable.
An STL-conforming hash table.
errorManip< error > abort(error &err)
A wordRe is a word, but can also have a regular expression for matching words.
A 1D vector of objects of type <T>, where the size of the vector is known and can be used for subscri...
void setSize(const label)
Reset size of PtrList. If extending the PtrList, new entries are.
compressionType
Enumeration for the format of data in the stream.
An Ostream is an abstract base class for all output systems (streams, files, token lists...
const Field< PointType > & localPoints() const
Return pointField of points in patch.
bool checkParallelSync(const bool report=false) const
Check whether all procs have all patches and in same order. Return.
label nEdges() const
Return number of edges in patch.
void updateMesh()
Correct polyBoundaryMesh after topology update.
defineTypeNameAndDebug(combustionModel, 0)
Ostream & decrIndent(Ostream &os)
Decrement the indent level.
Buffers for inter-processor communications streams (UOPstream, UIPstream).
void reduce(const List< UPstream::commsStruct > &comms, T &Value, const BinaryOp &bop, const int tag, const label comm)
word name(const complex &)
Return a string representation of a complex.
label size() const
Return the number of elements in the UPtrList.
static commsTypes defaultCommsType
Default commsType.
void setSize(const label)
Reset size of List.
static bool & parRun()
Is this a parallel run?
static label nProcs(const label communicator=0)
Number of processes in parallel run.
#define WarningInFunction
Report a warning using Foam::Warning.
A templated 1D list of pointers to objects of type <T>, where the size of the array is known and used...
void clearAddressing()
Clear addressing at this level and at patches.
prefixOSstream Pout(cout, "Pout")
label start() const
Return start label of this patch in the polyMesh face list.
Ostream & operator<<(Ostream &, const ensightPart &)
labelList findIndices(const keyType &, const bool usePatchGroups=true) const
Return patch indices for all matches. Optionally matches patchGroups.
regIOobject is an abstract class derived from IOobject to handle automatic object registration with t...
static autoPtr< polyPatch > New(const word &patchType, const word &name, const label size, const label start, const label index, const polyBoundaryMesh &bm)
Return a pointer to a new patch created on freestore from.
const doubleScalar e
Elementary charge.
Mesh consisting of general polyhedral cells.
A patch is a list of labels that address the faces in the global face list.
Ostream & incrIndent(Ostream &os)
Increment the indent level.
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
static void gatherList(const List< commsStruct > &comms, List< T > &Values, const int tag, const label comm)
Gather data but keep individual values separate.
label whichPatch(const label faceIndex) const
Return patch index for a given face label.
label findIndex(const keyType &) const
Return patch index for the first match, return -1 if not found.
void transfer(List< T > &)
Transfer the contents of the argument List into this list.
void reorder(const labelUList &, const bool validBoundary)
Reorders patches. Ordering does not have to be done in.
void setGroup(const word &groupName, const labelList &patchIDs)
Set/add group with patches.
bool isPattern() const
Should be treated as a match rather than a literal string.
const HashTable< labelList, word > & groupPatchIDs() const
Per patch group the patch indices.