src/mesh/cfMesh/utilities/meshes/polyMeshGenModifier/polyMeshGenModifierAddCellByCell.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | cfMesh: A library for mesh generation
   4    \\    /   O peration     |
   5     \\  /    A nd           | Author: Franjo Juretic (franjo.juretic@c-fields.com)
   6      \\/     M anipulation  | Copyright (C) Creative Fields, Ltd.
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of cfMesh.
  10
  11     cfMesh is free software; you can redistribute it and/or modify it
  12     under the terms of the GNU General Public License as published by the
  13     Free Software Foundation; either version 3 of the License, or (at your
  14     option) any later version.
  15
  16     cfMesh 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 cfMesh.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 Description
  25
  26 \*---------------------------------------------------------------------------*/
  27
  28 #include "polyMeshGenModifierAddCellByCell.H"
  29 #include "demandDrivenData.H"
  30
  31 namespace Foam
  32 {
  33
  34 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  35
  36 polyMeshGenModifierAddCellByCell::polyMeshGenModifierAddCellByCell
  37 (
  38     polyMeshGen& mesh
  39 )
  40 :
  41     polyMeshGenModifier(mesh),
  42     nFaces_(mesh.faces().size()),
  43     newFaces_(nFaces_),
  44     nCells_(mesh.cells().size()),
  45     newCells_(nCells_)
  46 {
  47     this->pointFaces();
  48     faceListPMG& faces = this->facesAccess();
  49     forAll(faces, faceI)
  50         newFaces_[faceI].transfer(faces[faceI]);
  51
  52     cellListPMG& cells = this->cellsAccess();
  53     forAll(cells, cellI)
  54         newCells_[cellI].transfer(cells[cellI]);
  55 };
  56
  57 // Destructor
  58 polyMeshGenModifierAddCellByCell::~polyMeshGenModifierAddCellByCell()
  59 {
  60     faceListPMG& faces = this->facesAccess();
  61     faces.setSize(nFaces_);
  62     forAll(faces, faceI)
  63         faces[faceI].transfer(newFaces_[faceI]);
  64
  65     cellListPMG& cells = this->cellsAccess();
  66     cells.setSize(newCells_.size());
  67     forAll(cells, cellI)
  68         cells[cellI].transfer(newCells_[cellI]);
  69 }
  70
  71 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  72
  73 void polyMeshGenModifierAddCellByCell::addCell(const faceList& cellFaces)
  74 {
  75     cell c(cellFaces.size());
  76
  77     VRWGraph& pointFaces = this->pointFaces();
  78
  79     forAll(cellFaces, faceI)
  80     {
  81         const face& f = cellFaces[faceI];
  82
  83         const label pointI = f[0];
  84
  85         label fLabel(-1);
  86         forAllRow(pointFaces, pointI, pfI)
  87         {
  88             const label faceI = pointFaces(pointI, pfI);
  89
  90             if( newFaces_[faceI] == f )
  91             {
  92                 fLabel = faceI;
  93                 break;
  94             }
  95         }
  96
  97         if( fLabel == -1 )
  98         {
  99             newFaces_.append(f);
 100             c[faceI] = nFaces_;
 101             forAll(f, pI)
 102                 pointFaces.append(f[pI], nFaces_);
 103
 104             ++nFaces_;
 105         }
 106         else
 107         {
 108             c[faceI] = fLabel;
 109         }
 110     }
 111
 112     newCells_.append(c);
 113     ++nCells_;
 114 }
 115
 116 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 117
 118 } // End namespace Foam
 119
 120 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //