src/dynamicMesh/fvMeshAdder/fvMeshAdder.C

   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"
  28 #include <dynamicMesh/faceCoupleInfo.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
  65 Foam::autoPtr<Foam::mapAddedPolyMesh> Foam::fvMeshAdder::add
  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!)
  76     autoPtr<mapAddedPolyMesh> mapPtr
  77     (
  78         polyMeshAdder::add
  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: ************************ //