Skip to content
Snippets Groups Projects
Commit e8254972 authored by mattijs's avatar mattijs
Browse files

ENH: GAMGAgglomeration: added profile printing

parent 2496b4ed
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
src/OpenFOAM/matrices/lduMatrix/solvers/GAMG/GAMGAgglomerations/GAMGAgglomeration/GAMGAgglomeration.C
+ 113
− 92
View file @ e8254972
...@@ -69,111 +69,132 @@ void Foam::GAMGAgglomeration::compactLevels(const label nCreatedLevels) ...@@ -69,111 +69,132 @@ void Foam::GAMGAgglomeration::compactLevels(const label nCreatedLevels)
procAgglomeratorPtr_().agglomerate(); procAgglomeratorPtr_().agglomerate();
if (debug) }
{
Info<< "GAMGAgglomeration:" << nl // Print a bit
<< " local agglomerator : " << type() << nl if (processorAgglomerate() && debug)
<< " processor agglomerator : " {
Info<< "GAMGAgglomeration:" << nl
<< " local agglomerator : " << type() << nl;
if (processorAgglomerate())
{
Info<< " processor agglomerator : "
<< procAgglomeratorPtr_().type() << nl << procAgglomeratorPtr_().type() << nl
<< nl; << nl;
}
Info<< setw(40) << "nCells" Info<< setw(36) << "nCells"
<< setw(24) << "nFaces/nCells" << setw(20) << "nFaces/nCells"
<< setw(24) << "nInterfaces" << setw(20) << "nInterfaces"
<< setw(24) << "nIntFaces/nCells" << nl << setw(20) << "nIntFaces/nCells"
<< setw(8) << "Level" << setw(12) << "profile"
<< setw(8) << "nProcs" << nl
<< " " << setw(8) << "Level"
<< setw(8) << "avg" << setw(8) << "nProcs"
<< setw(8) << "max" << " "
<< " " << setw(8) << "avg"
<< setw(8) << "avg" << setw(8) << "max"
<< setw(8) << "max" << " "
<< " " << setw(8) << "avg"
<< setw(8) << "avg" << setw(8) << "max"
<< setw(8) << "max" << " "
<< " " << setw(4) << "avg" << setw(8) << "avg"
<< " " << setw(4) << "max" << setw(8) << "max"
<< nl << " "
<< setw(8) << "-----" << setw(8) << "avg"
<< setw(8) << "------" << setw(8) << "max"
<< " " //<< " "
<< setw(8) << "---" << setw(12) << "avg"
<< setw(8) << "---" << nl
<< " " << setw(8) << "-----"
<< setw(8) << "---" << setw(8) << "------"
<< setw(8) << "---" << " "
<< " " << setw(4) << "---" << setw(8) << "---"
<< " " << setw(4) << "---" << setw(8) << "---"
<< nl; << " "
<< setw(8) << "---"
for (label levelI = 0; levelI <= size(); levelI++) << setw(8) << "---"
<< " "
<< setw(8) << "---"
<< setw(8) << "---"
<< " "
<< setw(8) << "---"
<< setw(8) << "---"
//<< " "
<< setw(12) << "---"
//<< " "
<< nl;
for (label levelI = 0; levelI <= size(); levelI++)
{
label nProcs = 0;
label nCells = 0;
scalar faceCellRatio = 0;
label nInterfaces = 0;
label nIntFaces = 0;
scalar ratio = 0.0;
scalar profile = 0.0;
if (hasMeshLevel(levelI))
{ {
label nProcs = 0; nProcs = 1;
label nCells = 0;
scalar faceCellRatio = 0;
label nInterfaces = 0;
label nIntFaces = 0;
scalar ratio = 0.0;
if (hasMeshLevel(levelI))
{
nProcs = 1;
const lduMesh& fineMesh = meshLevel(levelI); const lduMesh& fineMesh = meshLevel(levelI);
nCells = fineMesh.lduAddr().size(); nCells = fineMesh.lduAddr().size();
faceCellRatio = faceCellRatio =
scalar(fineMesh.lduAddr().lowerAddr().size())/nCells; scalar(fineMesh.lduAddr().lowerAddr().size())/nCells;
const lduInterfacePtrsList interfaces = const lduInterfacePtrsList interfaces =
fineMesh.interfaces(); fineMesh.interfaces();
forAll(interfaces, i) forAll(interfaces, i)
{
if (interfaces.set(i))
{ {
if (interfaces.set(i)) nInterfaces++;
{ nIntFaces += interfaces[i].faceCells().size();
nInterfaces++;
nIntFaces += interfaces[i].faceCells().size();
}
} }
ratio = scalar(nIntFaces)/nCells;
} }
ratio = scalar(nIntFaces)/nCells;
label totNprocs = returnReduce(nProcs, sumOp<label>()); profile = fineMesh.lduAddr().band().second();
label maxNCells = returnReduce(nCells, maxOp<label>());
label totNCells = returnReduce(nCells, sumOp<label>());
scalar maxFaceCellRatio =
returnReduce(faceCellRatio, maxOp<scalar>());
scalar totFaceCellRatio =
returnReduce(faceCellRatio, sumOp<scalar>());
label maxNInt = returnReduce(nInterfaces, maxOp<label>());
label totNInt = returnReduce(nInterfaces, sumOp<label>());
scalar maxRatio = returnReduce(ratio, maxOp<scalar>());
scalar totRatio = returnReduce(ratio, sumOp<scalar>());
Info<< setw(8) << levelI
<< setw(8) << totNprocs
<< setw(16) << totNCells/totNprocs
<< setw(8) << maxNCells
<< " "
<< setw(8) << setprecision(4) << totFaceCellRatio/totNprocs
<< setw(8) << setprecision(4) << maxFaceCellRatio
<< " "
<< setw(8) << scalar(totNInt)/totNprocs
<< setw(8) << maxNInt
<< " "
<< setw(8) << setprecision(4) << totRatio/totNprocs
<< setw(8) << setprecision(4) << maxRatio
<< nl;
} }
Info<< endl;
label totNprocs = returnReduce(nProcs, sumOp<label>());
label maxNCells = returnReduce(nCells, maxOp<label>());
label totNCells = returnReduce(nCells, sumOp<label>());
scalar maxFaceCellRatio =
returnReduce(faceCellRatio, maxOp<scalar>());
scalar totFaceCellRatio =
returnReduce(faceCellRatio, sumOp<scalar>());
label maxNInt = returnReduce(nInterfaces, maxOp<label>());
label totNInt = returnReduce(nInterfaces, sumOp<label>());
scalar maxRatio = returnReduce(ratio, maxOp<scalar>());
scalar totRatio = returnReduce(ratio, sumOp<scalar>());
scalar totProfile = returnReduce(profile, sumOp<scalar>());
Info<< setw(8) << levelI
<< setw(8) << totNprocs
<< " "
<< setw(8) << totNCells/totNprocs
<< setw(8) << maxNCells
<< " "
<< setw(8) << setprecision(4) << totFaceCellRatio/totNprocs
<< setw(8) << setprecision(4) << maxFaceCellRatio
<< " "
<< setw(8) << scalar(totNInt)/totNprocs
<< setw(8) << maxNInt
<< " "
<< setw(8) << setprecision(4) << totRatio/totNprocs
<< setw(8) << setprecision(4) << maxRatio
<< setw(12) << setprecision(4) << totProfile/totNprocs
<< nl;
} }
Info<< endl;
} }
} }
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment