fc0b980333
- since bounding boxes overlap, need to verify the splitting actually did something. Problem only really evident with higher tree depths.
158 lines
4.2 KiB
C
158 lines
4.2 KiB
C
/*---------------------------------------------------------------------------*\
|
|
========= |
|
|
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
|
\\ / O peration |
|
|
\\ / A nd | www.openfoam.com
|
|
\\/ M anipulation |
|
|
-------------------------------------------------------------------------------
|
|
Copyright (C) 2022 OpenCFD Ltd.
|
|
-------------------------------------------------------------------------------
|
|
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/>.
|
|
|
|
Application
|
|
Test-surfaceTree
|
|
|
|
Description
|
|
Simple tests for building indexedOctree etc.
|
|
|
|
\*---------------------------------------------------------------------------*/
|
|
|
|
#include "argList.H"
|
|
#include "clockTime.H"
|
|
#include "MeshedSurfaces.H"
|
|
#include "indexedOctree.H"
|
|
#include "AABBTree.H"
|
|
#include "treeDataPrimitivePatch.H"
|
|
#include "Random.H"
|
|
#include "ListListOps.H"
|
|
|
|
using namespace Foam;
|
|
|
|
|
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
|
// Main program:
|
|
|
|
int main(int argc, char *argv[])
|
|
{
|
|
argList::addNote
|
|
(
|
|
"Test octree building etc"
|
|
);
|
|
|
|
argList::noBanner();
|
|
argList::noParallel();
|
|
argList::addArgument("input", "The input surface file");
|
|
|
|
argList::addOption("maxLevel", "int", "The max level");
|
|
argList::addOption("leafSize", "int", "The min leaf size");
|
|
argList::addBoolOption("AABB", "AABBTree instead of indexedOctree");
|
|
|
|
argList args(argc, argv);
|
|
|
|
const auto importName = args.get<fileName>(1);
|
|
|
|
const bool useAABB = args.found("AABB");
|
|
|
|
label maxLevel = 10;
|
|
label leafSize = 10;
|
|
|
|
if (useAABB)
|
|
{
|
|
// May want different settings..
|
|
}
|
|
|
|
args.readIfPresent("maxLevel", maxLevel);
|
|
args.readIfPresent("leafSize", leafSize);
|
|
|
|
meshedSurface surf(importName);
|
|
|
|
Random rndGen(123456);
|
|
|
|
treeBoundBox overallBb
|
|
(
|
|
treeBoundBox(surf.box()).extend(rndGen, 1e-4, ROOTVSMALL)
|
|
);
|
|
|
|
Info<< "Surface with " << surf.size() << " faces, "
|
|
<< surf.nPoints() << " points" << endl;
|
|
|
|
|
|
clockTime timing;
|
|
|
|
if (useAABB)
|
|
{
|
|
AABBTree<face> tree
|
|
(
|
|
surf,
|
|
surf.points(),
|
|
true, // Equal bins
|
|
maxLevel, // maxLevel
|
|
leafSize // minLeafSize
|
|
);
|
|
|
|
Info<< "Built AABBTree, maxLevel:" << maxLevel
|
|
<< " minLeaf:" << leafSize
|
|
<< " with " << tree.boundBoxes().size()
|
|
<< " leaves - " << timing.elapsedTime() << 's' << endl;
|
|
|
|
{
|
|
OFstream os("AABBTree.obj");
|
|
tree.writeOBJ(os);
|
|
|
|
Info<< "Wrote " << os.name() << endl;
|
|
}
|
|
|
|
// Info<< "sizes: ";
|
|
// ListListOps::subSizes
|
|
// (
|
|
// tree.addressing(),
|
|
// identityOp{}
|
|
// ).writeList(Info) << endl;
|
|
}
|
|
else
|
|
{
|
|
indexedOctree<treeDataPrimitivePatch<meshedSurface>> tree
|
|
(
|
|
treeDataPrimitivePatch<meshedSurface>(surf, 1e-6),
|
|
overallBb,
|
|
maxLevel, // maxLevel
|
|
leafSize, // leafSize
|
|
3.0 // duplicity
|
|
);
|
|
|
|
Info<< "Built octree, maxLevel:" << maxLevel
|
|
<< " minLeaf:" << leafSize
|
|
<< " with " << tree.nodes().size()
|
|
<< " nodes, " << tree.nLeafs()
|
|
<< " leaves - " << timing.elapsedTime() << 's' << endl;
|
|
|
|
{
|
|
OFstream os("indexedOctree.obj");
|
|
tree.writeOBJ(os);
|
|
|
|
Info<< "Wrote " << os.name() << endl;
|
|
}
|
|
}
|
|
|
|
|
|
Info<< "\nEnd\n" << endl;
|
|
return 0;
|
|
}
|
|
|
|
|
|
// ************************************************************************* //
|