FreeFOAM The Cross-Platform CFD Toolkit
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) 1991-2010 OpenCFD Ltd.
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 <finiteVolume/fvMesh.H>
27 #include "fvMeshAdder.H"
29 #include <finiteVolume/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 // ************************ vim: set sw=4 sts=4 et: ************************ //