exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Better debug data for ggi and regionCouple

Browse source
This commit is contained in:
Hrvoje Jasak 2019-10-07 11:13:51 +01:00
parent 18de725d8f
commit c8cd512a24
2 changed files with 85 additions and 26 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

29
src/foam/meshes/polyMesh/polyPatches/constraint/ggi/ggiPolyPatch.C
View file

@ -1113,7 +1113,7 @@ void Foam::ggiPolyPatch::updateMesh()
void Foam::ggiPolyPatch::calcTransforms() const
{
// Simplest mixing plane: no transform or separation. HJ, 24/Oct/2008
// Simplest interface: no transform or separation. HJ, 24/Oct/2008
forwardT_.setSize(0);
reverseT_.setSize(0);
separation_.setSize(0);
@ -1167,6 +1167,7 @@ void Foam::ggiPolyPatch::calcTransforms() const
<< patchToPatch().uncoveredSlaveFaces().size()
<< " shadow().localPoints(): " << shadow().localPoints().size()
<< endl;
indirectPrimitivePatch::writeVTK
(
fvPath/fileName("uncoveredGgiFaces" + shadowName()),
@ -1180,31 +1181,7 @@ void Foam::ggiPolyPatch::calcTransforms() const
}
// Check for bridge overlap
if (!bridgeOverlap())
{
if
(
(
patchToPatch().uncoveredMasterFaces().size() > 0
&& !empty()
)
|| (
!shadow().empty()
&& patchToPatch().uncoveredSlaveFaces().size() > 0
)
)
{
FatalErrorIn("label ggiPolyPatch::calcTransforms() const")
<< "ggi patch " << name() << " with shadow "
<< shadowName() << " has "
<< patchToPatch().uncoveredMasterFaces().size()
<< " uncovered master faces and "
<< patchToPatch().uncoveredSlaveFaces().size()
<< " uncovered slave faces. Bridging is switched off. "
<< abort(FatalError);
}
}
else
if (bridgeOverlap())
{
InfoIn("label ggiPolyPatch::calcTransforms() const")
<< "ggi patch " << name() << " with shadow "

82
src/foam/meshes/polyMesh/polyPatches/constraint/regionCouple/regionCouplePolyPatch.C
View file

@ -35,6 +35,7 @@ Author
#include "polyPatchID.H"
#include "ZoneIDs.H"
#include "objectRegistry.H"
#include "indirectPrimitivePatch.H"
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
@ -1207,6 +1208,87 @@ void Foam::regionCouplePolyPatch::calcTransforms() const
forwardT_.setSize(0);
reverseT_.setSize(0);
separation_.setSize(0);
if (debug > 1 && master())
{
if
(
!empty()
&& patchToPatch().uncoveredMasterFaces().size() > 0
)
{
// Write uncovered master faces
Info<< "Writing uncovered master faces for patch "
<< name() << " as VTK." << endl;
const polyMesh& mesh = boundaryMesh().mesh();
fileName fvPath(mesh.time().path()/"VTK");
mkDir(fvPath);
indirectPrimitivePatch::writeVTK
(
fvPath/fileName("uncoveredRegionCoupleFaces" + name()),
IndirectList<face>
(
localFaces(),
patchToPatch().uncoveredMasterFaces()
),
localPoints()
);
}
if
(
!shadow().empty()
&& patchToPatch().uncoveredSlaveFaces().size() > 0
)
{
// Write uncovered master faces
Info<< "Writing uncovered shadow faces for patch "
<< shadowPatchName() << " as VTK." << endl;
const polyMesh& mesh = boundaryMesh().mesh();
fileName fvPath(mesh.time().path()/"VTK");
mkDir(fvPath);
Pout<< "Patch " << name()
<< " shadow().localFaces(): " << shadow().localFaces().size()
<< " patchToPatch().uncoveredSlaveFaces().size(): "
<< patchToPatch().uncoveredSlaveFaces().size()
<< " shadow().localPoints(): " << shadow().localPoints().size()
<< endl;
indirectPrimitivePatch::writeVTK
(
fvPath/fileName
(
"uncoveredRegionCoupleFaces"
+ shadowPatchName()
),
IndirectList<face>
(
shadow().localFaces(),
patchToPatch().uncoveredSlaveFaces()
),
shadow().localPoints()
);
}
// Check for bridge overlap
if (bridgeOverlap())
{
InfoInFunction
<< "regionCouple patch " << name() << " with shadow "
<< shadowPatchName() << " on region " << shadowRegionName()
<< " has "
<< patchToPatch().uncoveredMasterFaces().size()
<< " uncovered master faces and "
<< patchToPatch().uncoveredSlaveFaces().size()
<< " uncovered slave faces. Bridging is switched on. "
<< endl;
}
}
}