songtianlun/openfoam
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

ENH: use vtk::surfaceWriter in shortestPathSet debugging

Browse Source
This commit is contained in:
Mark Olesen 2019-02-07 09:45:13 +01:00 committed by Andrew Heather
parent 3a6e596a09
commit 2016f88eba
1 changed files with 29 additions and 35 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

64
src/sampling/sampledSet/shortestPath/shortestPathSet.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2017-2018 OpenCFD Ltd.
\\ / A nd | Copyright (C) 2017-2019 OpenCFD Ltd.
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -33,10 +33,8 @@ License
#include "fvMesh.H"
#include "volFields.H"
#include "OBJstream.H"
#include "surfaceWriter.H"
#include "meshedSurfRef.H"
#include "PatchTools.H"
#include "vtkSurfaceWriter.H"
#include "foamVtkSurfaceWriter.H"
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
@ -583,9 +581,9 @@ void Foam::shortestPathSet::genSamples
}
}
if
if
(
frontCellI != -1
frontCellI != -1
&& allCellInfo[frontCellI].distance() < minCellDistance
)
{
@ -840,20 +838,16 @@ Foam::shortestPathSet::shortestPathSet
{
if (debug)
{
fileName outputDir;
if (Pstream::master())
{
outputDir =
mesh.time().globalPath()
/ functionObject::outputPrefix
/ mesh.pointsInstance();
outputDir.clean();
mkDir(outputDir);
Info<< "shortestPathSet : Writing blocked faces to "
<< outputDir << endl;
}
fileName outputDir =
(
mesh.time().globalPath()
/ functionObject::outputPrefix
/ mesh.pointsInstance()
);
outputDir.clean();
const vtkSurfaceWriter surfWriter;
Info<< "shortestPathSet : Writing blocked faces to "
<< outputDir << endl;
const indirectPrimitivePatch setPatch
(
@ -864,8 +858,10 @@ Foam::shortestPathSet::shortestPathSet
),
mesh.points()
);
if (Pstream::parRun())
{
// Topological merge
labelList pointToGlobal;
labelList uniqueMeshPointLabels;
autoPtr<globalIndex> globalPoints;
@ -891,30 +887,28 @@ Foam::shortestPathSet::shortestPathSet
// Write
if (Pstream::master())
{
surfWriter.write
vtk::surfaceWriter writer
(
outputDir,
"blockedFace",
meshedSurfRef
(
mergedPoints,
mergedFaces
)
mergedPoints,
mergedFaces,
(outputDir / "blockedFace"),
false // serial only - already merged
);
writer.writeGeometry();
}
}
else
{
surfWriter.write
vtk::surfaceWriter writer
(
outputDir,
"blockedFace",
meshedSurfRef
(
setPatch.localPoints(),
setPatch.localFaces()
)
setPatch.localPoints(),
setPatch.localFaces(),
(outputDir / "blockedFace"),
false // serial only - redundant
);
writer.writeGeometry();
}
}