Commit 8898bc46 authored by graham's avatar graham
Browse files

ENH: Bigger bounds, adding limitBounds box.

parent bf755c17
......@@ -174,7 +174,6 @@ Foam::tmp<Foam::pointField> Foam::treeBoundBox::points() const
forAll(points, octant)
{
points[octant] = corner(octant);
}
return tPts;
......
......@@ -1219,11 +1219,13 @@ void Foam::conformalVoronoiMesh::insertBoundingPoints()
farPts -= geometryToConformTo_.globalBounds().midpoint();
// Scale the box up
farPts *= 2.0;
farPts *= 10.0;
// Shift corners of bounds back to be relative to midpoint
farPts += geometryToConformTo_.globalBounds().midpoint();
limitBounds_ = treeBoundBox(farPts);
forAll(farPts, fPI)
{
insertPoint(farPts[fPI], Vb::vtFar);
......@@ -1858,6 +1860,7 @@ Foam::conformalVoronoiMesh::conformalVoronoiMesh
cvMeshDict.subDict("motionControl")
),
startOfInternalPoints_(0),
limitBounds_(),
featureVertices_(),
featurePointLocations_(),
featurePointTreePtr_(),
......
......@@ -135,6 +135,9 @@ private:
// careful bookkeeping within several functions.
label startOfInternalPoints_;
//- Limiting bound box before infinity begins
treeBoundBox limitBounds_;
//- Store the feature constraining points to be reinserted after a
// triangulation clear
std::list<Vb> featureVertices_;
......
......@@ -172,7 +172,7 @@ public:
//- Does the Dual vertex form part of a processor patch
inline int parallelDualVertex() const
inline bool parallelDualVertex() const
{
return
(
......@@ -212,7 +212,7 @@ public:
}
// Is the Delaunay cell part of the final dual mesj, i.e. any vertex form
// Is the Delaunay cell part of the final dual mesh, i.e. any vertex form
// part of the internal or boundary definition
inline bool internalOrBoundaryDualVertex() const
{
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment