src/OpenFOAM/meshes/primitiveMesh/primitivePatch/patchZones.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 "patchZones.H"
  27
  28
  29 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  30
  31 defineTypeNameAndDebug(Foam::patchZones, 0);
  32
  33
  34 // * * * * * * * * * * * * * * * Global Functions  * * * * * * * * * * * * * //
  35
  36 // Gets labels of changed faces and propagates them to the edges. Returns
  37 // labels of edges changed.
  38 Foam::labelList Foam::patchZones::faceToEdge
  39 (
  40     const labelList& changedFaces,
  41     labelList& edgeRegion
  42 )
  43 {
  44     labelList changedEdges(pp_.nEdges(), -1);
  45     label changedI = 0;
  46
  47     forAll(changedFaces, i)
  48     {
  49         label faceI = changedFaces[i];
  50
  51         const labelList& fEdges = pp_.faceEdges()[faceI];
  52
  53         forAll(fEdges, fEdgeI)
  54         {
  55             label edgeI = fEdges[fEdgeI];
  56
  57             if (!borderEdge_[edgeI] && (edgeRegion[edgeI] == -1))
  58             {
  59                 edgeRegion[edgeI] = nZones_;
  60
  61                 changedEdges[changedI++] = edgeI;
  62             }
  63         }
  64     }
  65
  66     changedEdges.setSize(changedI);
  67
  68     return changedEdges;
  69 }
  70
  71
  72 // Reverse of faceToEdge: gets edges and returns faces
  73 Foam::labelList Foam::patchZones::edgeToFace(const labelList& changedEdges)
  74 {
  75     labelList changedFaces(pp_.size(), -1);
  76     label changedI = 0;
  77
  78     forAll(changedEdges, i)
  79     {
  80         label edgeI = changedEdges[i];
  81
  82         const labelList& eFaces = pp_.edgeFaces()[edgeI];
  83
  84         forAll(eFaces, eFaceI)
  85         {
  86             label faceI = eFaces[eFaceI];
  87
  88             if (operator[](faceI) == -1)
  89             {
  90                 operator[](faceI) = nZones_;
  91
  92                 changedFaces[changedI++] = faceI;
  93             }
  94         }
  95     }
  96
  97     changedFaces.setSize(changedI);
  98
  99     return changedFaces;
 100 }
 101
 102
 103 // Finds area, starting at faceI, delimited by borderEdge
 104 void Foam::patchZones::markZone(label faceI)
 105 {
 106     // List of faces whose faceZone has been set.
 107     labelList changedFaces(1, faceI);
 108     // List of edges whose faceZone has been set.
 109     labelList changedEdges;
 110
 111     // Zones on all edges.
 112     labelList edgeZone(pp_.nEdges(), -1);
 113
 114     while (true)
 115     {
 116         changedEdges = faceToEdge(changedFaces, edgeZone);
 117
 118         if (debug)
 119         {
 120             Info<< "From changedFaces:" << changedFaces.size()
 121                 << " to changedEdges:" << changedEdges.size()
 122                 << endl;
 123         }
 124
 125         if (changedEdges.empty())
 126         {
 127             break;
 128         }
 129
 130         changedFaces = edgeToFace(changedEdges);
 131
 132         if (debug)
 133         {
 134             Info<< "From changedEdges:" << changedEdges.size()
 135                 << " to changedFaces:" << changedFaces.size()
 136                 << endl;
 137         }
 138
 139         if (changedEdges.empty())
 140         {
 141             break;
 142         }
 143     }
 144 }
 145
 146
 147 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 148
 149 // Construct from components
 150 Foam::patchZones::patchZones
 151 (
 152     const primitivePatch& pp,
 153     const boolList& borderEdge
 154 )
 155 :
 156     labelList(pp.size(), -1),
 157     pp_(pp),
 158     borderEdge_(borderEdge),
 159     nZones_(0)
 160 {
 161     // Finds areas delimited by borderEdge (or 'real' edges).
 162     // Fills *this with zone number accordingly.
 163
 164     if (borderEdge.size() != pp_.nEdges())
 165     {
 166         FatalErrorIn
 167         (
 168             "patchZones::patchZones(const primitivePatch&, const boolList&)"
 169         )   << "borderEdge boolList not same size as number of edges" << endl
 170             << "borderEdge:" << borderEdge.size() << endl
 171             << "nEdges    :" << pp_.nEdges()
 172             << abort(FatalError);
 173     }
 174
 175     label faceI = 0;
 176
 177     while (true)
 178     {
 179         // Find first non-visited face
 180         for (; faceI < pp_.size(); faceI++)
 181         {
 182             if (operator[](faceI) == -1)
 183             {
 184                 operator[](faceI) = nZones_;
 185
 186                 markZone(faceI);
 187
 188                 break;
 189             }
 190         }
 191
 192         if (faceI == pp_.size())
 193         {
 194             // Finished.
 195             break;
 196         }
 197
 198         nZones_++;
 199     }
 200 }
 201
 202
 203 // ************************************************************************* //