From 9449d92fb9486f430210e568a7d90a3508f72fb2 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Thu, 7 Nov 2019 11:03:47 +0000
Subject: [PATCH] STYLE: shm: improved debug messages
---
.../medialAxisMeshMover.C | 42 +++++++++++++++----
.../meshRefinement/meshRefinementBaffles.C | 2 +-
2 files changed, 35 insertions(+), 9 deletions(-)
diff --git a/src/mesh/snappyHexMesh/externalDisplacementMeshMover/medialAxisMeshMover.C b/src/mesh/snappyHexMesh/externalDisplacementMeshMover/medialAxisMeshMover.C
index a70742cbf83..b5ce64a5f78 100644
--- a/src/mesh/snappyHexMesh/externalDisplacementMeshMover/medialAxisMeshMover.C
+++ b/src/mesh/snappyHexMesh/externalDisplacementMeshMover/medialAxisMeshMover.C
@@ -36,6 +36,7 @@ License
#include "OBJstream.H"
#include "PointData.H"
#include "zeroFixedValuePointPatchFields.H"
+#include "pointSet.H"
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
@@ -568,6 +569,33 @@ void Foam::medialAxisMeshMover::update(const dictionary& coeffDict)
maxInfo.shrink();
maxPoints.shrink();
+
+ if (debug)
+ {
+ mkDir(mesh().time().timePath());
+ OBJstream str(mesh().time().timePath()/"medialSurfacePoints.obj");
+
+ pointSet seedPoints
+ (
+ mesh(),
+ "medialSurfacePoints",
+ maxPoints
+ );
+
+ Info<< typeName
+ << " : Writing estimated medial surface:" << nl << incrIndent
+ << indent << "locations : " << str.name() << nl
+ << indent << "pointSet : " << seedPoints.name() << nl
+ << decrIndent << endl;
+
+ for (const auto& info : maxInfo)
+ {
+ str.write(info.origin());
+ }
+ seedPoints.write();
+ }
+
+
// Do all calculations
PointEdgeWave<pointEdgePoint> medialDistCalc
(
@@ -658,18 +686,16 @@ void Foam::medialAxisMeshMover::update(const dictionary& coeffDict)
if (debug)
{
Info<< typeName
- << " : Writing medial axis fields:" << nl
- << incrIndent
- << "ratio of medial distance to wall distance : "
+ << " : Writing medial axis fields:" << nl << incrIndent
+ << indent << "ratio of medial distance to wall distance : "
<< medialRatio_.name() << nl
- << "distance to nearest medial axis : "
+ << indent << "distance to nearest medial axis : "
<< medialDist_.name() << nl
- << "nearest medial axis location : "
+ << indent << "nearest medial axis location : "
<< medialVec_.name() << nl
- << "normal at nearest wall : "
+ << indent << "normal at nearest wall : "
<< dispVec_.name() << nl
- << decrIndent << nl
- << endl;
+ << decrIndent << endl;
dispVec_.write();
medialRatio_.write();
diff --git a/src/mesh/snappyHexMesh/meshRefinement/meshRefinementBaffles.C b/src/mesh/snappyHexMesh/meshRefinement/meshRefinementBaffles.C
index 440819950b2..eb0f8d5f350 100644
--- a/src/mesh/snappyHexMesh/meshRefinement/meshRefinementBaffles.C
+++ b/src/mesh/snappyHexMesh/meshRefinement/meshRefinementBaffles.C
@@ -4194,7 +4194,7 @@ Foam::autoPtr<Foam::mapPolyMesh> Foam::meshRefinement::removeLimitShells
if (debug&meshRefinement::MESH)
{
const_cast<Time&>(mesh_.time())++;
- Pout<< "Writing mesh after removing limitShells "
+ Pout<< "Writing mesh after removing limitShells"
<< " to time " << timeName() << endl;
write
(
--
GitLab