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: Test-GAMGAgglomeration: additional options

Browse Source
This commit is contained in:
mattijs
2011-11-08 12:56:33 +00:00
parent 063d3ca9a4
commit 827cc609d7
2 changed files with 58 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
applications/test/GAMGAgglomeration/Make/options
View File

@ -1,5 +1,7 @@
EXE_INC = \
-I$(LIB_SRC)/finiteVolume/lnInclude
-I$(LIB_SRC)/finiteVolume/lnInclude \
-I$(LIB_SRC)/meshTools/lnInclude
EXE_LIBS = \
-lfiniteVolume
-lfiniteVolume \
-lmeshTools

55
applications/test/GAMGAgglomeration/Test-GAMGAgglomeration.C
View File

@ -31,14 +31,31 @@ Description
#include "fvCFD.H"
#include "GAMGAgglomeration.H"
#include "OFstream.H"
#include "meshTools.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
// Main program:
int main(int argc, char *argv[])
{
argList::addBoolOption
(
"writeObj",
"write obj files of agglomeration"
);
argList::addBoolOption
(
"normalise",
"normalise agglomeration (0..1)"
);
#include "setRootCase.H"
#include "createTime.H"
bool writeObj = args.optionFound("writeObj");
bool normalise = args.optionFound("normalise");
#include "createMesh.H"
const fvSolution& sol = static_cast<const fvSolution&>(mesh);
@ -130,11 +147,47 @@ int main(int argc, char *argv[])
{
fld[cellI] = cellToCoarse[cellI];
}
fld /= max(fld);
if (normalise)
{
fld /= max(fld);
}
scalarAgglomeration.correctBoundaryConditions();
scalarAgglomeration.write();
}
if (writeObj)
{
OFstream str(runTime.path()/runTime.timeName()/"aggomeration.obj");
label vertI = 0;
// Write all mesh cc
forAll(mesh.cellCentres(), cellI)
{
meshTools::writeOBJ(str, mesh.cellCentres()[cellI]);
vertI++;
}
// Determine coarse cc
forAll(coarseToCell, coarseI)
{
const labelList& cellLabels = coarseToCell[coarseI];
point coarseCc = average
(
pointField(mesh.cellCentres(), cellLabels)
);
meshTools::writeOBJ(str, coarseCc);
vertI++;
forAll(cellLabels, i)
{
label cellI = cellLabels[i];
str << "l " << cellI+1 << ' ' << vertI << nl;
}
}
}
Info<< endl;
}