50 #ifndef primitiveMesh_H 51 #define primitiveMesh_H 82 label nInternalPoints_;
88 mutable label nInternal0Edges_;
91 mutable label nInternal1Edges_;
94 mutable label nInternalEdges_;
97 mutable label nEdges_;
100 label nInternalFaces_;
190 void calcCellShapes()
const;
193 void calcCellCells()
const;
196 void calcPointCells()
const;
199 void calcCells()
const;
202 void calcCellEdges()
const;
205 void calcPointPoints()
const;
210 void calcEdges(
const bool doFaceEdges)
const;
211 void clearOutEdges();
221 static label findFirstCommonElementFromSortedLists
318 const bool detailedReport,
327 const bool detailedReport,
346 const bool detailedReport,
347 const scalar minPyrVol,
379 const scalar warnFlatness,
536 static void calcCells
604 const bool report =
false,
611 const bool report =
false,
618 const bool report =
false,
625 const bool report =
false,
632 const bool report =
false,
646 const bool report =
false,
655 const bool report =
false,
662 const bool report =
false,
669 const bool report =
false,
676 const bool report =
false,
677 const scalar minPyrVol = -SMALL,
684 const bool report =
false,
691 const bool report =
false,
692 const scalar maxSin = 10,
702 const scalar warnFlatness,
711 const scalar reportDistSqr,
719 const scalar minLenSqr,
726 const bool report =
false,
740 virtual bool checkMesh(
const bool report =
false)
const;
765 scalar inflationFraction = 0
const vectorField & faceAreas() const
bool checkFaceAreas(const vectorField &faceAreas, const bool report, const bool detailedReport, labelHashSet *setPtr) const
Check for negative face areas.
const labelListList & cellPoints() const
ClassName("primitiveMesh")
void makeCellCentresAndVols(const vectorField &fCtrs, const vectorField &fAreas, vectorField &cellCtrs, scalarField &cellVols) const
static const unsigned cellsPerPoint_
Estimated number of cells per point.
virtual bool checkPoints(const bool report=false, labelHashSet *setPtr=NULL) const
Check for unused points.
const labelListList & cellCells() const
virtual const pointField & points() const =0
Return mesh points.
bool hasCellPoints() const
static scalar closedThreshold_
Static data to control mesh checking.
label findCell(const point &location) const
Find cell enclosing this location (-1 if not in mesh)
bool hasCellVolumes() const
void calcCellCentresAndVols() const
Calculate cell centres and volumes.
bool checkFacePyramids(const pointField &points, const vectorField &ctrs, const bool report, const bool detailedReport, const scalar minPyrVol, labelHashSet *setPtr) const
Check face pyramid volume.
void clearAddressing()
Clear topological data.
static const unsigned cellsPerEdge_
Estimated number of cells per edge.
const edgeList & edges() const
Return mesh edges. Uses calcEdges.
const labelListList & pointEdges() const
bool checkCellVolumes(const scalarField &vols, const bool report, const bool detailedReport, labelHashSet *setPtr) const
Check for negative cell volumes.
bool checkDuplicateFaces(const label, const Map< label > &, label &nBaffleFaces, labelHashSet *) const
Check if all points on face are shared with another face.
static const unsigned edgesPerFace_
Estimated number of edges per cell.
bool pointInCell(const point &p, label celli) const
Return true if the point is in the cell.
bool hasCellEdges() const
A simple container for copying or transferring objects of type <T>.
static scalar aspectThreshold_
Aspect ratio warning threshold.
bool checkClosedBoundary(const vectorField &, const bool, const PackedBoolList &) const
Check boundary for closedness.
virtual bool checkEdgeLength(const bool report, const scalar minLenSqr, labelHashSet *setPtr=NULL) const
Check edge length.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
static scalar planarCosAngle_
Threshold where faces are considered coplanar.
virtual bool checkCellsZipUp(const bool report=false, labelHashSet *setPtr=NULL) const
Check cell zip-up.
const cellList & cells() const
const labelListList & pointCells() const
const scalarField & cellVols
void printAllocated() const
Print a list of all the currently allocated mesh data.
const vectorField & cellCentres() const
const cellShapeList & cellShapes() const
Return cell shapes.
virtual const faceList & faces() const =0
Return faces.
bool checkFaceAngles(const pointField &points, const vectorField &faceAreas, const bool report, const scalar maxDeg, labelHashSet *setPtr) const
Check face angles.
label nInternal1Edges() const
Internal edges using 0 or 1 boundary point.
bool hasCellCentres() const
static const unsigned facesPerPoint_
Estimated number of faces per point.
const labelListList & faceEdges() const
bool checkConcaveCells(const vectorField &fAreas, const pointField &fCentres, const bool report, labelHashSet *setPtr) const
Check for concave cells by the planes of faces.
const labelListList & cellEdges() const
static const unsigned pointsPerCell_
Estimated number of points per cell.
bool hasPointFaces() const
bool hasPointCells() const
label findNearestCell(const point &location) const
Find the cell with the nearest cell centre to location.
Cell-face mesh analysis engine.
label nInternal0Edges() const
Internal edges (i.e. not on boundary face) using.
virtual bool checkUpperTriangular(const bool report=false, labelHashSet *setPtr=NULL) const
Check face ordering.
void makeFaceCentresAndAreas(const pointField &p, vectorField &fCtrs, vectorField &fAreas) const
primitiveMesh()
Construct null.
static scalar setNonOrthThreshold(const scalar)
Set the non-orthogonality warning threshold in degrees.
void reset(const label nPoints, const label nInternalFaces, const label nFaces, const label nCells)
Reset this primitiveMesh given the primitive array sizes.
static const unsigned edgesPerCell_
Estimated number of edges per cell.
bool checkClosedCells(const vectorField &faceAreas, const scalarField &cellVolumes, const bool report, labelHashSet *setPtr, labelHashSet *aspectSetPtr, const Vector< label > &meshD) const
Check cells for closedness.
static scalar skewThreshold_
Skewness warning threshold.
static scalar setSkewThreshold(const scalar)
Set the skewness warning threshold as percentage.
static const unsigned facesPerCell_
Estimated number of faces per cell.
virtual const pointField & oldPoints() const =0
Return old points for mesh motion.
static scalar nonOrthThreshold_
Non-orthogonality warning threshold in deg.
virtual bool checkFaceVertices(const bool report=false, labelHashSet *setPtr=NULL) const
Check uniqueness of face vertices.
virtual bool checkPointNearness(const bool report, const scalar reportDistSqr, labelHashSet *setPtr=NULL) const
Check for point-point-nearness,.
bool hasCellShapes() const
static const unsigned facesPerEdge_
Estimated number of faces per edge.
void clearOut()
Clear all geometry and addressing unnecessary for CFD.
virtual ~primitiveMesh()
Destructor.
static const unsigned pointsPerFace_
Estimated number of points per face.
const labelListList & pointFaces() const
static scalar setAspectThreshold(const scalar)
Set the aspect ratio warning threshold.
void calcFaceCentresAndAreas() const
Calculate face centres and areas.
bool hasFaceAreas() const
bool checkFaceOrthogonality(const vectorField &fAreas, const vectorField &cellCtrs, const bool report, labelHashSet *setPtr) const
Check for non-orthogonality.
label nInternalFaces() const
bool checkFaceFlatness(const pointField &points, const vectorField &faceCentres, const vectorField &faceAreas, const bool report, const scalar warnFlatness, labelHashSet *setPtr) const
Check face warpage.
static scalar setClosedThreshold(const scalar)
Set the closedness ratio warning threshold.
virtual bool checkFaceFaces(const bool report=false, labelHashSet *setPtr=NULL) const
Check face-face connectivity.
virtual bool checkMesh(const bool report=false) const
Check mesh for correctness. Returns false for no error.
bool hasEdgeCells() const
const labelListList & edgeFaces() const
bool checkFaceSkewness(const pointField &points, const vectorField &fCtrs, const vectorField &fAreas, const vectorField &cellCtrs, const bool report, labelHashSet *setPtr) const
Check face skewness.
A 1D vector of objects of type <T>, where the size of the vector is known and can be used for subscri...
label nInternalEdges() const
Internal edges using 0,1 or 2 boundary points.
tmp< scalarField > movePoints(const pointField &p, const pointField &oldP)
Move points, returns volumes swept by faces in motion.
bool pointInCellBB(const point &p, label celli, scalar inflationFraction=0) const
Return true if the point in the cell bounding box.
bool hasPointEdges() const
void clearGeom()
Clear geometry.
static const unsigned edgesPerPoint_
Estimated number of edges per point.
static bool calcPointOrder(label &nInternalPoints, labelList &pointMap, const faceList &, const label nInternalFaces, const label nPoints)
Helper function to calculate point ordering. Returns true.
void calcEdgeVectors() const
Calculate edge vectors.
const labelListList & edgeCells() const
const labelListList & pointPoints() const
bool checkCommonOrder(const label, const Map< label > &, labelHashSet *) const
Check that shared points are in consecutive order.
const scalarField & cellVolumes() const
bool hasPointPoints() const
virtual bool checkTopology(const bool report=false) const
Check mesh topology for correctness.
bool hasFaceEdges() const
bool hasCellCells() const
virtual const labelList & faceOwner() const =0
Face face-owner addresing.
bool hasEdgeFaces() const
virtual const labelList & faceNeighbour() const =0
Face face-neighbour addressing.
A class for managing temporary objects.
label nInternalPoints() const
Points not on boundary.
virtual bool checkGeometry(const bool report=false) const
Check mesh geometry (& implicitly topology) for correctness.
bool isInternalFace(const label faceIndex) const
Return true if given face label is internal to the mesh.
const vectorField & faceCentres() const
bool hasFaceCentres() const