87 <<
"Cannot find patch " <<
name
88 <<
" in the source mesh.\n"
89 <<
"Valid patch names are " <<
patches.names()
157 faceLabels[
n] = pp.
start()+j;
178 faceLabels[i] = pp.
start() + i;
194 label oldFacei = faceLabels[i];
196 if (reverseMap[oldFacei] >= 0)
198 newFaceLabels[newI++] = reverseMap[oldFacei];
214 label oldCelli = cellLabels[i];
216 if (reverseMap[oldCelli] >= 0)
218 newCellLabels.
insert(reverseMap[oldCelli]);
225 template<
class PatchType>
226 void changeFrontBackPatches
229 const word& frontPatchName,
230 const word& backPatchName
271 int main(
int argc,
char *argv[])
288 const ExtrudeMode
mode = ExtrudeModeNames.
read
296 Info<<
"Extruding from " << ExtrudeModeNames[
mode]
297 <<
" using model " << model().type() <<
endl;
300 Info<<
"Flipping normals before extruding" <<
endl;
304 Info<<
"Collapsing edges < " << mergeTol <<
" of bounding box" <<
endl;
308 Info<<
"Not collapsing any edges after extrusion" <<
endl;
328 if (flipNormals &&
mode == MESH)
331 <<
"Flipping normals not supported for extrusions from mesh."
336 sourceCasePath.expand();
347 if (sourcePatches.size() == 1)
349 frontPatchName = sourcePatches[0];
353 if (zoneNames.
size() == 1)
355 if (zoneNames[0] ==
"patchNames")
357 zoneNames = sourcePatches;
361 zoneNames =
wordList(sourcePatches.size(), zoneNames[0]);
364 else if (zoneNames.
size() && zoneNames.
size() != sourcePatches.size())
367 <<
"Number of zoneNames "
368 "does not equal the number of sourcePatches"
372 Info<<
"Extruding patches " << sourcePatches
373 <<
" on mesh " << sourceCasePath <<
nl
390 ? findPatchZones(mesh, zoneNames)
411 if (
mode == PATCH && flipNormals)
418 label patchFacei = patchFaces[i];
430 mesh.
faces()[patchFacei].reverseFace(),
467 extrudePatch.meshEdges
522 Info<<
"Adding overall " << nAdded <<
" processor patches." <<
endl;
548 <<
" and " << nbrProci
566 mesh.removeFvBoundary();
567 mesh.addFvPatches(newPatches,
true);
576 dict.
lookup(
"exposedPatchName") >> backPatchName;
585 pointField layer0Points(extrudePatch.nPoints());
586 pointField displacement(extrudePatch.nPoints());
587 forAll(displacement, pointi)
589 const vector& patchNormal = extrudePatchPointNormals[pointi];
592 layer0Points[pointi] = model()
594 extrudePatch.localPoints()[pointi],
599 point extrudePt = model()
601 extrudePatch.localPoints()[pointi],
605 displacement[pointi] = extrudePt - layer0Points[pointi];
611 if (
gMax(
mag(layer0Points-extrudePatch.localPoints())) > small)
613 Info<<
"Moving mesh to layer0 points since differ from original"
614 <<
" points - this can happen for wedge extrusions." <<
nl
618 forAll(extrudePatch.meshPoints(), i)
620 newPoints[extrudePatch.meshPoints()[i]] = layer0Points[i];
627 labelList nFaceLayers(extrudePatch.size(), model().nLayers());
630 labelList nPointLayers(extrudePatch.nPoints(), model().nLayers());
633 vectorField firstLayerDisp(displacement*model().sumThickness(1));
649 layerExtrude.setRefinement
666 forAll(layerExtrude.addedPoints(), pointi)
668 const labelList& pPoints = layerExtrude.addedPoints()[pointi];
671 label meshPointi = pPoints[pPointi];
677 extrudePatch.localPoints()[pointi],
678 extrudePatchPointNormals[pointi],
686 )[meshPointi] = modelPt;
695 forAll(backPatchFaces, patchFacei)
697 backPatchFaces[patchFacei] = layerFaces[patchFacei].
first();
698 frontPatchFaces[patchFacei] = layerFaces[patchFacei].
last();
709 runTimeExtruded.constant(),
718 layerExtrude.updateZones(meshFromMesh());
720 meshFromMesh().topoChange(map);
722 layerExtrude.topoChange
732 map().reverseFaceMap(),
737 map().reverseFaceMap(),
746 layerExtrude.addedCells
749 layerExtrude.layerFaces()
754 const labelList& aCells = addedCells[facei];
757 addedCellsSet.
insert(aCells[i]);
768 Info<<
"Extruding surfaceMesh read from file " << surfName <<
nl
779 faces[i] = fMesh[i].reverseFace();
783 Info<<
"Extruding surface with :" <<
nl
784 <<
" points : " << fMesh.points().size() <<
nl
785 <<
" faces : " << fMesh.size() <<
nl
786 <<
" normals[0] : " << fMesh.faceNormals()[0]
790 meshFromSurface.
reset
797 runTimeExtruded.constant(),
807 frontPatchName =
"originalPatch";
808 frontPatchFaces = findPatchFaces
810 meshFromSurface().boundaryMesh(),
813 backPatchName =
"otherSide";
814 backPatchFaces = findPatchFaces
816 meshFromSurface().boundaryMesh(),
832 const scalar mergeDim = mergeTol * bb.
minDim();
834 Info<<
"Mesh bounding box : " << bb <<
nl
835 <<
" with span : " << span <<
nl
836 <<
"Merge distance : " << mergeDim <<
nl
845 Pout<<
"Collapsing edges < " << mergeDim <<
" ..." <<
nl <<
endl;
858 const edge&
e = edges[edgeI];
864 Pout<<
"Merging edge " <<
e <<
" since length " << d
865 <<
" << " << mergeDim <<
nl;
868 collapsePointToLocation.set(
e[1],
points[
e[0]]);
876 collapser.consistentCollapse
880 collapsePointToLocation,
889 bool anyChange = collapser.setRefinement(allPointInfo, meshMod);
901 updateFaceLabels(map(), frontPatchFaces);
902 updateFaceLabels(map(), backPatchFaces);
903 updateCellSet(map(), addedCellsSet);
913 if (isType<extrudeModels::wedge>(model()))
915 changeFrontBackPatches<wedgePolyPatch>
922 else if (isType<extrudeModels::plane>(model()))
924 changeFrontBackPatches<emptyPolyPatch>
942 <<
"Cannot stitch front and back of extrusion since"
943 <<
" in 'mesh' mode (extrusion appended to mesh)."
947 Info<<
"Assuming full 360 degree axisymmetric case;"
948 <<
" stitching faces on patches "
949 << frontPatchName <<
" and "
950 << backPatchName <<
" together ..." <<
nl <<
endl;
952 if (frontPatchFaces.
size() != backPatchFaces.
size())
955 <<
"Differing number of faces on front ("
956 << frontPatchFaces.
size() <<
") and back ("
957 << backPatchFaces.
size() <<
")"
973 perfectStitcher.setRefinement
1003 updateCellSet(map(), addedCellsSet);
1019 cellSet addedCells(mesh,
"addedCells", addedCellsSet);
1020 Info<<
"Writing added cells to cellSet " << addedCells.name()
1022 if (!addedCells.write())
1025 << addedCells.name()
#define forAll(list, i)
Loop across all elements in list.
A 1D vector of objects of type <T> that resizes itself as necessary to accept the new objects.
bool insert(const Key &key)
Insert a new entry.
void transfer(HashTable< T, Key, Hash > &)
Transfer the contents of the argument table into this table.
label size() const
Return number of elements in table.
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
fileName relativeObjectPath() const
Return complete relativePath + object name.
autoPtr< IOobject > clone() const
Clone.
A List with indirect addressing.
void transfer(List< T > &)
Transfer the contents of the argument List into this list.
void size(const label)
Override size to be inconsistent with allocated storage.
static const List< word > & null()
Return a null List.
void setSize(const label)
Reset size of List.
Initialise the NamedEnum HashTable from the static list of names.
Enum read(Istream &) const
Read a word from Istream and return the corresponding.
autoPtr< T > set(const label, const word &key, T *)
Set element to pointer provided and return old element.
label findIndex(const word &key) const
Return the index of the given the key or -1 if not found.
void setSize(const label)
Reset size of PtrList. If extending the PtrList, new entries are.
A simple wrapper around bool so that it can be read as a word: true/false, on/off,...
Class to control time during OpenFOAM simulations that is also the top-level objectRegistry.
static word controlDictName
The default control dictionary name (normally "controlDict")
T & first()
Return the first element of the list.
label size() const
Return the number of elements in the UList.
T & last()
Return the last element of the list.
static bool & parRun()
Is this a parallel run?
static int myProcNo(const label communicator=0)
Number of this process (starting from masterNo() = 0)
label size() const
Return the number of elements in the UPtrList.
Adds layers of cells to outside of polyPatch. Can optionally create stand-alone extruded mesh (addToM...
static void calcSidePatch(const polyMesh &, const globalIndex &globalFaces, const labelListList &globalEdgeFaces, const indirectPrimitivePatch &pp, labelList &sidePatchID, label &nPatches, Map< label > &nbrProcToPatch, Map< label > &patchToNbrProc)
Boundary edges get extruded into boundary faces. Determine patch.
static labelListList globalEdgeFaces(const polyMesh &, const globalIndex &globalFaces, const indirectPrimitivePatch &pp)
Per patch edge the pp faces (in global indices) using it. Uses.
An auto-pointer similar to the STL auto_ptr but with automatic casting to a reference to the type and...
bool valid() const
Return true if the autoPtr valid (ie, the pointer is set)
void reset(T *=nullptr)
If object pointer already set, delete object and set to given.
A bounding box defined in terms of the points at its extremities.
scalar minDim() const
Smallest length/height/width dimension.
vector span() const
The bounding box span (from minimum to maximum)
A collection of cell labels.
A list of keyword definitions, which are a keyword followed by any number of values (e....
T lookupOrDefault(const word &, const T &, bool recursive=false, bool patternMatch=true) const
Find and return a T, if not found return the given default.
ITstream & lookup(const word &, bool recursive=false, bool patternMatch=true) const
Find and return an entry data stream.
Does polyTopoChanges to remove edges. Can remove faces due to edge collapse but can not remove cells ...
An edge is a list of two point labels. The functionality it provides supports the discretisation on a...
static autoPtr< extrudeModel > New(const dictionary &)
Select null constructed.
A class for handling file names.
word name() const
Return file name (part beyond last /)
fileName path() const
Return directory path name (part before last /)
Calculates a unique integer (label so might not have enough room - 2G max) for processor + local inde...
Calculates points shared by more than two processor patches or cyclic patches.
label index() const
Return the index of this patch in the boundaryMesh.
const word & name() const
Return name.
Calculates zone number for every face of patch.
Hack of attachDetach to couple patches when they perfectly align. Does not decouple....
Mesh consisting of general polyhedral cells.
const cellZoneList & cellZones() const
Return cell zones.
static word defaultRegion
Return the default region name.
virtual const faceList & faces() const
Return raw faces.
virtual const labelList & faceOwner() const
Return face owner.
virtual void topoChange(const polyTopoChangeMap &)
Update topology using the given map.
const polyBoundaryMesh & boundaryMesh() const
Return boundary mesh.
virtual const labelList & faceNeighbour() const
Return face neighbour.
void addPatches(const List< polyPatch * > &, const bool validBoundary=true)
Add boundary patches.
virtual const pointField & points() const
Return raw points.
void removeBoundary()
Remove boundary patches.
void setInstance(const fileName &)
Set the instance for mesh files.
const boundBox & bounds() const
Return mesh bounding box.
virtual void setPoints(const pointField &)
Reset the points.
A patch is a list of labels that address the faces in the global face list.
label start() const
Return start label of this patch in the polyMesh face list.
virtual autoPtr< polyPatch > clone(const polyBoundaryMesh &bm) const
Construct and return a clone, resetting the boundary mesh.
Class containing mesh-to-mesh mapping information after a change in polyMesh topology.
const labelList & reverseCellMap() const
Reverse cell map.
const labelList & reverseFaceMap() const
Reverse face map.
Direct mesh changes based on v1.3 polyTopoChange syntax.
const labelListList & pointEdges() const
const edgeList & edges() const
Return mesh edges. Uses calcEdges.
Neighbour processor patch.
virtual bool write(const bool write=true) const
Write using setting from DB.
A class for handling words, derived from string.
static const word null
An empty word.
label collapseEdge(triSurface &surf, const scalar minLen)
Keep collapsing all edges < minLen.
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
int main(int argc, char *argv[])
const fvPatchList & patches
errorManipArg< error, int > exit(error &err, const int errNo=1)
List< word > wordList
A List of words.
mode_t mode(const fileName &, const bool checkVariants=true, const bool followLink=true)
Return the file mode.
List< label > labelList
A List of labels.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
IOdictionary systemDict(const word &dictName, const argList &args, const objectRegistry &ob, const word ®ionName=polyMesh::defaultRegion)
Ostream & endl(Ostream &os)
Add newline and flush stream.
word name(const bool)
Return a word representation of a bool.
const word & regionName(const solver ®ion)
ListType renumber(const labelUList &oldToNew, const ListType &)
Renumber the values (not the indices) of a list.
void reduce(const List< UPstream::commsStruct > &comms, T &Value, const BinaryOp &bop, const int tag, const label comm)
T returnReduce(const T &Value, const BinaryOp &bop, const int tag=Pstream::msgType(), const label comm=UPstream::worldComm)
dimensioned< scalar > mag(const dimensioned< Type > &)
prefixOSstream Pout(cout, "Pout")
label findIndex(const ListType &, typename ListType::const_reference, const label start=0)
Find first occurrence of given element and return index,.
labelList identityMap(const label len)
Create identity map (map[i] == i) of given length.
Type gMax(const FieldField< Field, Type > &f)
wordList patchNames(nPatches)
Foam::argList args(argc, argv)