fvMeshAdder.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-2019 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 "fvMesh.H"
27 #include "fvMeshAdder.H"
28 #include "faceCoupleInfo.H"
29 #include "fvMesh.H"
30 
31 /* * * * * * * * * * * * * * * Static Member Data * * * * * * * * * * * * * */
32 
33 namespace Foam
34 {
35 defineTypeNameAndDebug(fvMeshAdder, 0);
36 }
37 
38 
39 // * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
40 
41 Foam::labelList Foam::fvMeshAdder::calcPatchMap
42 (
43  const label oldStart,
44  const label oldSize,
45  const labelList& oldToNew,
46  const polyPatch& newPatch,
47  const label unmappedValue
48 )
49 {
50  labelList newToOld(newPatch.size(), unmappedValue);
51 
52  label newStart = newPatch.start();
53  label newSize = newPatch.size();
54 
55  for (label i = 0; i < oldSize; i++)
56  {
57  label newFacei = oldToNew[oldStart+i];
58 
59  if (newFacei >= newStart && newFacei < newStart+newSize)
60  {
61  newToOld[newFacei-newStart] = i;
62  }
63  }
64  return newToOld;
65 }
66 
67 
68 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
69 
71 (
72  fvMesh& mesh0,
73  const fvMesh& mesh1,
74  const faceCoupleInfo& coupleInfo,
75  const bool validBoundary
76 )
77 {
78  // Store old mesh0 point maps
79  labelListList oldMeshPoints0;
80  const bool havePointMesh =
81  mesh0.foundObject<pointMesh>(pointMesh::typeName);
82  if (havePointMesh)
83  {
84  const polyBoundaryMesh& pbm0 = mesh0.boundaryMesh();
85  oldMeshPoints0.setSize(pbm0.size());
86  forAll(pbm0, patchi)
87  {
88  oldMeshPoints0[patchi] = pbm0[patchi].meshPoints();
89  }
90  }
91 
92  // Resulting merged mesh (polyMesh only!)
94  (
96  (
97  mesh0,
98  mesh1,
99  coupleInfo,
100  validBoundary
101  )
102  );
103 
104  // Adjust the fvMesh part.
105  const polyBoundaryMesh& patches = mesh0.boundaryMesh();
106 
107  fvBoundaryMesh& fvPatches = const_cast<fvBoundaryMesh&>(mesh0.boundary());
108  fvPatches.setSize(patches.size());
109  forAll(patches, patchi)
110  {
111  fvPatches.set(patchi, fvPatch::New(patches[patchi], fvPatches));
112  }
113 
114  // Do the mapping of the stored fields
115  // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
116  MapVolFields<scalar>(mapPtr, mesh0, mesh1);
117  MapVolFields<vector>(mapPtr, mesh0, mesh1);
118  MapVolFields<sphericalTensor>(mapPtr, mesh0, mesh1);
119  MapVolFields<symmTensor>(mapPtr, mesh0, mesh1);
120  MapVolFields<tensor>(mapPtr, mesh0, mesh1);
121 
122  MapSurfaceFields<scalar>(mapPtr, mesh0, mesh1);
123  MapSurfaceFields<vector>(mapPtr, mesh0, mesh1);
124  MapSurfaceFields<sphericalTensor>(mapPtr, mesh0, mesh1);
125  MapSurfaceFields<symmTensor>(mapPtr, mesh0, mesh1);
126  MapSurfaceFields<tensor>(mapPtr, mesh0, mesh1);
127 
128  if (havePointMesh)
129  {
130  // Recreate point mesh
131  const pointMesh& pointMesh0 = pointMesh::New(mesh0);
132 
133  MapPointFields<scalar>(mapPtr, pointMesh0, oldMeshPoints0, mesh1);
134  MapPointFields<vector>(mapPtr, pointMesh0, oldMeshPoints0, mesh1);
135  MapPointFields<sphericalTensor>
136  (
137  mapPtr,
138  pointMesh0,
139  oldMeshPoints0,
140  mesh1
141  );
142  MapPointFields<symmTensor>(mapPtr, pointMesh0, oldMeshPoints0, mesh1);
143  MapPointFields<tensor>(mapPtr, pointMesh0, oldMeshPoints0, mesh1);
144  }
145 
146  MapDimFields<scalar>(mapPtr, mesh0, mesh1);
147  MapDimFields<vector>(mapPtr, mesh0, mesh1);
148  MapDimFields<sphericalTensor>(mapPtr, mesh0, mesh1);
149  MapDimFields<symmTensor>(mapPtr, mesh0, mesh1);
150  MapDimFields<tensor>(mapPtr, mesh0, mesh1);
151 
152  return mapPtr;
153 }
154 
155 
156 // ************************************************************************* //
const polyBoundaryMesh & boundaryMesh() const
Return boundary mesh.
Definition: polyMesh.H:434
#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
bool set(const label) const
Is element set.
Definition: PtrListI.H:65
bool foundObject(const word &name) const
Is the named Type found?
patches[0]
Container for information needed to couple to meshes. When constructed from two meshes and a list of ...
Mesh representing a set of points created from polyMesh.
Definition: pointMesh.H:48
List< label > labelList
A List of labels.
Definition: labelList.H:56
static autoPtr< mapAddedPolyMesh > add(fvMesh &mesh0, const fvMesh &mesh1, const faceCoupleInfo &coupleInfo, const bool validBoundary=true)
Inplace add mesh to fvMesh. Maps all stored fields. Returns map.
Definition: fvMeshAdder.C:71
void setSize(const label)
Reset size of PtrList. If extending the PtrList, new entries are.
Definition: PtrList.C:131
Foam::polyBoundaryMesh.
defineTypeNameAndDebug(combustionModel, 0)
label size() const
Return the number of elements in the UPtrList.
Definition: UPtrListI.H:29
void setSize(const label)
Reset size of List.
Definition: List.C:281
label patchi
Foam::fvBoundaryMesh.
Mesh data needed to do the Finite Volume discretisation.
Definition: fvMesh.H:78
An auto-pointer similar to the STL auto_ptr but with automatic casting to a reference to the type and...
Definition: PtrList.H:52
static autoPtr< fvPatch > New(const polyPatch &, const fvBoundaryMesh &)
Return a pointer to a new patch created on freestore from polyPatch.
Definition: fvPatchNew.C:32
static autoPtr< polyMesh > add(const IOobject &io, const polyMesh &mesh0, const polyMesh &mesh1, const faceCoupleInfo &coupleInfo, autoPtr< mapAddedPolyMesh > &mapPtr)
Add two polyMeshes. Returns new polyMesh and map construct.
Namespace for OpenFOAM.
const fvBoundaryMesh & boundary() const
Return reference to boundary mesh.
Definition: fvMesh.C:540