47 #ifndef Foam_meshRefinement_H 48 #define Foam_meshRefinement_H 71 class mapDistributePolyMesh;
73 class decompositionMethod;
74 class refinementSurfaces;
75 class refinementFeatures;
78 class fvMeshDistribute;
80 class localPointRegion;
161 const scalar mergeDistance_;
164 const bool overwrite_;
167 const word oldInstance_;
218 template<
class GeoField>
219 static void addPatchFields(
fvMesh&,
const word& patchFieldType);
222 template<
class GeoField>
235 label globalFaceCount(
const labelList& elems)
const;
242 void calcCellCellRays
261 const label exposedPatchi
270 static bool markForRefine
272 const label markValue,
273 const label nAllowRefine,
280 void markFeatureCellLevel
288 label markFeatureRefinement
291 const label nAllowRefine,
298 label markInternalDistanceToFeatureRefinement
300 const label nAllowRefine,
306 label markInternalRefinement
308 const label nAllowRefine,
315 label unmarkInternalRefinement
329 label markSurfaceRefinement
331 const label nAllowRefine,
340 void collectGapCandidates
353 const scalar planarCos,
366 const label nAllowRefine,
377 label markSurfaceGapRefinement
379 const scalar planarCos,
380 const label nAllowRefine,
391 const point& nearPoint,
406 const bool useSurfaceNormal,
408 const point& nearPoint,
427 void selectGapCandidates
451 label markInternalGapRefinement
453 const scalar planarCos,
454 const bool spreadGapSize,
455 const label nAllowRefine,
463 label markSmallFeatureRefinement
465 const scalar planarCos,
466 const label nAllowRefine,
477 label markSurfaceFieldRefinement
479 const label nAllowRefine,
492 const scalar tol = 1
e-6
499 const scalar minCosAngle,
500 const scalar lengthScale,
510 label markSurfaceCurvatureRefinement
512 const scalar curvature,
513 const label nAllowRefine,
523 const scalar planarCos,
524 const label nAllowRefine,
526 const label surfaceLevel,
528 const vector& surfaceNormal,
541 label markProximityRefinement
543 const scalar curvature,
549 const label nAllowRefine,
558 label markProximityRefinementWave
560 const scalar planarCos,
562 const label nAllowRefine,
576 const bool allowBoundary,
583 void getIntersections
596 void getIntersections
607 void getBafflePatches
609 const label nErodeCellZones,
614 const bool exitIfLeakPath,
625 const label nBufferLayers,
639 const label ownPatch,
640 const label neiPatch,
659 void markBoundaryFace
686 const scalar minFaceArea,
687 const scalar maxNonOrtho,
694 const scalar volFraction,
701 void markFacesOnProblemCells
704 const bool removeEdgeConnectedCells,
717 const bitSet& isBlockedFace,
745 const label defaultRegion
750 void markFacesOnProblemCellsGeometric
769 const scalar freeStandingAngle
779 void findCellZoneGeometric
790 void findCellZoneInsideWalk
800 void findCellZoneInsideWalk
809 bool calcRegionToZone
811 const label backgroundZoneID,
812 const label surfZoneI,
813 const label ownRegion,
814 const label neiRegion,
823 void findCellZoneTopo
825 const label backgroundZoneID,
837 const label nErodeCellZones,
838 const label backgroundZoneID,
848 const label nGrowCellZones,
849 const label backgroundZoneID,
858 void makeConsistentFaceIndex
868 const bool allowFreeStandingZoneFaces,
869 const label nErodeCellZones,
870 const label backgroundZoneID,
874 const bool exitIfLeakPath,
887 const bitSet& isMasterFace,
891 const bitSet& meshFlipMap,
896 void allocateInterRegionFaceZone
905 void handleSnapProblems
908 const bool useTopologicalSnapDetection,
909 const bool removeEdgeConnectedCells,
931 void calcPatchNumMasterFaces
933 const bitSet& isMasterFace,
948 void consistentOrientation
950 const bitSet& isMasterFace,
1004 return mergeDistance_;
1016 return oldInstance_;
1020 const refinementSurfaces&
surfaces()
const 1026 const refinementFeatures&
features()
const 1032 const shellSurfaces&
shells()
const 1040 return limitShells_;
1058 return faceToCoupledPatch_;
1065 const List<Tuple2<mapType, labelList>>&
userFaceData()
const 1067 return userFaceData_;
1072 return userFaceData_;
1085 autoPtr<mapDistributePolyMesh>
balance 1087 const bool keepZoneFaces,
1088 const bool keepBaffles,
1090 decompositionMethod& decomposer,
1091 fvMeshDistribute& distributor
1101 static autoPtr<indirectPrimitivePatch>
makePatch 1115 const pointMesh& pMesh,
1125 const polyMesh&
mesh,
1126 const bitSet& isMasterEdge,
1135 template<
class Type>
1138 const polyMesh&
mesh,
1139 const bitSet& isMasterEdge,
1143 const Field<Type>& data,
1163 const scalar planarCos,
1176 const scalar curvature,
1177 const scalar planarAngle,
1179 const bool featureRefinement,
1180 const bool featureDistanceRefinement,
1181 const bool internalRefinement,
1182 const bool surfaceRefinement,
1183 const bool curvatureRefinement,
1184 const bool smallFeatureRefinement,
1185 const bool gapRefinement,
1186 const bool bigGapRefinement,
1187 const bool spreadGapSize,
1188 const label maxGlobalCells,
1189 const label maxLocalCells
1208 const bitSet& isOutsideFace,
1216 const bitSet& isOutsideFace,
1224 const scalar planarAngle,
1227 const label growIter
1281 const scalar maxLoadUnbalance,
1282 const label maxCellUnbalance
1292 const scalar maxLoadUnbalance,
1293 const label maxCellUnbalance
1303 const scalar maxLoadUnbalance,
1304 const label maxCellUnbalance
1310 const label maxGlobalCells,
1311 const label maxLocalCells,
1330 const bool handleSnapProblems,
1334 const bool useTopologicalSnapDetection,
1335 const bool removeEdgeConnectedCells,
1337 const label nErodeCellZones,
1345 const bool exitIfLeakPath,
1353 const bool useTopologicalSnapDetection,
1354 const bool removeEdgeConnectedCells,
1356 const scalar planarAngle,
1369 const label nBufferLayers,
1370 const label nErodeCellZones,
1377 const bool exitIfLeakPath,
1384 const label nBufferLayers,
1385 const label nErodeCellZones,
1464 const bool doInternalZones,
1465 const bool doBaffleZones
1475 const bool allowFreeStandingZoneFaces,
1476 const label nErodeCellZones,
1480 const bool exitIfLeakPath,
1492 const label insertPatchi,
1512 const word& masterPatch,
1513 const word& slavePatch,
1522 label& masterPatchID,
1523 label& slavePatchID,
1540 void nearestIntersection
1567 const vector& perturbVec,
1576 const vector& perturbVec,
1585 const vector& perturbVec,
1588 const label nRegions,
1592 const bool exitIfLeakPath,
1605 const bool exitIfLeakPath,
1678 const scalar minCos,
1679 const scalar concaveCos,
1680 const label mergeSize,
1689 const scalar minCos,
1690 const scalar concaveCos,
1731 const scalar minCos,
1823 const bitSet& isMasterElem,
1837 template<
class EnumContainer>
1840 const EnumContainer& namedEnum,
1845 template<
class Type>
1849 const word& keyword,
1852 const Type& deflt =
Zero 1859 const word& keyword,
1868 const word& keyword,
void printMeshInfo(const bool, const string &) const
Print some mesh stats.
label addFaceZone(const word &fzName, const word &masterPatch, const word &slavePatch, const surfaceZonesInfo::faceZoneType &fzType)
Add/lookup faceZone and update information. Return index of.
void markOutsideFaces(const labelList &cellLevel, const labelList &neiLevel, const labelList &refineCell, bitSet &isOutsideFace) const
Mark faces on interface between set and rest.
const refinementSurfaces & surfaces() const
Reference to surface search engines.
const labelIOList & zoneIDs
const labelList patchIDs(pbm.patchSet(polyPatchNames, false, true).sortedToc())
static const dictionary & subDict(const dictionary &dict, const word &keyword, const bool noExit, enum keyType::option matchOpt=keyType::REGEX)
Wrapper around dictionary::subDict which does not exit.
dimensioned< Type > sum(const DimensionedField< Type, GeoMesh > &f1)
A class for handling file names.
const List< Tuple2< mapType, labelList > > & userFaceData() const
Additional face data that is maintained across.
label addPointZone(const word &name)
Add pointZone if does not exist. Return index of zone.
autoPtr< mapPolyMesh > mergePoints(const labelList &pointToDuplicate)
Merge duplicate points.
autoPtr< mapDistributePolyMesh > refineAndBalance(const string &msg, decompositionMethod &decomposer, fvMeshDistribute &distributor, const labelList &cellsToRefine, const scalar maxLoadUnbalance, const label maxCellUnbalance)
Refine some cells and rebalance.
label mergePatchFacesUndo(const scalar minCos, const scalar concaveCos, const labelList &patchIDs, const dictionary &motionDict, const labelList &preserveFaces, const meshRefinement::FaceMergeType mergeType)
Merge coplanar faces. preserveFaces is != -1 for faces.
const word & oldInstance() const
(points)instance of mesh upon construction
label mergePatchFaces(const scalar minCos, const scalar concaveCos, const label mergeSize, const labelList &patchIDs, const meshRefinement::FaceMergeType mergeType)
Merge coplanar faces if sets are of size mergeSize.
label countHits() const
Count number of intersections (local)
A list of keyword definitions, which are a keyword followed by a number of values (eg...
static void removeFiles(const polyMesh &)
Helper: remove all relevant files from mesh instance.
labelList growFaceCellFace(const labelUList &set) const
List< edge > edgeList
List of edge.
bool getFaceZoneInfo(const word &fzName, label &masterPatchID, label &slavePatchID, surfaceZonesInfo::faceZoneType &fzType) const
Lookup faceZone information. Return false if no information.
static label findRegions(const polyMesh &, const vector &perturbVec, const pointField &locationsInMesh, const pointField &locationsOutsideMesh, const label nRegions, labelList &cellRegion, const boolList &blockedFace, const bool exitIfLeakPath, const refPtr< coordSetWriter > &leakPathFormatter)
Find regions points are in.
bool overwrite() const
Overwrite the mesh?
A 1D array of objects of type <T>, where the size of the vector is known and used for subscript bound...
bool isNormalGap(const scalar planarCos, const label level0, const vector &point0, const vector &normal0, const label level1, const vector &point1, const vector &normal1) const
Is local topology a small gap normal to the test vector.
static label findRegion(const polyMesh &, const labelList &cellRegion, const vector &perturbVec, const point &p)
Find region point is in. Uses optional perturbation to re-test.
autoPtr< mapPolyMesh > doRemoveCells(const labelList &cellsToRemove, const labelList &exposedFaces, const labelList &exposedPatchIDs, removeCells &cellRemover)
Remove cells. Put exposedFaces into exposedPatchIDs.
static ITstream & lookup(const dictionary &dict, const word &keyword, const bool noExit, enum keyType::option matchOpt=keyType::REGEX)
Wrapper around dictionary::lookup which does not exit.
label splitFacesUndo(const labelList &splitFaces, const labelPairList &splits, const dictionary &motionDict, labelList &duplicateFace, List< labelPair > &baffles)
Split faces along diagonal. Maintain mesh quality. Return.
const shellSurfaces & shells() const
Reference to refinement shells (regions)
autoPtr< mapPolyMesh > removeGapCells(const scalar planarAngle, const labelList &minSurfaceLevel, const labelList &globalToMasterPatch, const label growIter)
Detect gapRefinement cells and remove them.
Given list of cells to remove, insert all the topology changes.
autoPtr< mapPolyMesh > createBaffles(const labelList &ownPatch, const labelList &neiPatch)
Create baffle for every internal face where ownPatch != -1.
labelList intersectedFaces() const
Get faces with intersection.
static writeType writeLevel()
Get/set write level.
Class containing mesh-to-mesh mapping information after a mesh distribution where we send parts of me...
label countFaceDirs(const bitSet &isOutsideFace, const label celli) const
Count number of faces on cell that are in set.
Forwards and collection of common point field types.
static bitSet getMasterEdges(const polyMesh &mesh, const labelList &meshEdges)
Determine master edge for subset of edges. If coupled.
labelList getZones(const List< surfaceZonesInfo::faceZoneType > &fzTypes) const
Get zones of given type.
labelList refineCandidates(const pointField &keepPoints, const scalar curvature, const scalar planarAngle, const bool featureRefinement, const bool featureDistanceRefinement, const bool internalRefinement, const bool surfaceRefinement, const bool curvatureRefinement, const bool smallFeatureRefinement, const bool gapRefinement, const bool bigGapRefinement, const bool spreadGapSize, const label maxGlobalCells, const label maxLocalCells) const
Calculate list of cells to refine.
static List< labelPair > subsetBaffles(const polyMesh &mesh, const labelList &zoneIDs, const List< labelPair > &baffles)
Subset baffles according to zones.
set value to -1 any face that was refined
static void calculateEdgeWeights(const polyMesh &mesh, const bitSet &isMasterEdge, const labelList &meshPoints, const edgeList &edges, scalarField &edgeWeights, scalarField &invSumWeight)
Helper: calculate edge weights (1/length)
An enumeration wrapper for classification of a location as being inside/outside of a volume...
labelList intersectedPoints() const
Get points on surfaces with intersection and boundary faces.
autoPtr< mapPolyMesh > directionalRefine(const string &msg, const direction cmpt, const labelList &cellsToRefine)
Directionally refine in direction cmpt.
Class to control time during OpenFOAM simulations that is also the top-level objectRegistry.
A class for managing references or pointers (no reference counting)
Class containing mesh-to-mesh mapping information after a change in polyMesh topology.
labelList collectFaces(const labelList &candidateFaces, const labelHashSet &set) const
Takes mesh with 'baffles' (= boundary faces sharing points). Determines for selected points on bounda...
const fvMesh & mesh() const
Reference to mesh.
Container for data on surfaces used for surface-driven refinement. Contains all the data about the le...
static label addPatch(fvMesh &, const word &name, const dictionary &)
Helper:add patch to mesh. Update all registered fields.
label mergeEdgesUndo(const scalar minCos, const dictionary &motionDict)
Merge edges, maintain mesh quality. Return global number.
static const Enum< writeType > writeTypeNames
Encapsulates queries for features.
static label appendPatch(fvMesh &, const label insertPatchi, const word &, const dictionary &)
Helper:append patch to end of mesh.
void doSplitFaces(const labelList &splitFaces, const labelPairList &splits, polyTopoChange &meshMod) const
Split faces into two.
void mergeFreeStandingBaffles(const snapParameters &snapParams, const bool useTopologicalSnapDetection, const bool removeEdgeConnectedCells, const scalarField &perpendicularAngle, const scalar planarAngle, const dictionary &motionDict, Time &runTime, const labelList &globalToMasterPatch, const labelList &globalToSlavePatch, const pointField &locationsInMesh, const pointField &locationsOutsideMesh)
Merge free-standing baffles.
List< T > values(const HashTable< T, Key, Hash > &tbl, const bool doSort=false)
List of values from HashTable, optionally sorted.
void distribute(const mapDistributePolyMesh &)
Update local numbering for mesh redistribution.
static void updateList(const labelList &newToOld, const T &nullValue, List< T > &elems)
Helper: reorder list according to map.
A list of faces which address into the list of points.
Sends/receives parts of mesh+fvfields to neighbouring processors. Used in load balancing.
vectorField pointField
pointField is a vectorField.
const dimensionedScalar e
Elementary charge.
FaceMergeType
Enumeration for what to do with co-planar patch faces on a single.
word name(const expressions::valueTypeCode typeCode)
A word representation of a valueTypeCode. Empty for INVALID.
ClassName("meshRefinement")
Runtime type information.
Refinement of (split) hexes using polyTopoChange.
A class for handling words, derived from Foam::string.
Field< scalar > scalarField
Specialisation of Field<T> for scalar.
Simple container to keep together snap specific information.
autoPtr< mapPolyMesh > removeLeakCells(const labelList &globalToMasterPatch, const labelList &globalToSlavePatch, const pointField &locationsInMesh, const pointField &locationsOutsideMesh, const labelList &selectedSurfaces)
Remove minimum amount of cells to break any leak from.
void baffleAndSplitMesh(const bool handleSnapProblems, const snapParameters &snapParams, const bool useTopologicalSnapDetection, const bool removeEdgeConnectedCells, const scalarField &perpendicularAngle, const label nErodeCellZones, const dictionary &motionDict, Time &runTime, const labelList &globalToMasterPatch, const labelList &globalToSlavePatch, const pointField &locationsInMesh, const wordList ®ionsInMesh, const pointField &locationsOutsideMesh, const bool exitIfLeakPath, const refPtr< coordSetWriter > &leakPathFormatter)
Split off unreachable areas of mesh.
autoPtr< mapPolyMesh > blockLeakFaces(const labelList &globalToMasterPatch, const labelList &globalToSlavePatch, const pointField &locationsInMesh, const wordList &zonesInMesh, const pointField &locationsOutsideMesh, const labelList &selectedSurfaces)
Baffle faces to break any leak from inside to outside.
autoPtr< mapPolyMesh > mergeZoneBaffles(const bool doInternalZones, const bool doBaffleZones)
Merge all baffles on faceZones.
Base class for writing coordSet(s) and tracks with fields.
static void checkCoupledFaceZones(const polyMesh &)
Helper function: check that face zones are synced.
void testSyncBoundaryFaceList(const scalar mergeDistance, const string &, const UList< T > &, const UList< T > &) const
Compare two lists over all boundary faces.
void growSet(const labelList &neiLevel, const bitSet &isOutsideFace, labelList &refineCell, label &nRefine) const
Add one layer of cells to set.
const labelList & surfaceIndex() const
Per start-end edge the index of the surface hit.
autoPtr< mapPolyMesh > removeLimitShells(const label nBufferLayers, const label nErodeCellZones, const labelList &globalToMasterPatch, const labelList &globalToSlavePatch, const pointField &locationsInMesh, const wordList ®ionsInMesh, const pointField &locationsOutsideMesh)
Remove cells from limitRegions if level -1.
Abstract base class for domain decomposition.
Encapsulates queries for volume refinement ('refine all cells within shell').
autoPtr< mapDistributePolyMesh > balance(const bool keepZoneFaces, const bool keepBaffles, const scalarField &cellWeights, decompositionMethod &decomposer, fvMeshDistribute &distributor)
Redecompose according to cell count.
Container with cells to refine. Refinement given as single direction.
A HashTable similar to std::unordered_map.
autoPtr< mapPolyMesh > doRemovePoints(removePoints &pointRemover, const boolList &pointCanBeDeleted)
labelList directionalRefineCandidates(const label maxGlobalCells, const label maxLocalCells, const labelList ¤tLevel, const direction dir) const
Calculate list of cells to directionally refine.
A HashTable to objects of type <T> with a labelPair key. The hashing is based on labelPair (FixedList...
const refinementFeatures & features() const
Reference to feature edge mesh.
word timeName() const
Replacement for Time::timeName() that returns oldInstance (if overwrite_)
mapType
Enumeration for how userdata is to be mapped upon refinement.
const hexRef8 & meshCutter() const
Reference to meshcutting engine.
static int readFlags(const EnumContainer &namedEnum, const wordList &words)
Helper: convert wordList into bit pattern using provided Enum.
autoPtr< mapDistributePolyMesh > balanceAndRefine(const string &msg, decompositionMethod &decomposer, fvMeshDistribute &distributor, const labelList &cellsToRefine, const scalar maxLoadUnbalance, const label maxCellUnbalance)
Balance before refining some cells.
A Vector of values with scalar precision, where scalar is float/double depending on the compilation f...
void checkData()
Debugging: check that all faces still obey start()>end()
constexpr auto end(C &c) -> decltype(c.end())
Return iterator to the end of the container c.
Database for solution data, solver performance and other reduced data.
const Map< label > & faceToCoupledPatch() const
For faces originating from processor faces store the original.
label addMeshedPatch(const word &name, const dictionary &)
Add patch originating from meshing. Update meshedPatches_.
labelList meshedPatches() const
Get patchIDs for patches added in addMeshedPatch.
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
autoPtr< mapPolyMesh > doRestorePoints(removePoints &pointRemover, const labelList &facesToRestore)
static autoPtr< indirectPrimitivePatch > makePatch(const polyMesh &, const labelList &)
Create patch from set of patches.
void updateMesh(const mapPolyMesh &, const labelList &changedFaces)
Update for external change to mesh. changedFaces are in new mesh.
void updateIntersections(const labelList &changedFaces)
Find any intersection of surface. Store in surfaceIndex_.
scalar mergeDistance() const
autoPtr< mapPolyMesh > splitMeshRegions(const labelList &globalToMasterPatch, const labelList &globalToSlavePatch, const pointField &locationsInMesh, const pointField &locationsOutsideMesh, const bool exitIfLeakPath, const refPtr< coordSetWriter > &leakPathFormatter)
Split mesh. Keep part containing point. Return empty map if.
gmvFile<< "tracers "<< particles.size()<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().x()<< ' ';}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().y()<< ' ';}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().z()<< ' ';}gmvFile<< nl;for(const word &name :lagrangianScalarNames){ IOField< scalar > fld(IOobject(name, runTime.timeName(), cloud::prefix, mesh, IOobject::MUST_READ, IOobject::NO_WRITE))
Helper class which maintains intersections of (changing) mesh with (static) surfaces.
autoPtr< mapPolyMesh > createZoneBaffles(const labelList &zoneIDs, List< labelPair > &baffles, labelList &originatingFaceZone)
Create baffles for faces on faceZones. Return created baffles.
autoPtr< mapPolyMesh > refine(const labelList &cellsToRefine)
Refine some cells.
void setInstance(const fileName &)
Set instance of all local IOobjects.
static void collectAndPrint(const UList< point > &points, const UList< T > &data)
Print list according to (collected and) sorted coordinate.
autoPtr< mapPolyMesh > dupNonManifoldPoints()
Find boundary points that connect to more than one cell.
void storeData(const labelList &pointsToStore, const labelList &facesToStore, const labelList &cellsToStore)
Signal points/face/cells for which to store data.
autoPtr< mapPolyMesh > dupNonManifoldBoundaryPoints()
Find boundary points that are on faceZones of type boundary.
A bitSet stores bits (elements with only two states) in packed internal format and supports a variety...
static void testSyncPointList(const string &msg, const polyMesh &mesh, const List< scalar > &fld)
have slaves (upon refinement) from master
const shellSurfaces & limitShells() const
Reference to limit shells (regions)
static bitSet getMasterPoints(const polyMesh &mesh, const labelList &meshPoints)
Determine master point for subset of points. If coupled.
Mesh data needed to do the Finite Volume discretisation.
Direct mesh changes based on v1.3 polyTopoChange syntax.
void getZoneFaces(const labelList &zoneIDs, labelList &faceZoneID, labelList &ownPatch, labelList &neiPatch, labelList &nBaffles) const
Get per-face information (faceZone, master/slave patch)
option
Enumeration for the data type and search/match modes (bitmask)
const std::string patch
OpenFOAM patch number as a std::string.
void selectIntersectedFaces(const labelList &surfaces, boolList &isBlockedFace) const
Faces currently on boundary or intersected by surface.
autoPtr< mapPolyMesh > mergeBaffles(const List< labelPair > &, const Map< label > &faceToPatch)
Merge baffles. Gets pairs of faces and boundary faces to move.
static const Enum< debugType > debugTypeNames
void selectSeparatedCoupledFaces(boolList &) const
Select coupled faces that are not collocated.
writeType
Enumeration for what to write. Used as a bit-pattern.
faceZoneType
What to do with faceZone faces.
debugType
Enumeration for what to debug. Used as a bit-pattern.
Pointer management similar to std::unique_ptr, with some additional methods and type checking...
Removes selected points from mesh and updates faces using these points.
Mesh consisting of general polyhedral cells.
bool isGap(const scalar, const vector &, const vector &, const vector &, const vector &) const
Is local topology a small gap?
A subset of mesh faces organised as a primitive patch.
Contains information about location on a triSurface.
bool write() const
Write mesh and all data.
List< label > labelList
A List of labels.
static tmp< pointVectorField > makeDisplacementField(const pointMesh &pMesh, const labelList &adaptPatchIDs)
Helper function to make a pointVectorField with correct.
static T gAverage(const bitSet &isMasterElem, const UList< T > &values)
Helper: calculate average.
void dumpIntersections(const fileName &prefix) const
Debug: Write intersection information to OBJ format.
static void weightedSum(const polyMesh &mesh, const bitSet &isMasterEdge, const labelList &meshPoints, const edgeList &edges, const scalarField &edgeWeights, const Field< Type > &data, Field< Type > &sum)
Helper: weighted sum (over all subset of mesh points) by.
labelList detectLeakCells(const boolList &isBlockedFace, const labelList &leakFaces, const labelList &seedCells) const
Return list of cells to block by walking from the seedCells.
mode_t mode(const fileName &name, const bool followLink=true)
Return the file mode, normally following symbolic links.
void dumpRefinementLevel() const
Write refinement level as volScalarFields for postprocessing.
labelList countEdgeFaces(const uindirectPrimitivePatch &pp) const
Count number of faces per patch edge. Parallel consistent.
const volScalarField & p0
An input stream of tokens.
static label findCell(const polyMesh &, const vector &perturbVec, const point &p)
Find cell point is in. Uses optional perturbation to re-test.
uindirectPrimitivePatch pp(UIndirectList< face >(mesh.faces(), faceLabels), mesh.points())
static constexpr const zero Zero
Global zero (0)