zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

ENH: refinementFeatures: use extendedEdgeMesh feature points

Browse Source
This commit is contained in:
mattijs
2013-12-20 16:32:34 +00:00
parent b75620e3f6
commit 386696d55e
3 changed files with 11 additions and 30 deletions
Download Patch File Download Diff File Expand all files Collapse all files

7
src/mesh/autoMesh/autoHexMesh/meshRefinement/meshRefinementRefine.C
View File

@ -35,7 +35,6 @@ License
#include "fvMeshDistribute.H"
#include "polyTopoChange.H"
#include "mapDistributePolyMesh.H"
#include "featureEdgeMesh.H"
#include "Cloud.H"
//#include "globalIndex.H"
#include "OBJstream.H"
@ -358,7 +357,7 @@ void Foam::meshRefinement::markFeatureCellLevel
forAll(features_, featI)
{
const featureEdgeMesh& featureMesh = features_[featI];
const edgeMesh& featureMesh = features_[featI];
const label featureLevel = features_.levels()[featI][0];
const labelListList& pointEdges = featureMesh.pointEdges();
@ -505,7 +504,7 @@ void Foam::meshRefinement::markFeatureCellLevel
label featI = startTp.i();
label pointI = startTp.j();
const featureEdgeMesh& featureMesh = features_[featI];
const edgeMesh& featureMesh = features_[featI];
const labelList& pEdges = featureMesh.pointEdges()[pointI];
// Now shoot particles down all pEdges.
@ -561,7 +560,7 @@ void Foam::meshRefinement::markFeatureCellLevel
label featI = tp.i();
label pointI = tp.j();
const featureEdgeMesh& featureMesh = features_[featI];
const edgeMesh& featureMesh = features_[featI];
const labelList& pEdges = featureMesh.pointEdges()[pointI];
// Particle now at pointI. Check connected edges to see which one

32
src/mesh/autoMesh/autoHexMesh/refinementFeatures/refinementFeatures.C
View File

@ -186,7 +186,7 @@ void Foam::refinementFeatures::read
set(featI, new extendedFeatureEdgeMesh(featObj, eeMesh));
}
const edgeMesh& eMesh = operator[](featI);
const extendedEdgeMesh& eMesh = operator[](featI);
//eMesh.mergePoints(meshRefiner_.mergeDistance());
@ -258,13 +258,9 @@ void Foam::refinementFeatures::read
}
void Foam::refinementFeatures::buildTrees
(
const label featI,
const labelList& featurePoints
)
void Foam::refinementFeatures::buildTrees(const label featI)
{
const edgeMesh& eMesh = operator[](featI);
const extendedEdgeMesh& eMesh = operator[](featI);
const pointField& points = eMesh.points();
const edgeList& edges = eMesh.edges();
@ -299,6 +295,9 @@ void Foam::refinementFeatures::buildTrees
)
);
labelList featurePoints(identity(eMesh.nonFeatureStart()));
pointTrees_.set
(
featI,
@ -408,24 +407,7 @@ Foam::refinementFeatures::refinementFeatures
// Search engines
forAll(*this, i)
{
const extendedEdgeMesh& eMesh = operator[](i);
const labelListList& pointEdges = eMesh.pointEdges();
DynamicList<label> featurePoints;
forAll(pointEdges, pointI)
{
if (pointEdges[pointI].size() > 2)
{
featurePoints.append(pointI);
}
}
Info<< "Detected " << featurePoints.size()
<< " featurePoints out of " << pointEdges.size()
<< " points on feature " << operator[](i).name()
<< endl;
buildTrees(i, featurePoints);
buildTrees(i);
}
}

2
src/mesh/autoMesh/autoHexMesh/refinementFeatures/refinementFeatures.H
View File

@ -76,7 +76,7 @@ private:
void read(const objectRegistry&, const PtrList<dictionary>&);
//- Build edge tree and feature point tree
void buildTrees(const label, const labelList&);
void buildTrees(const label);
//- Find shell level higher than ptLevel
void findHigherLevel