All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
patchWriter.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration | Website: https://openfoam.org
5  \\ / A nd | Copyright (C) 2011-2020 OpenFOAM Foundation
6  \\/ M anipulation |
7 -------------------------------------------------------------------------------
8 License
9  This file is part of OpenFOAM.
10 
11  OpenFOAM is free software: you can redistribute it and/or modify it
12  under the terms of the GNU General Public License as published by
13  the Free Software Foundation, either version 3 of the License, or
14  (at your option) any later version.
15 
16  OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
17  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
18  FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
19  for more details.
20 
21  You should have received a copy of the GNU General Public License
22  along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
23 
24 \*---------------------------------------------------------------------------*/
25 
26 #include "patchWriter.H"
27 #include "vtkWriteFieldOps.H"
28 
29 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
30 
32 (
33  const vtkMesh& vMesh,
34  const bool binary,
35  const bool nearCellValue,
36  const fileName& fName,
37  const labelList& patchIDs
38 )
39 :
40  vMesh_(vMesh),
41  binary_(binary),
42  nearCellValue_(nearCellValue),
43  fName_(fName),
44  patchIDs_(patchIDs),
45  os_(fName.c_str())
46 {
47  const fvMesh& mesh = vMesh_.mesh();
48  const polyBoundaryMesh& patches = mesh.boundaryMesh();
49 
50  // Write header
51  if (patchIDs_.size() == 1)
52  {
53  vtkWriteOps::writeHeader(os_, binary_, patches[patchIDs_[0]].name());
54  }
55  else
56  {
57  vtkWriteOps::writeHeader(os_, binary_, "patches");
58  }
59  os_ << "DATASET POLYDATA" << std::endl;
60 
61  // Write topology
62  nPoints_ = 0;
63  nFaces_ = 0;
64  label nFaceVerts = 0;
65 
66  forAll(patchIDs_, i)
67  {
68  const polyPatch& pp = patches[patchIDs_[i]];
69 
70  nPoints_ += pp.nPoints();
71  nFaces_ += pp.size();
72 
73  forAll(pp, facei)
74  {
75  nFaceVerts += pp[facei].size() + 1;
76  }
77  }
78 
79  os_ << "POINTS " << nPoints_ << " float" << std::endl;
80 
81  DynamicList<floatScalar> ptField(3*nPoints_);
82 
83  forAll(patchIDs_, i)
84  {
85  const polyPatch& pp = patches[patchIDs_[i]];
86 
87  vtkWriteOps::insert(pp.localPoints(), ptField);
88  }
89  vtkWriteOps::write(os_, binary_, ptField);
90 
91  os_ << "POLYGONS " << nFaces_ << ' ' << nFaceVerts << std::endl;
92 
93  DynamicList<label> vertLabels(nFaceVerts);
94 
95  label offset = 0;
96 
97  forAll(patchIDs_, i)
98  {
99  const polyPatch& pp = patches[patchIDs_[i]];
100 
101  forAll(pp, facei)
102  {
103  const face& f = pp.localFaces()[facei];
104 
105  vertLabels.append(f.size());
106  vtkWriteOps::insert(f + offset, vertLabels);
107  }
108  offset += pp.nPoints();
109  }
110  vtkWriteOps::write(os_, binary_, vertLabels);
111 }
112 
113 
114 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
115 
117 {
118  const fvMesh& mesh = vMesh_.mesh();
119 
120  DynamicList<floatScalar> fField(nFaces_);
121 
122  os_ << "patchID 1 " << nFaces_ << " float" << std::endl;
123 
124  forAll(patchIDs_, i)
125  {
126  label patchi = patchIDs_[i];
127 
128  const polyPatch& pp = mesh.boundaryMesh()[patchi];
129 
130  if (!isA<emptyPolyPatch>(pp))
131  {
132  vtkWriteOps::insert(scalarField(pp.size(), patchi), fField);
133  }
134  }
135  vtkWriteOps::write(os_, binary_, fField);
136 }
137 
138 
139 // ************************************************************************* //
const polyBoundaryMesh & boundaryMesh() const
Return boundary mesh.
Definition: polyMesh.H:434
const fvMesh & mesh() const
Access either mesh or submesh.
Definition: vtkMesh.H:100
#define forAll(list, i)
Loop across all elements in list.
Definition: UList.H:434
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
Definition: label.H:59
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251
void writePatchIDs()
Write cellIDs.
void insert(const scalar, DynamicList< floatScalar > &)
Append scalar to given DynamicList.
void writeHeader(std::ostream &, const bool isBinary, const std::string &title)
Write header.
Field< scalar > scalarField
Specialisation of Field<T> for scalar.
void write(std::ostream &os, const bool binary, List< floatScalar > &fField)
Write floats ascii or binary.
List< label > labelList
A List of labels.
Definition: labelList.H:56
word name(const complex &)
Return a string representation of a complex.
Definition: complex.C:47
label patchi
patchWriter(const vtkMesh &, const bool binary, const bool nearCellValue, const fileName &, const labelList &patchIDs)
Construct from components.