FreeFOAM The Cross-Platform CFD Toolkit
Hosted by SourceForge:
Get FreeFOAM at SourceForge.net.
            Fast, secure and Free Open Source software downloads

polyMeshAdder Class Reference

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

#include <dynamicMesh/polyMeshAdder.H>


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 60 of file polyMeshAdder.H.

Inheritance diagram for polyMeshAdder:

List of all members.

Static Public Member Functions

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

Member Function Documentation

Foam::autoPtr< Foam::polyMesh > add ( const IOobject &   io,
const polyMesh &   mesh0,
const polyMesh &   mesh1,
const faceCoupleInfo &   coupleInfo,
autoPtr< mapAddedPolyMesh > &   mapPtr  
) [static]
Foam::Map< Foam::label > findSharedPoints ( const polyMesh &   mesh,
const scalar   mergeTol  
) [static]

Find topologically&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; pointField mergedPoints; bool hasMerged = Foam::mergePoints ( pointField ( mesh.points(), sharedPointLabels ), mergeDist, false, sharedToMerged, mergedPoints );

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 (hasMerged) { 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 1785 of file polyMeshAdder.C.

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

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

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