src/mesh/cfMesh/meshLibrary/utilities/meshes/triSurf/triSurfFacets.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 "triSurfFacets.H"
  29 #include "pointIOField.H"
  30 #include "IOobjectList.H"
  31 #include "pointSet.H"
  32
  33 namespace Foam
  34 {
  35
  36 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  37
  38 triSurfFacets::triSurfFacets()
  39 :
  40     triangles_(),
  41     patches_(),
  42     facetSubsets_()
  43 {}
  44
  45 triSurfFacets::triSurfFacets(const LongList<labelledTri>& triangles)
  46 :
  47     triangles_(triangles),
  48     patches_(1),
  49     facetSubsets_()
  50 {
  51     forAll(triangles_, triI)
  52         triangles_[triI].region() = 0;
  53
  54     patches_[0].name() = "patch";
  55 }
  56
  57 triSurfFacets::triSurfFacets
  58 (
  59     const LongList<labelledTri>& triangles,
  60     const geometricSurfacePatchList& patches
  61 )
  62 :
  63     triangles_(triangles),
  64     patches_(patches),
  65     facetSubsets_()
  66 {}
  67
  68 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  69 // Destructor
  70 triSurfFacets::~triSurfFacets()
  71 {}
  72
  73 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  74
  75 label triSurfFacets::addFacetSubset(const word& subsetName)
  76 {
  77     label id = facetSubsetIndex(subsetName);
  78     if( id >= 0 )
  79     {
  80         Warning << "Point subset " << subsetName << " already exists!" << endl;
  81         return id;
  82     }
  83
  84     id = 0;
  85     forAllConstIter(Map<meshSubset>, facetSubsets_, it)
  86         id = Foam::max(id, it.key()+1);
  87
  88     facetSubsets_.insert
  89     (
  90         id,
  91         meshSubset(subsetName, meshSubset::FACESUBSET)
  92     );
  93
  94     return id;
  95 }
  96
  97 void triSurfFacets::removeFacetSubset(const label subsetID)
  98 {
  99     if( facetSubsets_.find(subsetID) == facetSubsets_.end() )
 100         return;
 101
 102     facetSubsets_.erase(subsetID);
 103 }
 104
 105 word triSurfFacets::facetSubsetName(const label subsetID) const
 106 {
 107     Map<meshSubset>::const_iterator it = facetSubsets_.find(subsetID);
 108     if( it == facetSubsets_.end() )
 109     {
 110         Warning << "Subset " << subsetID << " is not a facet subset" << endl;
 111         return word();
 112     }
 113
 114     return it().name();
 115 }
 116
 117 label triSurfFacets::facetSubsetIndex(const word& subsetName) const
 118 {
 119     forAllConstIter(Map<meshSubset>, facetSubsets_, it)
 120     {
 121         if( it().name() == subsetName )
 122             return it.key();
 123     }
 124
 125     return -1;
 126 }
 127
 128 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 129
 130 } // End namespace Foam
 131
 132 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //