fvMeshAdder.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration |
5  \\ / A nd | Copyright (C) 2011 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 // * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
32 
33 //- Calculate map from new patch faces to old patch faces. -1 where
34 // could not map.
35 Foam::labelList Foam::fvMeshAdder::calcPatchMap
36 (
37  const label oldStart,
38  const label oldSize,
39  const labelList& oldToNew,
40  const polyPatch& newPatch,
41  const label unmappedValue
42 )
43 {
44  labelList newToOld(newPatch.size(), unmappedValue);
45 
46  label newStart = newPatch.start();
47  label newSize = newPatch.size();
48 
49  for (label i = 0; i < oldSize; i++)
50  {
51  label newFaceI = oldToNew[oldStart+i];
52 
53  if (newFaceI >= newStart && newFaceI < newStart+newSize)
54  {
55  newToOld[newFaceI-newStart] = i;
56  }
57  }
58  return newToOld;
59 }
60 
61 
62 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
63 
64 // Inplace add mesh1 to mesh0
66 (
67  fvMesh& mesh0,
68  const fvMesh& mesh1,
69  const faceCoupleInfo& coupleInfo,
70  const bool validBoundary
71 )
72 {
73  mesh0.clearOut();
74 
75  // Resulting merged mesh (polyMesh only!)
77  (
79  (
80  mesh0,
81  mesh1,
82  coupleInfo,
83  validBoundary
84  )
85  );
86 
87  // Adjust the fvMesh part.
88  const polyBoundaryMesh& patches = mesh0.boundaryMesh();
89 
90  fvBoundaryMesh& fvPatches = const_cast<fvBoundaryMesh&>(mesh0.boundary());
91  fvPatches.setSize(patches.size());
92  forAll(patches, patchI)
93  {
94  fvPatches.set(patchI, fvPatch::New(patches[patchI], fvPatches));
95  }
96 
97  // Do the mapping of the stored fields
98  // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
99  fvMeshAdder::MapVolFields<scalar>(mapPtr, mesh0, mesh1);
100  fvMeshAdder::MapVolFields<vector>(mapPtr, mesh0, mesh1);
101  fvMeshAdder::MapVolFields<sphericalTensor>(mapPtr, mesh0, mesh1);
102  fvMeshAdder::MapVolFields<symmTensor>(mapPtr, mesh0, mesh1);
103  fvMeshAdder::MapVolFields<tensor>(mapPtr, mesh0, mesh1);
104 
105  fvMeshAdder::MapSurfaceFields<scalar>(mapPtr, mesh0, mesh1);
106  fvMeshAdder::MapSurfaceFields<vector>(mapPtr, mesh0, mesh1);
107  fvMeshAdder::MapSurfaceFields<sphericalTensor>(mapPtr, mesh0, mesh1);
108  fvMeshAdder::MapSurfaceFields<symmTensor>(mapPtr, mesh0, mesh1);
109  fvMeshAdder::MapSurfaceFields<tensor>(mapPtr, mesh0, mesh1);
110 
111  return mapPtr;
112 }
113 
114 
115 // ************************************************************************* //
void setSize(const label)
Reset size of PtrList. If extending the PtrList, new entries are.
Definition: PtrList.C:142
static autoPtr< fvPatch > New(const polyPatch &, const fvBoundaryMesh &)
Return a pointer to a new patch created on freestore from polyPatch.
Definition: fvPatchNew.C:33
Mesh data needed to do the Finite Volume discretisation.
Definition: fvMesh.H:78
label size() const
Return the number of elements in the PtrList.
Definition: PtrListI.H:32
bool set(const label) const
Is element set.
Definition: PtrListI.H:107
Foam::fvBoundaryMesh.
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
Container for information needed to couple to meshes. When constructed from two meshes and a geometri...
patches[0]
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.
Foam::polyBoundaryMesh.
#define forAll(list, i)
Definition: UList.H:421
const polyBoundaryMesh & boundaryMesh() const
Return boundary mesh.
Definition: polyMesh.H:421
List< label > labelList
A List of labels.
Definition: labelList.H:56
void clearOut()
Clear all geometry and addressing.
Definition: fvMesh.C:231
An auto-pointer similar to the STL auto_ptr but with automatic casting to a reference to the type and...
Definition: PtrList.H:117
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:66
const fvBoundaryMesh & boundary() const
Return reference to boundary mesh.
Definition: fvMesh.C:552