mirror of
https://develop.openfoam.com/Development/openfoam.git
synced 2025-11-28 03:28:01 +00:00
ENH: Added new processor level of detail (LOD) to optimise parallel data transfers
This commit is contained in:
Andrew Heather
@ -231,6 +231,11 @@ triSurface/triSurfaceTools/pointToPointPlanarInterpolation.C
|
|||||||
|
|
||||||
twoDPointCorrector/twoDPointCorrector.C
|
twoDPointCorrector/twoDPointCorrector.C
|
||||||
|
|
||||||
|
processorLOD/processorLOD/processorLOD.C
|
||||||
|
processorLOD/box/box.C
|
||||||
|
processorLOD/cellBox/cellBox.C
|
||||||
|
processorLOD/faceBox/faceBox.C
|
||||||
|
|
||||||
AMI=AMIInterpolation
|
AMI=AMIInterpolation
|
||||||
$(AMI)/AMIInterpolation/AMIInterpolationName.C
|
$(AMI)/AMIInterpolation/AMIInterpolationName.C
|
||||||
$(AMI)/AMIInterpolation/AMIPatchToPatchInterpolation.C
|
$(AMI)/AMIInterpolation/AMIPatchToPatchInterpolation.C
|
||||||
|
|||||||
635
src/meshTools/processorLOD/box/box.C
Normal file
635
src/meshTools/processorLOD/box/box.C
Normal file
@ -0,0 +1,635 @@
|
|||||||
|
/*---------------------------------------------------------------------------*\
|
||||||
|
========= |
|
||||||
|
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||||
|
\\ / O peration |
|
||||||
|
\\ / A nd | Copyright (C) 2017 OpenCFD Ltd.
|
||||||
|
\\/ M anipulation |
|
||||||
|
-------------------------------------------------------------------------------
|
||||||
|
License
|
||||||
|
This file is part of OpenFOAM.
|
||||||
|
|
||||||
|
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||||
|
under the terms of the GNU General Public License as published by
|
||||||
|
the Free Software Foundation, either version 3 of the License, or
|
||||||
|
(at your option) any later version.
|
||||||
|
|
||||||
|
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||||
|
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||||
|
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||||
|
for more details.
|
||||||
|
|
||||||
|
You should have received a copy of the GNU General Public License
|
||||||
|
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||||
|
|
||||||
|
\*---------------------------------------------------------------------------*/
|
||||||
|
|
||||||
|
#include "box.H"
|
||||||
|
#include "mapDistribute.H"
|
||||||
|
#include "OFstream.H"
|
||||||
|
#include "meshTools.H"
|
||||||
|
|
||||||
|
const Foam::label Foam::processorLODs::box::DROP = 0;
|
||||||
|
const Foam::label Foam::processorLODs::box::REFINE = 1;
|
||||||
|
const Foam::label Foam::processorLODs::box::FIXED = 2;
|
||||||
|
const Foam::label Foam::processorLODs::box::nStartUpIter = 2;
|
||||||
|
|
||||||
|
namespace Foam
|
||||||
|
{
|
||||||
|
namespace processorLODs
|
||||||
|
{
|
||||||
|
defineTypeNameAndDebug(box, 0);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * Protected Member Functions * * * * * * * * * * * //
|
||||||
|
|
||||||
|
void Foam::processorLODs::box::writeBoxes
|
||||||
|
(
|
||||||
|
const List<DynamicList<treeBoundBox>>& fixedBoxes,
|
||||||
|
const label iter
|
||||||
|
) const
|
||||||
|
{
|
||||||
|
static label time = 0;
|
||||||
|
|
||||||
|
OFstream os
|
||||||
|
(
|
||||||
|
"processor" + Foam::name(Pstream::myProcNo())
|
||||||
|
+ "_time" + Foam::name(time)
|
||||||
|
+ "_iter" + Foam::name(iter) + ".obj"
|
||||||
|
);
|
||||||
|
|
||||||
|
label verti = 0;
|
||||||
|
for (label proci = 0; proci < Pstream::nProcs(); ++proci)
|
||||||
|
{
|
||||||
|
if (proci == Pstream::myProcNo())
|
||||||
|
{
|
||||||
|
continue;
|
||||||
|
}
|
||||||
|
|
||||||
|
const DynamicList<treeBoundBox>& procBoxes = fixedBoxes[proci];
|
||||||
|
forAll(procBoxes, boxi)
|
||||||
|
{
|
||||||
|
const treeBoundBox& bb = procBoxes[boxi];
|
||||||
|
|
||||||
|
// Write the points
|
||||||
|
const pointField pts(bb.points());
|
||||||
|
for (const point& p : pts)
|
||||||
|
{
|
||||||
|
meshTools::writeOBJ(os, p);
|
||||||
|
}
|
||||||
|
|
||||||
|
// Write the box faces
|
||||||
|
for (const face& f : bb.faces)
|
||||||
|
{
|
||||||
|
os << 'f';
|
||||||
|
for (const label fpi : f)
|
||||||
|
{
|
||||||
|
os << ' ' << fpi + verti + 1;
|
||||||
|
}
|
||||||
|
os << nl;
|
||||||
|
}
|
||||||
|
verti += pts.size();
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
++time;
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
void Foam::processorLODs::box::setRefineFlags
|
||||||
|
(
|
||||||
|
const label refineIter,
|
||||||
|
const label nTgtObjects,
|
||||||
|
List<labelHashSet>& fixedSendElems,
|
||||||
|
List<List<labelList>>& localTgtElems,
|
||||||
|
List<labelList>& refineFlags,
|
||||||
|
labelList& nElems
|
||||||
|
) const
|
||||||
|
{
|
||||||
|
PstreamBuffers pBufs(Pstream::commsTypes::nonBlocking);
|
||||||
|
|
||||||
|
// Identify src boxes that can be refined and send to all remote procs
|
||||||
|
for (label proci = 0; proci < Pstream::nProcs(); ++proci)
|
||||||
|
{
|
||||||
|
if (proci != Pstream::myProcNo())
|
||||||
|
{
|
||||||
|
UOPstream toProc(proci, pBufs);
|
||||||
|
toProc << nObjectsOfType_ << boxes_[proci] << newToOld_[proci];
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
pBufs.finishedSends();
|
||||||
|
|
||||||
|
// Test each remote src bb with local tgt objects to identify which remote
|
||||||
|
// src boxes can/should be refined
|
||||||
|
for (label proci = 0; proci < Pstream::nProcs(); ++proci)
|
||||||
|
{
|
||||||
|
if (proci == Pstream::myProcNo())
|
||||||
|
{
|
||||||
|
// Not refining boxes I send to myself - will be sending all local
|
||||||
|
// elements
|
||||||
|
continue;
|
||||||
|
}
|
||||||
|
|
||||||
|
// Receive the subset of changing src bound boxes for proci
|
||||||
|
UIPstream fromProc(proci, pBufs);
|
||||||
|
const label nObjects = readLabel(fromProc);
|
||||||
|
List<treeBoundBox> remoteSrcBoxes(fromProc);
|
||||||
|
const List<label> newToOld(fromProc);
|
||||||
|
|
||||||
|
if (remoteSrcBoxes.empty())
|
||||||
|
{
|
||||||
|
continue;
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
labelList& procRefineFlags = refineFlags[proci];
|
||||||
|
procRefineFlags.setSize(remoteSrcBoxes.size(), FIXED);
|
||||||
|
|
||||||
|
if (scalar(nTgtObjects)/scalar(nObjects) < 0.1)
|
||||||
|
{
|
||||||
|
// Sending less than 10% of objects of this type
|
||||||
|
// - shortcut by sending all
|
||||||
|
fixedSendElems[proci].insert(identity(nTgtObjects));
|
||||||
|
nElems[proci] = nTgtObjects;
|
||||||
|
|
||||||
|
continue;
|
||||||
|
}
|
||||||
|
|
||||||
|
label nElem = 0;
|
||||||
|
List<labelList>& localProcTgtElems = localTgtElems[proci];
|
||||||
|
List<labelList> newLocalProcTgtElems(remoteSrcBoxes.size());
|
||||||
|
|
||||||
|
forAll(remoteSrcBoxes, srcBoxi)
|
||||||
|
{
|
||||||
|
const treeBoundBox& remSrcBb = remoteSrcBoxes[srcBoxi];
|
||||||
|
DynamicList<label> selectedElems(maxObjectsPerLeaf_);
|
||||||
|
|
||||||
|
if (refineIter > nStartUpIter)
|
||||||
|
{
|
||||||
|
// Iterate over cached subset of tgt elements
|
||||||
|
const label oldBoxi = newToOld[srcBoxi];
|
||||||
|
const labelList& tgtBoxElems = localProcTgtElems[oldBoxi];
|
||||||
|
|
||||||
|
for (const label tgtObji : tgtBoxElems)
|
||||||
|
{
|
||||||
|
if (calcTgtBox(tgtObji).overlaps(remSrcBb))
|
||||||
|
{
|
||||||
|
selectedElems.append(tgtObji);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
}
|
||||||
|
else
|
||||||
|
{
|
||||||
|
// Iterating over all target elements
|
||||||
|
for (label tgtObji = 0; tgtObji < nTgtObjects; ++tgtObji)
|
||||||
|
{
|
||||||
|
if (calcTgtBox(tgtObji).overlaps(remSrcBb))
|
||||||
|
{
|
||||||
|
selectedElems.append(tgtObji);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
nElem += selectedElems.size();
|
||||||
|
|
||||||
|
if
|
||||||
|
(
|
||||||
|
proci == Pstream::myProcNo()
|
||||||
|
|| selectedElems.size() < maxObjectsPerLeaf_
|
||||||
|
)
|
||||||
|
{
|
||||||
|
procRefineFlags[srcBoxi] = FIXED;
|
||||||
|
fixedSendElems[proci].insert(selectedElems);
|
||||||
|
}
|
||||||
|
else
|
||||||
|
{
|
||||||
|
procRefineFlags[srcBoxi] = REFINE;
|
||||||
|
if (refineIter >= nStartUpIter)
|
||||||
|
{
|
||||||
|
newLocalProcTgtElems[srcBoxi].transfer(selectedElems);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
localProcTgtElems.transfer(newLocalProcTgtElems);
|
||||||
|
nElems[proci] = nElem;
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
void Foam::processorLODs::box::refineBox
|
||||||
|
(
|
||||||
|
const label boxi,
|
||||||
|
const label refineIter,
|
||||||
|
const label nSrcElem,
|
||||||
|
const treeBoundBox& origBox,
|
||||||
|
DynamicList<treeBoundBox>& procBoxes,
|
||||||
|
DynamicList<labelList>& procBoxElems,
|
||||||
|
DynamicList<label>& procNewToOld
|
||||||
|
) const
|
||||||
|
{
|
||||||
|
// Create the new boxes
|
||||||
|
|
||||||
|
if (refineIter == nStartUpIter)
|
||||||
|
{
|
||||||
|
// Start caching the addressing
|
||||||
|
for (direction octant = 0; octant < 8; ++octant)
|
||||||
|
{
|
||||||
|
const treeBoundBox subBb(origBox.subBbox(octant));
|
||||||
|
|
||||||
|
// Identify the src elements for each box
|
||||||
|
DynamicList<label> newElems(nSrcElem/2);
|
||||||
|
|
||||||
|
for (label srcElemi = 0; srcElemi < nSrcElem; ++srcElemi)
|
||||||
|
{
|
||||||
|
if (subBb.overlaps(calcSrcBox(srcElemi)))
|
||||||
|
{
|
||||||
|
newElems.append(srcElemi);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
if (newElems.size())
|
||||||
|
{
|
||||||
|
procBoxes.append(subBb);
|
||||||
|
procBoxElems.append(newElems);
|
||||||
|
procNewToOld.append(boxi);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
}
|
||||||
|
else
|
||||||
|
{
|
||||||
|
for (direction octant = 0; octant < 8; ++octant)
|
||||||
|
{
|
||||||
|
const treeBoundBox subBb(origBox.subBbox(octant));
|
||||||
|
|
||||||
|
for (label srcElemi = 0; srcElemi < nSrcElem; ++srcElemi)
|
||||||
|
{
|
||||||
|
if (subBb.overlaps(calcSrcBox(srcElemi)))
|
||||||
|
{
|
||||||
|
procBoxes.append(subBb);
|
||||||
|
break;
|
||||||
|
}
|
||||||
|
}
|
||||||
|
}
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
void Foam::processorLODs::box::refineBox
|
||||||
|
(
|
||||||
|
const label boxi,
|
||||||
|
const labelList& srcAddr,
|
||||||
|
const treeBoundBox& origBox,
|
||||||
|
DynamicList<treeBoundBox>& procBoxes,
|
||||||
|
DynamicList<labelList>& procBoxElems,
|
||||||
|
DynamicList<label>& procNewToOld
|
||||||
|
) const
|
||||||
|
{
|
||||||
|
// Create the new boxes
|
||||||
|
for (direction octant = 0; octant < 8; ++octant)
|
||||||
|
{
|
||||||
|
const treeBoundBox subBb(origBox.subBbox(octant));
|
||||||
|
|
||||||
|
// Identify the src elements for each box
|
||||||
|
DynamicList<label> newElems(srcAddr.size()/2);
|
||||||
|
|
||||||
|
for (const label srcElemi : srcAddr)
|
||||||
|
{
|
||||||
|
if (subBb.overlaps(calcSrcBox(srcElemi)))
|
||||||
|
{
|
||||||
|
newElems.append(srcElemi);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
// Only keeping new box if it overlaps src objects
|
||||||
|
if (newElems.size())
|
||||||
|
{
|
||||||
|
procBoxes.append(subBb);
|
||||||
|
procBoxElems.append(newElems);
|
||||||
|
procNewToOld.append(boxi);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
bool Foam::processorLODs::box::doRefineBoxes
|
||||||
|
(
|
||||||
|
const label refineIter,
|
||||||
|
const label nSrcFaces,
|
||||||
|
const List<labelList>& refineFlags,
|
||||||
|
const labelList& nElems,
|
||||||
|
List<DynamicList<treeBoundBox>>& fixedBoxes
|
||||||
|
)
|
||||||
|
{
|
||||||
|
PstreamBuffers pBufs(Pstream::commsTypes::nonBlocking);
|
||||||
|
|
||||||
|
// Send refine info back for list of evolving src boxes
|
||||||
|
forAll(refineFlags, proci)
|
||||||
|
{
|
||||||
|
if (proci != Pstream::myProcNo())
|
||||||
|
{
|
||||||
|
UOPstream toProc(proci, pBufs);
|
||||||
|
toProc << nElems[proci] << refineFlags[proci];
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
pBufs.finishedSends();
|
||||||
|
|
||||||
|
// Receive refine refinement actions from remote procs and use to refine
|
||||||
|
// local src boxes
|
||||||
|
bool refineBoxes = false;
|
||||||
|
label nElem = 0;
|
||||||
|
List<DynamicList<label>> newToOld(Pstream::nProcs());
|
||||||
|
|
||||||
|
|
||||||
|
for (label proci = 0; proci < Pstream::nProcs(); ++proci)
|
||||||
|
{
|
||||||
|
if (proci == Pstream::myProcNo())
|
||||||
|
{
|
||||||
|
continue;
|
||||||
|
}
|
||||||
|
|
||||||
|
UIPstream fromProc(proci, pBufs);
|
||||||
|
nElem += readLabel(fromProc);
|
||||||
|
const labelList refineFlags(fromProc);
|
||||||
|
|
||||||
|
const List<treeBoundBox>& procBoxes = boxes_[proci];
|
||||||
|
DynamicList<treeBoundBox> newProcBoxes(procBoxes.size());
|
||||||
|
DynamicList<labelList> newProcBoxElems(procBoxes.size());
|
||||||
|
newToOld[proci].setCapacity(boxes_[proci].size());
|
||||||
|
DynamicList<label>& newProcNewToOld = newToOld[proci];
|
||||||
|
|
||||||
|
|
||||||
|
forAll(procBoxes, boxi)
|
||||||
|
{
|
||||||
|
if (refineFlags[boxi] == DROP)
|
||||||
|
{
|
||||||
|
// Skip box
|
||||||
|
}
|
||||||
|
else if (refineFlags[boxi] == REFINE)
|
||||||
|
{
|
||||||
|
if (refineIter > nStartUpIter)
|
||||||
|
{
|
||||||
|
// Can use locally cached info to speed up intersection
|
||||||
|
// tests
|
||||||
|
refineBox
|
||||||
|
(
|
||||||
|
boxi,
|
||||||
|
boxSrcElems_[proci][boxi],
|
||||||
|
procBoxes[boxi],
|
||||||
|
newProcBoxes,
|
||||||
|
newProcBoxElems,
|
||||||
|
newProcNewToOld
|
||||||
|
);
|
||||||
|
}
|
||||||
|
else
|
||||||
|
{
|
||||||
|
refineBox
|
||||||
|
(
|
||||||
|
boxi,
|
||||||
|
refineIter,
|
||||||
|
nSrcFaces,
|
||||||
|
procBoxes[boxi],
|
||||||
|
newProcBoxes,
|
||||||
|
newProcBoxElems,
|
||||||
|
newProcNewToOld
|
||||||
|
);
|
||||||
|
}
|
||||||
|
|
||||||
|
refineBoxes = true;
|
||||||
|
}
|
||||||
|
else if (refineFlags[boxi] == FIXED)
|
||||||
|
{
|
||||||
|
fixedBoxes[proci].append(procBoxes[boxi]);
|
||||||
|
}
|
||||||
|
else
|
||||||
|
{
|
||||||
|
FatalErrorInFunction
|
||||||
|
<< "Unhandled refine action " << refineFlags[boxi]
|
||||||
|
<< abort(FatalError);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
// Only keeping boxes that are to be refined
|
||||||
|
boxes_[proci].transfer(newProcBoxes);
|
||||||
|
boxSrcElems_[proci].transfer(newProcBoxElems);
|
||||||
|
newToOld_[proci].transfer(newProcNewToOld);
|
||||||
|
}
|
||||||
|
|
||||||
|
return returnReduce(refineBoxes, orOp<bool>());
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
Foam::autoPtr<Foam::mapDistribute> Foam::processorLODs::box::createMap
|
||||||
|
(
|
||||||
|
const label nSrcElems,
|
||||||
|
const label nTgtElems
|
||||||
|
)
|
||||||
|
{
|
||||||
|
// Store elements to send - will be used to build the mapDistribute
|
||||||
|
List<labelHashSet> fixedSendElems(Pstream::nProcs());
|
||||||
|
|
||||||
|
// List of local tgt elems to optimise searching fot tgt elements inside
|
||||||
|
// remote src boxes
|
||||||
|
List<List<labelList>> localTgtElems(Pstream::nProcs());
|
||||||
|
|
||||||
|
// Storage of boxes per processor - only useful for debugging
|
||||||
|
List<DynamicList<treeBoundBox>> fixedBoxes(Pstream::nProcs());
|
||||||
|
|
||||||
|
// Iterate to subdivide src boxes
|
||||||
|
label refinementIter = 1;
|
||||||
|
bool refineSrcBoxes = true;
|
||||||
|
while (refineSrcBoxes && (refinementIter <= nRefineIterMax_))
|
||||||
|
{
|
||||||
|
// Per processor refinement info
|
||||||
|
List<labelList> refineFlags(Pstream::nProcs());
|
||||||
|
labelList nElems(Pstream::nProcs(), 0);
|
||||||
|
|
||||||
|
// Assess how many target elements intersect the source bounding boxes
|
||||||
|
// and use the info to flag how the source boxes should be refined
|
||||||
|
setRefineFlags
|
||||||
|
(
|
||||||
|
refinementIter,
|
||||||
|
nTgtElems,
|
||||||
|
fixedSendElems,
|
||||||
|
localTgtElems,
|
||||||
|
refineFlags,
|
||||||
|
nElems
|
||||||
|
);
|
||||||
|
|
||||||
|
// Refine the source bounding boxes
|
||||||
|
refineSrcBoxes =
|
||||||
|
doRefineBoxes
|
||||||
|
(
|
||||||
|
refinementIter,
|
||||||
|
nSrcElems,
|
||||||
|
refineFlags,
|
||||||
|
nElems,
|
||||||
|
fixedBoxes
|
||||||
|
);
|
||||||
|
|
||||||
|
++refinementIter;
|
||||||
|
|
||||||
|
if (debug > 1)
|
||||||
|
{
|
||||||
|
// Include any boxes that are still evolving
|
||||||
|
List<DynamicList<treeBoundBox>> allBoxes(fixedBoxes);
|
||||||
|
forAll(allBoxes, proci)
|
||||||
|
{
|
||||||
|
allBoxes[proci].append(boxes_[proci]);
|
||||||
|
}
|
||||||
|
writeBoxes(allBoxes, refinementIter);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
if (debug)
|
||||||
|
{
|
||||||
|
Pout<< "Local src boxes after " << refinementIter-1 << " iterations:"
|
||||||
|
<< nl;
|
||||||
|
|
||||||
|
forAll(fixedBoxes, proci)
|
||||||
|
{
|
||||||
|
// Include any boxes that are still evolving in box count
|
||||||
|
label nBox = fixedBoxes[proci].size() + boxes_[proci].size();
|
||||||
|
Pout<< " proc:" << proci << " nBoxes:" << nBox << nl;
|
||||||
|
}
|
||||||
|
Pout<< endl;
|
||||||
|
}
|
||||||
|
|
||||||
|
// Convert send elems into a List<labelList>
|
||||||
|
List<labelList> sendElems(Pstream::nProcs());
|
||||||
|
forAll(localTgtElems, proci)
|
||||||
|
{
|
||||||
|
if (proci == Pstream::myProcNo() && nSrcElems)
|
||||||
|
{
|
||||||
|
sendElems[proci] = identity(nTgtElems);
|
||||||
|
}
|
||||||
|
else
|
||||||
|
{
|
||||||
|
labelHashSet& allIDs = fixedSendElems[proci];
|
||||||
|
|
||||||
|
const List<labelList>& procBoxElems = localTgtElems[proci];
|
||||||
|
|
||||||
|
for (const labelList& elems: procBoxElems)
|
||||||
|
{
|
||||||
|
allIDs.insert(elems);
|
||||||
|
}
|
||||||
|
|
||||||
|
sendElems[proci] = allIDs.toc();
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
fixedSendElems.clear();
|
||||||
|
localTgtElems.clear();
|
||||||
|
|
||||||
|
if (debug)
|
||||||
|
{
|
||||||
|
Pout<< "Local objects: " << nTgtElems << " Send map:" << nl
|
||||||
|
<< tab << "proc" << tab << "objects" << endl;
|
||||||
|
|
||||||
|
forAll(sendElems, proci)
|
||||||
|
{
|
||||||
|
Pout<< tab << proci << tab << sendElems[proci].size()
|
||||||
|
<< endl;
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
return createLODMap(sendElems);
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
Foam::autoPtr<Foam::mapDistribute> Foam::processorLODs::box::createLODMap
|
||||||
|
(
|
||||||
|
List<labelList>& sendElems
|
||||||
|
) const
|
||||||
|
{
|
||||||
|
// Send over how many objects I need to receive
|
||||||
|
const label localProci = Pstream::myProcNo();
|
||||||
|
labelListList sendSizes(Pstream::nProcs());
|
||||||
|
sendSizes[localProci].setSize(Pstream::nProcs());
|
||||||
|
forAll(sendElems, proci)
|
||||||
|
{
|
||||||
|
sendSizes[localProci][proci] = sendElems[proci].size();
|
||||||
|
}
|
||||||
|
Pstream::gatherList(sendSizes);
|
||||||
|
Pstream::scatterList(sendSizes);
|
||||||
|
|
||||||
|
|
||||||
|
// Determine order of receiving
|
||||||
|
labelListList constructMap(Pstream::nProcs());
|
||||||
|
|
||||||
|
// My local segment first
|
||||||
|
constructMap[localProci] = identity(sendElems[localProci].size());
|
||||||
|
|
||||||
|
label segmenti = constructMap[localProci].size();
|
||||||
|
forAll(constructMap, proci)
|
||||||
|
{
|
||||||
|
if (proci != localProci)
|
||||||
|
{
|
||||||
|
// What I need to receive is what other processor is sending to me
|
||||||
|
label nRecv = sendSizes[proci][localProci];
|
||||||
|
constructMap[proci].setSize(nRecv);
|
||||||
|
|
||||||
|
for (label& addr : constructMap[proci])
|
||||||
|
{
|
||||||
|
addr = segmenti++;
|
||||||
|
}
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
autoPtr<mapDistribute> mapPtr
|
||||||
|
(
|
||||||
|
new mapDistribute
|
||||||
|
(
|
||||||
|
segmenti, // size after construction
|
||||||
|
std::move(sendElems),
|
||||||
|
std::move(constructMap)
|
||||||
|
)
|
||||||
|
);
|
||||||
|
|
||||||
|
return mapPtr;
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
Foam::processorLODs::box::box
|
||||||
|
(
|
||||||
|
const UList<point>& srcPoints,
|
||||||
|
const UList<point>& tgtPoints,
|
||||||
|
const label maxObjectsPerLeaf,
|
||||||
|
const label nObjectsOfType,
|
||||||
|
const label nRefineIterMax
|
||||||
|
)
|
||||||
|
:
|
||||||
|
processorLOD(maxObjectsPerLeaf, nObjectsOfType),
|
||||||
|
srcPoints_(srcPoints),
|
||||||
|
tgtPoints_(tgtPoints),
|
||||||
|
boxes_(Pstream::nProcs()),
|
||||||
|
nRefineIterMax_(nRefineIterMax),
|
||||||
|
newToOld_(Pstream::nProcs()),
|
||||||
|
boxSrcElems_(Pstream::nProcs())
|
||||||
|
{
|
||||||
|
// Initialise each processor with a single box large enough to include all
|
||||||
|
// of its local src points
|
||||||
|
if (srcPoints_.size())
|
||||||
|
{
|
||||||
|
forAll(boxes_, proci)
|
||||||
|
{
|
||||||
|
List<treeBoundBox>& procBoxes = boxes_[proci];
|
||||||
|
|
||||||
|
// Note: inflate to ease overlap operations and to handle 2-D
|
||||||
|
// axis-aligned objects
|
||||||
|
treeBoundBox srcBb(srcPoints_);
|
||||||
|
srcBb.inflate(0.01);
|
||||||
|
|
||||||
|
DynamicList<treeBoundBox> newProcBoxes(1);
|
||||||
|
newProcBoxes.append(srcBb);
|
||||||
|
procBoxes.transfer(newProcBoxes);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
// ************************************************************************* //
|
||||||
206
src/meshTools/processorLOD/box/box.H
Normal file
206
src/meshTools/processorLOD/box/box.H
Normal file
@ -0,0 +1,206 @@
|
|||||||
|
/*---------------------------------------------------------------------------*\
|
||||||
|
========= |
|
||||||
|
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||||
|
\\ / O peration |
|
||||||
|
\\ / A nd | Copyright (C) 2017 OpeCFD Ltd.
|
||||||
|
\\/ M anipulation |
|
||||||
|
-------------------------------------------------------------------------------
|
||||||
|
License
|
||||||
|
This file is part of OpenFOAM.
|
||||||
|
|
||||||
|
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||||
|
under the terms of the GNU General Public License as published by
|
||||||
|
the Free Software Foundation, either version 3 of the License, or
|
||||||
|
(at your option) any later version.
|
||||||
|
|
||||||
|
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||||
|
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||||
|
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||||
|
for more details.
|
||||||
|
|
||||||
|
You should have received a copy of the GNU General Public License
|
||||||
|
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||||
|
|
||||||
|
Class
|
||||||
|
Foam::processorLODs::box
|
||||||
|
|
||||||
|
Description
|
||||||
|
Creates the parallel distribution map by describing the source and target
|
||||||
|
objects using box shapes.
|
||||||
|
|
||||||
|
A single box is created for the source object, which is then split using
|
||||||
|
2x2x2 refinement based on the number of remote target objects that overlap.
|
||||||
|
The refinement is local between communicating processor pairs, where the
|
||||||
|
refinement continues until a threshold number of remote target objects per
|
||||||
|
source box is achieved.
|
||||||
|
|
||||||
|
\*---------------------------------------------------------------------------*/
|
||||||
|
|
||||||
|
#ifndef processorLODs_box
|
||||||
|
#define processorLODs_box
|
||||||
|
|
||||||
|
#include "processorLOD.H"
|
||||||
|
#include "treeBoundBox.H"
|
||||||
|
#include "HashSet.H"
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
namespace Foam
|
||||||
|
{
|
||||||
|
|
||||||
|
namespace processorLODs
|
||||||
|
{
|
||||||
|
|
||||||
|
/*---------------------------------------------------------------------------*\
|
||||||
|
Class Box Declaration
|
||||||
|
\*---------------------------------------------------------------------------*/
|
||||||
|
|
||||||
|
class box
|
||||||
|
:
|
||||||
|
public processorLOD
|
||||||
|
{
|
||||||
|
protected:
|
||||||
|
|
||||||
|
// Protected data
|
||||||
|
|
||||||
|
// Flags to indicate what to do with a box
|
||||||
|
|
||||||
|
//- Drop/discard
|
||||||
|
static const label DROP;
|
||||||
|
|
||||||
|
//- Refine
|
||||||
|
static const label REFINE;
|
||||||
|
|
||||||
|
//- Fixed - do not touch
|
||||||
|
static const label FIXED;
|
||||||
|
|
||||||
|
|
||||||
|
//- Reference to the source points
|
||||||
|
const UList<point>& srcPoints_;
|
||||||
|
|
||||||
|
//- Reference to the target points
|
||||||
|
const UList<point>& tgtPoints_;
|
||||||
|
|
||||||
|
//- Per processor, the list of src bound boxes
|
||||||
|
// Note: using treeBoundBox to take advantage of subBbox() method
|
||||||
|
List<List<treeBoundBox>> boxes_;
|
||||||
|
|
||||||
|
//- Maximum number of refinement iterations
|
||||||
|
label nRefineIterMax_;
|
||||||
|
|
||||||
|
|
||||||
|
// Caching controls
|
||||||
|
|
||||||
|
//- Number of iterations before element indices are cached
|
||||||
|
static const label nStartUpIter;
|
||||||
|
|
||||||
|
//- Addressing per proc of new to old bound boxes
|
||||||
|
List<labelList> newToOld_;
|
||||||
|
|
||||||
|
//- Indices of elements in the src boxes
|
||||||
|
List<List<labelList>> boxSrcElems_;
|
||||||
|
|
||||||
|
|
||||||
|
// Private Member Functions
|
||||||
|
|
||||||
|
//- Helper function to write the boxes in OBJ format
|
||||||
|
void writeBoxes
|
||||||
|
(
|
||||||
|
const List<DynamicList<treeBoundBox>>& fixedBoxes,
|
||||||
|
const label iter
|
||||||
|
) const;
|
||||||
|
|
||||||
|
|
||||||
|
virtual boundBox calcSrcBox(const label srcObji) const = 0;
|
||||||
|
virtual boundBox calcTgtBox(const label tgtObji) const = 0;
|
||||||
|
|
||||||
|
//- Set the box refinement flags
|
||||||
|
void setRefineFlags
|
||||||
|
(
|
||||||
|
const label refineIter,
|
||||||
|
const label nTgtObjects,
|
||||||
|
List<labelHashSet>& fixedSendElems,
|
||||||
|
List<List<labelList>>& localTgtElems,
|
||||||
|
List<labelList>& refineFlags,
|
||||||
|
labelList& nElems
|
||||||
|
) const;
|
||||||
|
|
||||||
|
void refineBox
|
||||||
|
(
|
||||||
|
const label boxi,
|
||||||
|
const label refineIter,
|
||||||
|
const label nSrcElem,
|
||||||
|
const treeBoundBox& origBox,
|
||||||
|
DynamicList<treeBoundBox>& procBoxes,
|
||||||
|
DynamicList<labelList>& procBoxElems,
|
||||||
|
DynamicList<label>& procNewToOld
|
||||||
|
) const;
|
||||||
|
|
||||||
|
void refineBox
|
||||||
|
(
|
||||||
|
const label boxi,
|
||||||
|
const labelList& srcAddr,
|
||||||
|
const treeBoundBox& origBox,
|
||||||
|
DynamicList<treeBoundBox>& procBoxes,
|
||||||
|
DynamicList<labelList>& procBoxElems,
|
||||||
|
DynamicList<label>& procNewToOld
|
||||||
|
) const;
|
||||||
|
|
||||||
|
|
||||||
|
//- Apply the box refinements
|
||||||
|
// \return true if still refining
|
||||||
|
bool doRefineBoxes
|
||||||
|
(
|
||||||
|
const label refineIter,
|
||||||
|
const label nSrcFaces,
|
||||||
|
const List<labelList>& refineFlags,
|
||||||
|
const labelList& nElems,
|
||||||
|
List<DynamicList<treeBoundBox>>& fixedBoxes
|
||||||
|
);
|
||||||
|
|
||||||
|
autoPtr<mapDistribute> createMap
|
||||||
|
(
|
||||||
|
const label nSrcElems,
|
||||||
|
const label nTgtElems
|
||||||
|
);
|
||||||
|
|
||||||
|
//- Use the current list of send elements to create the mapDistribute
|
||||||
|
autoPtr<mapDistribute> createLODMap
|
||||||
|
(
|
||||||
|
List<labelList>& sendElems
|
||||||
|
) const;
|
||||||
|
|
||||||
|
|
||||||
|
public:
|
||||||
|
|
||||||
|
//- Runtime type information
|
||||||
|
TypeName("box");
|
||||||
|
|
||||||
|
//- Construct from list of points
|
||||||
|
box
|
||||||
|
(
|
||||||
|
const UList<point>& srcPoints,
|
||||||
|
const UList<point>& tgtPoints,
|
||||||
|
const label maxObjectsPerLeaf,
|
||||||
|
const label nObjectsOfType,
|
||||||
|
const label nRefineIterMax = 100
|
||||||
|
);
|
||||||
|
|
||||||
|
//- Destructor
|
||||||
|
virtual ~box() = default;
|
||||||
|
};
|
||||||
|
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
} // End namespace processorLODs
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
} // End namespace Foam
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
#endif
|
||||||
|
|
||||||
|
// ************************************************************************* //
|
||||||
110
src/meshTools/processorLOD/cellBox/cellBox.C
Normal file
110
src/meshTools/processorLOD/cellBox/cellBox.C
Normal file
@ -0,0 +1,110 @@
|
|||||||
|
/*---------------------------------------------------------------------------*\
|
||||||
|
========= |
|
||||||
|
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||||
|
\\ / O peration |
|
||||||
|
\\ / A nd | Copyright (C) 2017 OpenCFD Ltd.
|
||||||
|
\\/ M anipulation |
|
||||||
|
-------------------------------------------------------------------------------
|
||||||
|
License
|
||||||
|
This file is part of OpenFOAM.
|
||||||
|
|
||||||
|
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||||
|
under the terms of the GNU General Public License as published by
|
||||||
|
the Free Software Foundation, either version 3 of the License, or
|
||||||
|
(at your option) any later version.
|
||||||
|
|
||||||
|
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||||
|
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||||
|
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||||
|
for more details.
|
||||||
|
|
||||||
|
You should have received a copy of the GNU General Public License
|
||||||
|
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||||
|
|
||||||
|
\*---------------------------------------------------------------------------*/
|
||||||
|
|
||||||
|
#include "cellBox.H"
|
||||||
|
#include "mapDistribute.H"
|
||||||
|
|
||||||
|
namespace Foam
|
||||||
|
{
|
||||||
|
namespace processorLODs
|
||||||
|
{
|
||||||
|
defineTypeNameAndDebug(cellBox, 0);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * Protected Member Functions * * * * * * * * * * * //
|
||||||
|
|
||||||
|
Foam::boundBox Foam::processorLODs::cellBox::calcSrcBox
|
||||||
|
(
|
||||||
|
const label srcObji
|
||||||
|
) const
|
||||||
|
{
|
||||||
|
const UList<label>& cellFaces = srcCells_[srcObji];
|
||||||
|
|
||||||
|
boundBox bb(srcPoints_, srcFaces_[cellFaces[0]], false);
|
||||||
|
for (label i = 1; i < cellFaces.size(); ++i)
|
||||||
|
{
|
||||||
|
bb.add(srcPoints_, srcFaces_[cellFaces[i]]);
|
||||||
|
}
|
||||||
|
|
||||||
|
return bb;
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
Foam::boundBox Foam::processorLODs::cellBox::calcTgtBox
|
||||||
|
(
|
||||||
|
const label tgtObji
|
||||||
|
) const
|
||||||
|
{
|
||||||
|
const UList<label>& cellFaces = tgtCells_[tgtObji];
|
||||||
|
|
||||||
|
boundBox bb(tgtPoints_, tgtFaces_[cellFaces[0]], false);
|
||||||
|
for (label i = 1; i < cellFaces.size(); ++i)
|
||||||
|
{
|
||||||
|
bb.add(tgtPoints_, tgtFaces_[cellFaces[i]]);
|
||||||
|
}
|
||||||
|
|
||||||
|
return bb;
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
Foam::processorLODs::cellBox::cellBox
|
||||||
|
(
|
||||||
|
const cellList& srcCells,
|
||||||
|
const faceList& srcFaces,
|
||||||
|
const UList<point>& srcPoints,
|
||||||
|
const cellList& tgtCells,
|
||||||
|
const faceList& tgtFaces,
|
||||||
|
const UList<point>& tgtPoints,
|
||||||
|
const label maxObjectsPerLeaf,
|
||||||
|
const label nObjectsOfType,
|
||||||
|
const label nRefineIterMax
|
||||||
|
)
|
||||||
|
:
|
||||||
|
faceBox
|
||||||
|
(
|
||||||
|
srcFaces,
|
||||||
|
srcPoints,
|
||||||
|
tgtFaces,
|
||||||
|
tgtPoints,
|
||||||
|
maxObjectsPerLeaf,
|
||||||
|
nObjectsOfType,
|
||||||
|
nRefineIterMax
|
||||||
|
),
|
||||||
|
srcCells_(srcCells),
|
||||||
|
tgtCells_(tgtCells)
|
||||||
|
{}
|
||||||
|
|
||||||
|
|
||||||
|
Foam::autoPtr<Foam::mapDistribute> Foam::processorLODs::cellBox::map()
|
||||||
|
{
|
||||||
|
return createMap(srcCells_.size(), tgtCells_.size());
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
// ************************************************************************* //
|
||||||
120
src/meshTools/processorLOD/cellBox/cellBox.H
Normal file
120
src/meshTools/processorLOD/cellBox/cellBox.H
Normal file
@ -0,0 +1,120 @@
|
|||||||
|
/*---------------------------------------------------------------------------*\
|
||||||
|
========= |
|
||||||
|
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||||
|
\\ / O peration |
|
||||||
|
\\ / A nd | Copyright (C) 2017 OpeCFD Ltd.
|
||||||
|
\\/ M anipulation |
|
||||||
|
-------------------------------------------------------------------------------
|
||||||
|
License
|
||||||
|
This file is part of OpenFOAM.
|
||||||
|
|
||||||
|
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||||
|
under the terms of the GNU General Public License as published by
|
||||||
|
the Free Software Foundation, either version 3 of the License, or
|
||||||
|
(at your option) any later version.
|
||||||
|
|
||||||
|
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||||
|
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||||
|
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||||
|
for more details.
|
||||||
|
|
||||||
|
You should have received a copy of the GNU General Public License
|
||||||
|
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||||
|
|
||||||
|
Class
|
||||||
|
Foam::processorLODs::cellBox
|
||||||
|
|
||||||
|
Description
|
||||||
|
Creates the parallel distribution map by describing the source and target
|
||||||
|
objects using box shapes.
|
||||||
|
|
||||||
|
A single box is created for the source object, which is then split using
|
||||||
|
2x2x2 refinement based on the number of remote target objects that overlap.
|
||||||
|
The refinement is local between communicating processor pairs, where the
|
||||||
|
refinement continues until a threshold number of remote target objects per
|
||||||
|
source box is achieved.
|
||||||
|
|
||||||
|
\*---------------------------------------------------------------------------*/
|
||||||
|
|
||||||
|
#ifndef processorLODs_cellBox
|
||||||
|
#define processorLODs_cellBox
|
||||||
|
|
||||||
|
#include "faceBox.H"
|
||||||
|
#include "cellList.H"
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
namespace Foam
|
||||||
|
{
|
||||||
|
|
||||||
|
namespace processorLODs
|
||||||
|
{
|
||||||
|
|
||||||
|
/*---------------------------------------------------------------------------*\
|
||||||
|
Class cellBox Declaration
|
||||||
|
\*---------------------------------------------------------------------------*/
|
||||||
|
|
||||||
|
class cellBox
|
||||||
|
:
|
||||||
|
public faceBox
|
||||||
|
{
|
||||||
|
private:
|
||||||
|
|
||||||
|
// Private data
|
||||||
|
|
||||||
|
//- Reference to the source face list
|
||||||
|
const cellList& srcCells_;
|
||||||
|
|
||||||
|
//- Reference to the target face list
|
||||||
|
const cellList& tgtCells_;
|
||||||
|
|
||||||
|
|
||||||
|
// Private Member Functions
|
||||||
|
|
||||||
|
virtual boundBox calcSrcBox(const label srcObji) const;
|
||||||
|
virtual boundBox calcTgtBox(const label tgtObji) const;
|
||||||
|
|
||||||
|
|
||||||
|
public:
|
||||||
|
|
||||||
|
//- Runtime type information
|
||||||
|
TypeName("box");
|
||||||
|
|
||||||
|
//- Construct from list of points
|
||||||
|
cellBox
|
||||||
|
(
|
||||||
|
const cellList& srcCells,
|
||||||
|
const faceList& srcFaces,
|
||||||
|
const UList<point>& srcPoints,
|
||||||
|
const cellList& tgtCells,
|
||||||
|
const faceList& tgtFaces,
|
||||||
|
const UList<point>& tgtPoints,
|
||||||
|
const label maxObjectsPerLeaf,
|
||||||
|
const label nObjectsOfType,
|
||||||
|
const label nRefineIterMax = 100
|
||||||
|
);
|
||||||
|
|
||||||
|
//- Destructor
|
||||||
|
virtual ~cellBox() = default;
|
||||||
|
|
||||||
|
|
||||||
|
// Member Functions
|
||||||
|
|
||||||
|
//- Return the parallel distribution map
|
||||||
|
virtual autoPtr<mapDistribute> map();
|
||||||
|
};
|
||||||
|
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
} // End namespace processorLODs
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
} // End namespace Foam
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
#endif
|
||||||
|
|
||||||
|
// ************************************************************************* //
|
||||||
82
src/meshTools/processorLOD/faceBox/faceBox.C
Normal file
82
src/meshTools/processorLOD/faceBox/faceBox.C
Normal file
@ -0,0 +1,82 @@
|
|||||||
|
/*---------------------------------------------------------------------------*\
|
||||||
|
========= |
|
||||||
|
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||||
|
\\ / O peration |
|
||||||
|
\\ / A nd | Copyright (C) 2017 OpenCFD Ltd.
|
||||||
|
\\/ M anipulation |
|
||||||
|
-------------------------------------------------------------------------------
|
||||||
|
License
|
||||||
|
This file is part of OpenFOAM.
|
||||||
|
|
||||||
|
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||||
|
under the terms of the GNU General Public License as published by
|
||||||
|
the Free Software Foundation, either version 3 of the License, or
|
||||||
|
(at your option) any later version.
|
||||||
|
|
||||||
|
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||||
|
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||||
|
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||||
|
for more details.
|
||||||
|
|
||||||
|
You should have received a copy of the GNU General Public License
|
||||||
|
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||||
|
|
||||||
|
\*---------------------------------------------------------------------------*/
|
||||||
|
|
||||||
|
#include "faceBox.H"
|
||||||
|
#include "mapDistribute.H"
|
||||||
|
|
||||||
|
namespace Foam
|
||||||
|
{
|
||||||
|
namespace processorLODs
|
||||||
|
{
|
||||||
|
defineTypeNameAndDebug(faceBox, 0);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * Protected Member Functions * * * * * * * * * * * //
|
||||||
|
|
||||||
|
Foam::boundBox Foam::processorLODs::faceBox::calcSrcBox
|
||||||
|
(
|
||||||
|
const label srcObji
|
||||||
|
) const
|
||||||
|
{
|
||||||
|
return boundBox(srcPoints_, srcFaces_[srcObji], false);
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
Foam::boundBox Foam::processorLODs::faceBox::calcTgtBox
|
||||||
|
(
|
||||||
|
const label tgtObji
|
||||||
|
) const
|
||||||
|
{
|
||||||
|
return boundBox(tgtPoints_, tgtFaces_[tgtObji], false);
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
Foam::processorLODs::faceBox::faceBox
|
||||||
|
(
|
||||||
|
const faceList& srcFaces,
|
||||||
|
const UList<point>& srcPoints,
|
||||||
|
const faceList& tgtFaces,
|
||||||
|
const UList<point>& tgtPoints,
|
||||||
|
const label maxObjectsPerLeaf,
|
||||||
|
const label nObjectsOfType,
|
||||||
|
const label nRefineIterMax
|
||||||
|
)
|
||||||
|
:
|
||||||
|
box(srcPoints, tgtPoints, maxObjectsPerLeaf, nObjectsOfType),
|
||||||
|
srcFaces_(srcFaces),
|
||||||
|
tgtFaces_(tgtFaces)
|
||||||
|
{}
|
||||||
|
|
||||||
|
|
||||||
|
Foam::autoPtr<Foam::mapDistribute> Foam::processorLODs::faceBox::map()
|
||||||
|
{
|
||||||
|
return createMap(srcFaces_.size(), tgtFaces_.size());
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
// ************************************************************************* //
|
||||||
117
src/meshTools/processorLOD/faceBox/faceBox.H
Normal file
117
src/meshTools/processorLOD/faceBox/faceBox.H
Normal file
@ -0,0 +1,117 @@
|
|||||||
|
/*---------------------------------------------------------------------------*\
|
||||||
|
========= |
|
||||||
|
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||||
|
\\ / O peration |
|
||||||
|
\\ / A nd | Copyright (C) 2017 OpeCFD Ltd.
|
||||||
|
\\/ M anipulation |
|
||||||
|
-------------------------------------------------------------------------------
|
||||||
|
License
|
||||||
|
This file is part of OpenFOAM.
|
||||||
|
|
||||||
|
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||||
|
under the terms of the GNU General Public License as published by
|
||||||
|
the Free Software Foundation, either version 3 of the License, or
|
||||||
|
(at your option) any later version.
|
||||||
|
|
||||||
|
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||||
|
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||||
|
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||||
|
for more details.
|
||||||
|
|
||||||
|
You should have received a copy of the GNU General Public License
|
||||||
|
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||||
|
|
||||||
|
Class
|
||||||
|
Foam::processorLODs::faceBox
|
||||||
|
|
||||||
|
Description
|
||||||
|
Creates the parallel distribution map by describing the source and target
|
||||||
|
objects using box shapes.
|
||||||
|
|
||||||
|
A single box is created for the source object, which is then split using
|
||||||
|
2x2x2 refinement based on the number of remote target objects that overlap.
|
||||||
|
The refinement is local between communicating processor pairs, where the
|
||||||
|
refinement continues until a threshold number of remote target objects per
|
||||||
|
source box is achieved.
|
||||||
|
|
||||||
|
\*---------------------------------------------------------------------------*/
|
||||||
|
|
||||||
|
#ifndef processorLODs_faceBox
|
||||||
|
#define processorLODs_faceBox
|
||||||
|
|
||||||
|
#include "box.H"
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
namespace Foam
|
||||||
|
{
|
||||||
|
|
||||||
|
namespace processorLODs
|
||||||
|
{
|
||||||
|
|
||||||
|
/*---------------------------------------------------------------------------*\
|
||||||
|
Class faceBox Declaration
|
||||||
|
\*---------------------------------------------------------------------------*/
|
||||||
|
|
||||||
|
class faceBox
|
||||||
|
:
|
||||||
|
public box
|
||||||
|
{
|
||||||
|
protected:
|
||||||
|
|
||||||
|
// Protected data
|
||||||
|
|
||||||
|
//- Reference to the source face list
|
||||||
|
const faceList& srcFaces_;
|
||||||
|
|
||||||
|
//- Reference to the target face list
|
||||||
|
const faceList& tgtFaces_;
|
||||||
|
|
||||||
|
|
||||||
|
// Protected Member Functions
|
||||||
|
|
||||||
|
virtual boundBox calcSrcBox(const label srcObji) const;
|
||||||
|
virtual boundBox calcTgtBox(const label tgtObji) const;
|
||||||
|
|
||||||
|
|
||||||
|
public:
|
||||||
|
|
||||||
|
//- Runtime type information
|
||||||
|
TypeName("box");
|
||||||
|
|
||||||
|
//- Construct from list of points
|
||||||
|
faceBox
|
||||||
|
(
|
||||||
|
const faceList& srcFaces,
|
||||||
|
const UList<point>& srcPoints,
|
||||||
|
const faceList& tgtFaces,
|
||||||
|
const UList<point>& tgtPoints,
|
||||||
|
const label maxObjectsPerLeaf,
|
||||||
|
const label nObjectsOfType,
|
||||||
|
const label nRefineIterMax = 100
|
||||||
|
);
|
||||||
|
|
||||||
|
//- Destructor
|
||||||
|
virtual ~faceBox() = default;
|
||||||
|
|
||||||
|
|
||||||
|
// Member Functions
|
||||||
|
|
||||||
|
//- Return the parallel distribution map
|
||||||
|
virtual autoPtr<mapDistribute> map();
|
||||||
|
};
|
||||||
|
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
} // End namespace processorLODs
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
} // End namespace Foam
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
#endif
|
||||||
|
|
||||||
|
// ************************************************************************* //
|
||||||
46
src/meshTools/processorLOD/processorLOD/processorLOD.C
Normal file
46
src/meshTools/processorLOD/processorLOD/processorLOD.C
Normal file
@ -0,0 +1,46 @@
|
|||||||
|
/*---------------------------------------------------------------------------*\
|
||||||
|
========= |
|
||||||
|
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||||
|
\\ / O peration |
|
||||||
|
\\ / A nd | Copyright (C) 2017 OpenCFD Ltd.
|
||||||
|
\\/ M anipulation |
|
||||||
|
-------------------------------------------------------------------------------
|
||||||
|
License
|
||||||
|
This file is part of OpenFOAM.
|
||||||
|
|
||||||
|
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||||
|
under the terms of the GNU General Public License as published by
|
||||||
|
the Free Software Foundation, either version 3 of the License, or
|
||||||
|
(at your option) any later version.
|
||||||
|
|
||||||
|
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||||
|
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||||
|
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||||
|
for more details.
|
||||||
|
|
||||||
|
You should have received a copy of the GNU General Public License
|
||||||
|
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||||
|
|
||||||
|
\*---------------------------------------------------------------------------*/
|
||||||
|
|
||||||
|
#include "processorLOD.H"
|
||||||
|
|
||||||
|
namespace Foam
|
||||||
|
{
|
||||||
|
defineTypeNameAndDebug(processorLOD, 0);
|
||||||
|
}
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
Foam::processorLOD::processorLOD
|
||||||
|
(
|
||||||
|
const label maxObjectsPerLeaf,
|
||||||
|
const label nObjectsOfType
|
||||||
|
)
|
||||||
|
:
|
||||||
|
maxObjectsPerLeaf_(maxObjectsPerLeaf),
|
||||||
|
nObjectsOfType_(nObjectsOfType)
|
||||||
|
{}
|
||||||
|
|
||||||
|
|
||||||
|
// ************************************************************************* //
|
||||||
97
src/meshTools/processorLOD/processorLOD/processorLOD.H
Normal file
97
src/meshTools/processorLOD/processorLOD/processorLOD.H
Normal file
@ -0,0 +1,97 @@
|
|||||||
|
/*---------------------------------------------------------------------------*\
|
||||||
|
========= |
|
||||||
|
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||||
|
\\ / O peration |
|
||||||
|
\\ / A nd | Copyright (C) 2017 OpenCFD Ltd.
|
||||||
|
\\/ M anipulation |
|
||||||
|
-------------------------------------------------------------------------------
|
||||||
|
License
|
||||||
|
This file is part of OpenFOAM.
|
||||||
|
|
||||||
|
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||||
|
under the terms of the GNU General Public License as published by
|
||||||
|
the Free Software Foundation, either version 3 of the License, or
|
||||||
|
(at your option) any later version.
|
||||||
|
|
||||||
|
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||||
|
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||||
|
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||||
|
for more details.
|
||||||
|
|
||||||
|
You should have received a copy of the GNU General Public License
|
||||||
|
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||||
|
|
||||||
|
Class
|
||||||
|
Foam::processorLOD
|
||||||
|
|
||||||
|
Description
|
||||||
|
Base class to generate a parallel distribution map for sending sufficient
|
||||||
|
target objects to cover a description of the source object, based on
|
||||||
|
processor Level Of Detail (LOD) shapes
|
||||||
|
|
||||||
|
\*---------------------------------------------------------------------------*/
|
||||||
|
|
||||||
|
#ifndef processorLOD_H
|
||||||
|
#define processorLOD_H
|
||||||
|
|
||||||
|
#include "autoPtr.H"
|
||||||
|
#include "typeInfo.H"
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
namespace Foam
|
||||||
|
{
|
||||||
|
|
||||||
|
class mapDistribute;
|
||||||
|
|
||||||
|
/*---------------------------------------------------------------------------*\
|
||||||
|
Class processorLOD Declaration
|
||||||
|
\*---------------------------------------------------------------------------*/
|
||||||
|
|
||||||
|
class processorLOD
|
||||||
|
{
|
||||||
|
|
||||||
|
protected:
|
||||||
|
|
||||||
|
// Protected data
|
||||||
|
|
||||||
|
//- Maximum number of objects per leaf
|
||||||
|
label maxObjectsPerLeaf_;
|
||||||
|
|
||||||
|
//- Number of objects of this type, e.g. number of faces/cells on this
|
||||||
|
//- processor
|
||||||
|
label nObjectsOfType_;
|
||||||
|
|
||||||
|
|
||||||
|
public:
|
||||||
|
|
||||||
|
//- Runtime type information
|
||||||
|
TypeName("processorLOD");
|
||||||
|
|
||||||
|
//- Construct from components
|
||||||
|
processorLOD
|
||||||
|
(
|
||||||
|
const label maxObjectsPerLeaf,
|
||||||
|
const label nObjectsOfType
|
||||||
|
);
|
||||||
|
|
||||||
|
//- Destructor
|
||||||
|
virtual ~processorLOD() = default;
|
||||||
|
|
||||||
|
|
||||||
|
// Member Functions
|
||||||
|
|
||||||
|
//- Return the parallel distribution map
|
||||||
|
virtual autoPtr<mapDistribute> map() = 0;
|
||||||
|
};
|
||||||
|
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
} // End namespace Foam
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
#endif
|
||||||
|
|
||||||
|
// ************************************************************************* //
|
||||||
Reference in New Issue
Block a user