All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
polyMeshAdder Class Reference

Adds two meshes without using any polyMesh morphing. More...

Inheritance diagram for polyMeshAdder:

Static Public Member Functions

static autoPtr< polyMeshadd (const IOobject &io, const polyMesh &mesh0, const polyMesh &mesh1, const faceCoupleInfo &coupleInfo, autoPtr< mapAddedPolyMesh > &mapPtr)
 Add two polyMeshes. Returns new polyMesh and map construct. More...
 
static autoPtr< mapAddedPolyMeshadd (polyMesh &mesh0, const polyMesh &mesh1, const faceCoupleInfo &coupleInfo, const bool validBoundary=true)
 Inplace add mesh to polyMesh. Returns map construct. More...
 
static Map< labelfindSharedPoints (const polyMesh &, const scalar mergeTol)
 Find topologically and geometrically shared points. More...
 
static void mergePoints (const polyMesh &, const Map< label > &pointToMaster, polyTopoChange &meshMod)
 Helper: Merge points. More...
 

Detailed Description

Adds two meshes without using any polyMesh morphing.

Gets faces to couple as faceCoupleInfo which is list of faces on both meshes. Returns map from last mesh addition.

Source files

Definition at line 59 of file polyMeshAdder.H.

Member Function Documentation

◆ add() [1/2]

Foam::autoPtr< Foam::polyMesh > add ( const IOobject io,
const polyMesh mesh0,
const polyMesh mesh1,
const faceCoupleInfo coupleInfo,
autoPtr< mapAddedPolyMesh > &  mapPtr 
)
static

Add two polyMeshes. Returns new polyMesh and map construct.

Definition at line 1465 of file polyMeshAdder.C.

References polyMesh::addPatches(), Foam::DelaunayMeshTools::allPoints(), polyMesh::boundaryMesh(), faceCoupleInfo::cutFaces(), Foam::identity(), primitiveMesh::nCells(), primitiveMesh::nFaces(), primitiveMesh::nPoints(), autoPtr< T >::reset(), and UPtrList< T >::size().

Referenced by fvMeshAdder::add().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ add() [2/2]

Foam::autoPtr< Foam::mapAddedPolyMesh > add ( polyMesh mesh0,
const polyMesh mesh1,
const faceCoupleInfo coupleInfo,
const bool  validBoundary = true 
)
static

◆ findSharedPoints()

Foam::Map< Foam::label > findSharedPoints ( const polyMesh mesh,
const scalar  mergeTol 
)
static

Find topologically and geometrically shared points.

  • should only be called for parallel correct mesh (since uses mesh.globalData)
  • returns Map from point to master point (all in mesh point labels) for any sets of points that need to be merged.

Old: geometric merging. Causes problems for two close shared points.

labelList sharedToMerged; label nUnique = Foam::mergePoints ( pointField ( mesh.points(), sharedPointLabels ), mergeDist, false, sharedToMerged );

Find out which sets of points get merged and create a map from mesh point to unique point.

Map<label> pointToMaster(10*sharedToMerged.size());

if (nUnique < sharedPointLabels.size()) { labelListList mergeSets ( invertOneToMany ( sharedToMerged.size(), sharedToMerged ) );

label nMergeSets = 0;

forAll(mergeSets, setI) { const labelList& mergeSet = mergeSets[setI];

if (mergeSet.size() > 1) { Take as master the shared point with the lowest mesh point label. (rather arbitrarily - could use max or any other one of the points)

nMergeSets++;

label masterI = labelMax;

forAll(mergeSet, i) { label sharedI = mergeSet[i];

masterI = min(masterI, sharedPointLabels[sharedI]); }

forAll(mergeSet, i) { label sharedI = mergeSet[i];

pointToMaster.insert(sharedPointLabels[sharedI], masterI); } } }

if (debug) { Pout<< "polyMeshAdder : merging:" << pointToMaster.size() << " into " << nMergeSets << " sets." << endl; } }

Definition at line 2002 of file polyMeshAdder.C.

References Foam::abort(), Foam::endl(), Foam::FatalError, FatalErrorInFunction, Foam::findIndex(), forAll, forAllConstIter(), polyMesh::globalData(), HashTable< T, Key, Hash >::insert(), Foam::invertOneToMany(), Foam::labelMax, Foam::mergePoints(), polyMeshAdder::mergePoints(), Foam::min(), polyMesh::points(), List< T >::setSize(), globalMeshData::sharedPointAddr(), globalMeshData::sharedPointLabels(), and List< T >::size().

Referenced by polyMeshAdder::add().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ mergePoints()

void mergePoints ( const polyMesh mesh,
const Map< label > &  pointToMaster,
polyTopoChange meshMod 
)
static

Helper: Merge points.

  • Gets map from point to destination point
  • Removes all points that don't map to themselves
  • Modifies all faces that use the points to be removed.

Definition at line 2214 of file polyMeshAdder.C.

References polyMesh::boundaryMesh(), HashTableCore::end(), f(), polyMesh::faceNeighbour(), polyMesh::faceOwner(), polyMesh::faces(), polyMesh::faceZones(), HashTable< T, Key, Hash >::find(), faceZone::flipMap(), forAll, polyMesh::points(), polyTopoChange::removePoint(), polyTopoChange::setAction(), faceZone::whichFace(), polyBoundaryMesh::whichPatch(), and ZoneMesh< ZoneType, MeshType >::whichZone().

Referenced by polyMeshAdder::findSharedPoints(), and fvMeshDistribute::printCoupleInfo().

Here is the call graph for this function:
Here is the caller graph for this function:

The documentation for this class was generated from the following files: