src/OpenFOAM/meshes/polyMesh/syncTools/syncTools.C

   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 "syncTools.H"
  27
  28 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  29
  30 // Determines for every point whether it is coupled and if so sets only one.
  31 Foam::PackedBoolList Foam::syncTools::getMasterPoints(const polyMesh& mesh)
  32 {
  33     PackedBoolList isMasterPoint(mesh.nPoints());
  34     PackedBoolList donePoint(mesh.nPoints());
  35
  36     const globalMeshData& globalData = mesh.globalData();
  37     const labelList& meshPoints = globalData.coupledPatch().meshPoints();
  38     const labelListList& slaves = globalData.globalPointSlaves();
  39     const labelListList& transformedSlaves =
  40             globalData.globalPointTransformedSlaves();
  41
  42     forAll(meshPoints, coupledPointI)
  43     {
  44         label meshPointI = meshPoints[coupledPointI];
  45         if
  46         (
  47             (
  48                 slaves[coupledPointI].size()
  49               + transformedSlaves[coupledPointI].size()
  50             )
  51           > 0
  52         )
  53         {
  54             isMasterPoint[meshPointI] = true;
  55         }
  56         donePoint[meshPointI] = true;
  57     }
  58
  59
  60     // Do all other points
  61     // ~~~~~~~~~~~~~~~~~~~
  62
  63     forAll(donePoint, pointI)
  64     {
  65         if (!donePoint[pointI])
  66         {
  67             isMasterPoint[pointI] = true;
  68         }
  69     }
  70
  71     return isMasterPoint;
  72 }
  73
  74
  75 // Determines for every edge whether it is coupled and if so sets only one.
  76 Foam::PackedBoolList Foam::syncTools::getMasterEdges(const polyMesh& mesh)
  77 {
  78     PackedBoolList isMasterEdge(mesh.nEdges());
  79     PackedBoolList doneEdge(mesh.nEdges());
  80
  81     const globalMeshData& globalData = mesh.globalData();
  82     const labelList& meshEdges = globalData.coupledPatchMeshEdges();
  83     const labelListList& slaves = globalData.globalEdgeSlaves();
  84     const labelListList& transformedSlaves =
  85         globalData.globalEdgeTransformedSlaves();
  86
  87     forAll(meshEdges, coupledEdgeI)
  88     {
  89         label meshEdgeI = meshEdges[coupledEdgeI];
  90         if
  91         (
  92             (
  93                 slaves[coupledEdgeI].size()
  94               + transformedSlaves[coupledEdgeI].size()
  95             )
  96           > 0
  97         )
  98         {
  99             isMasterEdge[meshEdgeI] = true;
 100         }
 101         doneEdge[meshEdgeI] = true;
 102     }
 103
 104
 105     // Do all other edges
 106     // ~~~~~~~~~~~~~~~~~~
 107
 108     forAll(doneEdge, edgeI)
 109     {
 110         if (!doneEdge[edgeI])
 111         {
 112             isMasterEdge[edgeI] = true;
 113         }
 114     }
 115
 116     return isMasterEdge;
 117 }
 118
 119
 120 // Determines for every face whether it is coupled and if so sets only one.
 121 Foam::PackedBoolList Foam::syncTools::getMasterFaces(const polyMesh& mesh)
 122 {
 123     PackedBoolList isMasterFace(mesh.nFaces(), 1);
 124
 125     const polyBoundaryMesh& patches = mesh.boundaryMesh();
 126
 127     forAll(patches, patchI)
 128     {
 129         if (patches[patchI].coupled())
 130         {
 131             const coupledPolyPatch& pp =
 132                 refCast<const coupledPolyPatch>(patches[patchI]);
 133
 134             if (!pp.owner())
 135             {
 136                 forAll(pp, i)
 137                 {
 138                     isMasterFace.unset(pp.start()+i);
 139                 }
 140             }
 141         }
 142     }
 143
 144     return isMasterFace;
 145 }
 146
 147
 148 // ************************************************************************* //