diff --git a/applications/utilities/mesh/generation/foamyMesh/foamyHexMesh/foamyHexMesh.C b/applications/utilities/mesh/generation/foamyMesh/foamyHexMesh/foamyHexMesh.C index caff0c5474..bbb1c2630d 100644 --- a/applications/utilities/mesh/generation/foamyMesh/foamyHexMesh/foamyHexMesh.C +++ b/applications/utilities/mesh/generation/foamyMesh/foamyHexMesh/foamyHexMesh.C @@ -6,7 +6,7 @@ \\/ M anipulation | ------------------------------------------------------------------------------- Copyright (C) 2011-2016 OpenFOAM Foundation - Copyright (C) 2016-2020 OpenCFD Ltd. + Copyright (C) 2016-2022 OpenCFD Ltd. ------------------------------------------------------------------------------- License This file is part of OpenFOAM. @@ -40,7 +40,7 @@ Description #include "IOdictionary.H" #include "searchableSurfaces.H" #include "conformalVoronoiMesh.H" -#include "vtkSetWriter.H" +#include "vtkCoordSetWriter.H" using namespace Foam; @@ -106,8 +106,10 @@ int main(int argc, char *argv[]) foamyHexMeshDict.getOrDefault("singleRegionName", true) ); + autoPtr setWriterPtr(new coordSetWriters::vtkWriter()); + // Write some stats - allGeometry.writeStats(List(0), Info); + allGeometry.writeStats(List(), Info); // Check topology allGeometry.checkTopology(true); // Check geometry @@ -115,7 +117,7 @@ int main(int argc, char *argv[]) ( 100.0, // max size ratio 1e-9, // intersection tolerance - autoPtr>(new vtkSetWriter()), + setWriterPtr, 0.01, // min triangle quality true ); diff --git a/applications/utilities/mesh/generation/snappyHexMesh/snappyHexMesh.C b/applications/utilities/mesh/generation/snappyHexMesh/snappyHexMesh.C index 1d18e68fb9..12959d2361 100644 --- a/applications/utilities/mesh/generation/snappyHexMesh/snappyHexMesh.C +++ b/applications/utilities/mesh/generation/snappyHexMesh/snappyHexMesh.C @@ -6,7 +6,7 @@ \\/ M anipulation | ------------------------------------------------------------------------------- Copyright (C) 2011-2016 OpenFOAM Foundation - Copyright (C) 2015-2021 OpenCFD Ltd. + Copyright (C) 2015-2022 OpenCFD Ltd. ------------------------------------------------------------------------------- License This file is part of OpenFOAM. @@ -51,7 +51,7 @@ Description #include "refinementParameters.H" #include "snapParameters.H" #include "layerParameters.H" -#include "vtkSetWriter.H" +#include "vtkCoordSetWriter.H" #include "faceSet.H" #include "motionSmoother.H" #include "polyTopoChange.H" @@ -854,19 +854,25 @@ int main(int argc, char *argv[]) const bool keepPatches(meshDict.getOrDefault("keepPatches", false)); - // format to be used for writing lines - const word setFormat - ( - meshDict.getOrDefault + // Writer for writing lines + autoPtr setFormatter; + { + const word setFormat ( - "setFormat", - vtkSetWriter::typeName - ) - ); - const autoPtr> setFormatter - ( - writer::New(setFormat) - ); + meshDict.getOrDefault + ( + "setFormat", + coordSetWriters::vtkWriter::typeName // Default: "vtk" + ) + ); + + setFormatter = coordSetWriter::New + ( + setFormat, + meshDict.subOrEmptyDict("formatOptions").optionalSubDict(setFormat) + ); + } + const scalar maxSizeRatio ( diff --git a/applications/utilities/mesh/manipulation/checkMesh/checkGeometry.C b/applications/utilities/mesh/manipulation/checkMesh/checkGeometry.C index b0c4a33676..ce2c8f83ce 100644 --- a/applications/utilities/mesh/manipulation/checkMesh/checkGeometry.C +++ b/applications/utilities/mesh/manipulation/checkMesh/checkGeometry.C @@ -11,7 +11,7 @@ #include "checkTools.H" #include "functionObject.H" -#include "vtkSetWriter.H" +#include "vtkCoordSetWriter.H" #include "vtkSurfaceWriter.H" #include "cyclicACMIPolyPatch.H" @@ -537,7 +537,7 @@ Foam::label Foam::checkGeometry const polyMesh& mesh, const bool allGeometry, autoPtr& surfWriter, - const autoPtr>& setWriter + autoPtr& setWriter ) { label noFailedChecks = 0; @@ -594,7 +594,7 @@ Foam::label Foam::checkGeometry << nonAlignedPoints.name() << endl; nonAlignedPoints.instance() = mesh.pointsInstance(); nonAlignedPoints.write(); - if (setWriter) + if (setWriter && setWriter->enabled()) { mergeAndWrite(*setWriter, nonAlignedPoints); } @@ -628,7 +628,7 @@ Foam::label Foam::checkGeometry << " non closed cells to set " << cells.name() << endl; cells.instance() = mesh.pointsInstance(); cells.write(); - if (surfWriter) + if (surfWriter && surfWriter->enabled()) { mergeAndWrite(*surfWriter, cells); } @@ -644,7 +644,7 @@ Foam::label Foam::checkGeometry << aspectCells.name() << endl; aspectCells.instance() = mesh.pointsInstance(); aspectCells.write(); - if (surfWriter) + if (surfWriter && surfWriter->enabled()) { mergeAndWrite(*surfWriter, aspectCells); } @@ -665,7 +665,7 @@ Foam::label Foam::checkGeometry << " zero area faces to set " << faces.name() << endl; faces.instance() = mesh.pointsInstance(); faces.write(); - if (surfWriter) + if (surfWriter && surfWriter->enabled()) { mergeAndWrite(*surfWriter, faces); } @@ -687,7 +687,7 @@ Foam::label Foam::checkGeometry << " zero volume cells to set " << cells.name() << endl; cells.instance() = mesh.pointsInstance(); cells.write(); - if (surfWriter) + if (surfWriter && surfWriter->enabled()) { mergeAndWrite(*surfWriter, cells); } @@ -710,7 +710,7 @@ Foam::label Foam::checkGeometry << " non-orthogonal faces to set " << faces.name() << endl; faces.instance() = mesh.pointsInstance(); faces.write(); - if (surfWriter) + if (surfWriter && surfWriter->enabled()) { mergeAndWrite(*surfWriter, faces); } @@ -732,7 +732,7 @@ Foam::label Foam::checkGeometry << faces.name() << endl; faces.instance() = mesh.pointsInstance(); faces.write(); - if (surfWriter) + if (surfWriter && surfWriter->enabled()) { mergeAndWrite(*surfWriter, faces); } @@ -754,7 +754,7 @@ Foam::label Foam::checkGeometry << " skew faces to set " << faces.name() << endl; faces.instance() = mesh.pointsInstance(); faces.write(); - if (surfWriter) + if (surfWriter && surfWriter->enabled()) { mergeAndWrite(*surfWriter, faces); } @@ -778,7 +778,7 @@ Foam::label Foam::checkGeometry << faces.name() << endl; faces.instance() = mesh.pointsInstance(); faces.write(); - if (surfWriter) + if (surfWriter && surfWriter->enabled()) { mergeAndWrite(*surfWriter, faces); } @@ -811,7 +811,7 @@ Foam::label Foam::checkGeometry << "decomposition tets to set " << faces.name() << endl; faces.instance() = mesh.pointsInstance(); faces.write(); - if (surfWriter) + if (surfWriter && surfWriter->enabled()) { mergeAndWrite(*surfWriter, faces); } @@ -836,7 +836,7 @@ Foam::label Foam::checkGeometry << endl; points.instance() = mesh.pointsInstance(); points.write(); - if (setWriter) + if (setWriter && setWriter->enabled()) { mergeAndWrite(*setWriter, points); } @@ -859,7 +859,7 @@ Foam::label Foam::checkGeometry << " apart) points to set " << nearPoints.name() << endl; nearPoints.instance() = mesh.pointsInstance(); nearPoints.write(); - if (setWriter) + if (setWriter && setWriter->enabled()) { mergeAndWrite(*setWriter, nearPoints); } @@ -883,7 +883,7 @@ Foam::label Foam::checkGeometry << endl; faces.instance() = mesh.pointsInstance(); faces.write(); - if (surfWriter) + if (surfWriter && surfWriter->enabled()) { mergeAndWrite(*surfWriter, faces); } @@ -906,7 +906,7 @@ Foam::label Foam::checkGeometry << " warped faces to set " << faces.name() << endl; faces.instance() = mesh.pointsInstance(); faces.write(); - if (surfWriter) + if (surfWriter && surfWriter->enabled()) { mergeAndWrite(*surfWriter, faces); } @@ -927,7 +927,7 @@ Foam::label Foam::checkGeometry << " under-determined cells to set " << cells.name() << endl; cells.instance() = mesh.pointsInstance(); cells.write(); - if (surfWriter) + if (surfWriter && surfWriter->enabled()) { mergeAndWrite(*surfWriter, cells); } @@ -947,7 +947,7 @@ Foam::label Foam::checkGeometry << " concave cells to set " << cells.name() << endl; cells.instance() = mesh.pointsInstance(); cells.write(); - if (surfWriter) + if (surfWriter && surfWriter->enabled()) { mergeAndWrite(*surfWriter, cells); } @@ -968,7 +968,7 @@ Foam::label Foam::checkGeometry << faces.name() << endl; faces.instance() = mesh.pointsInstance(); faces.write(); - if (surfWriter) + if (surfWriter && surfWriter->enabled()) { mergeAndWrite(*surfWriter, faces); } @@ -989,7 +989,7 @@ Foam::label Foam::checkGeometry << faces.name() << endl; faces.instance() = mesh.pointsInstance(); faces.write(); - if (surfWriter) + if (surfWriter && surfWriter->enabled()) { mergeAndWrite(*surfWriter, faces); } diff --git a/applications/utilities/mesh/manipulation/checkMesh/checkGeometry.H b/applications/utilities/mesh/manipulation/checkMesh/checkGeometry.H index e69ab01e0f..03d6da2f54 100644 --- a/applications/utilities/mesh/manipulation/checkMesh/checkGeometry.H +++ b/applications/utilities/mesh/manipulation/checkMesh/checkGeometry.H @@ -7,7 +7,7 @@ namespace Foam // Forward Declarations class polyMesh; class wedgePolyPatch; - template class writer; + class coordSetWriter; class surfaceWriter; label findOppositeWedge(const polyMesh&, const wedgePolyPatch&); @@ -47,6 +47,6 @@ namespace Foam const polyMesh& mesh, const bool allGeometry, autoPtr& surfWriter, - const autoPtr>& setWriter + autoPtr& setWriter ); } diff --git a/applications/utilities/mesh/manipulation/checkMesh/checkMesh.C b/applications/utilities/mesh/manipulation/checkMesh/checkMesh.C index 2455b93e4f..76db56ca70 100644 --- a/applications/utilities/mesh/manipulation/checkMesh/checkMesh.C +++ b/applications/utilities/mesh/manipulation/checkMesh/checkMesh.C @@ -74,7 +74,7 @@ Usage #include "Time.H" #include "fvMesh.H" #include "globalMeshData.H" -#include "vtkSetWriter.H" +#include "vtkCoordSetWriter.H" #include "vtkSurfaceWriter.H" #include "IOdictionary.H" #include "regionProperties.H" @@ -262,11 +262,11 @@ int main(int argc, char *argv[]) autoPtr surfWriter; - autoPtr> setWriter; + autoPtr setWriter; if (writeSets) { surfWriter = surfaceWriter::New(surfaceFormat); - setWriter = writer::New(vtkSetWriter::typeName); + setWriter = coordSetWriter::New(coordSetWriters::vtkWriter::typeName); } diff --git a/applications/utilities/mesh/manipulation/checkMesh/checkTools.C b/applications/utilities/mesh/manipulation/checkMesh/checkTools.C index 00c1f4b0f2..7e461b37a9 100644 --- a/applications/utilities/mesh/manipulation/checkMesh/checkTools.C +++ b/applications/utilities/mesh/manipulation/checkMesh/checkTools.C @@ -41,7 +41,7 @@ License #include "faceSet.H" #include "cellSet.H" #include "Time.H" -#include "writer.H" +#include "coordSetWriter.H" #include "surfaceWriter.H" #include "syncTools.H" #include "globalIndex.H" @@ -417,7 +417,7 @@ void Foam::mergeAndWrite void Foam::mergeAndWrite ( - const writer& writer, + coordSetWriter& writer, const pointSet& set ) { @@ -441,39 +441,26 @@ void Foam::mergeAndWrite } - // Write with scalar pointID + // Write with pointID if (Pstream::master()) { - scalarField scalarPointIDs(mergedIDs.size()); - forAll(mergedIDs, i) - { - scalarPointIDs[i] = 1.0*mergedIDs[i]; - } + coordSet coords(set.name(), "distance", mergedPts, mag(mergedPts)); - coordSet points(set.name(), "distance", mergedPts, mag(mergedPts)); + // Output. E.g. pointSet p0 -> postProcessing/