src/mesh/cfMesh/meshLibrary/utilities/triSurfaceTools/triSurfacePartitioner/triSurfacePartitioner.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 "triSurfacePartitioner.H"
  29 #include "demandDrivenData.H"
  30
  31 # ifdef DEBUGPartitioner
  32 #include <sstream>
  33 # endif
  34
  35 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  36
  37 namespace Foam
  38 {
  39
  40 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  41
  42 triSurfacePartitioner::triSurfacePartitioner
  43 (
  44     const triSurf& surface
  45 )
  46 :
  47     surface_(surface),
  48     corners_(),
  49     cornerPatches_(),
  50     patchPatches_(surface.patches().size()),
  51     edgeGroups_(),
  52     edgeGroupEdgeGroups_(),
  53     patchesEdgeGroups_(),
  54     edgeGroupsCorners_()
  55 {
  56     calculatePatchAddressing();
  57 }
  58
  59 triSurfacePartitioner::~triSurfacePartitioner()
  60 {}
  61
  62 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  63
  64 const labelList& triSurfacePartitioner::corners() const
  65 {
  66     return corners_;
  67 }
  68
  69 const List<DynList<label> >& triSurfacePartitioner::cornerPatches() const
  70 {
  71     return cornerPatches_;
  72 }
  73
  74 const List<labelHashSet>& triSurfacePartitioner::patchPatches() const
  75 {
  76     return patchPatches_;
  77 }
  78
  79 const labelList& triSurfacePartitioner::edgeGroups() const
  80 {
  81     return edgeGroups_;
  82 }
  83
  84 const List<labelHashSet>& triSurfacePartitioner::edgeGroupEdgeGroups() const
  85 {
  86     return edgeGroupEdgeGroups_;
  87 }
  88
  89 void triSurfacePartitioner::edgeGroupsSharedByPatches
  90 (
  91     const label patch1,
  92     const label patch2,
  93     DynList<label>& edgeGroups
  94 ) const
  95 {
  96     edgeGroups.clear();
  97
  98     std::pair<label, label> pp
  99     (
 100         Foam::min(patch1, patch2),
 101         Foam::max(patch1, patch2)
 102     );
 103
 104     std::map<std::pair<label, label>, labelHashSet>::const_iterator it =
 105         patchesEdgeGroups_.find(pp);
 106
 107     if( it != patchesEdgeGroups_.end() )
 108     {
 109         const labelHashSet& eGroups = it->second;
 110
 111         forAllConstIter(labelHashSet, eGroups, it)
 112             edgeGroups.append(it.key());
 113     }
 114 }
 115
 116 void triSurfacePartitioner::cornersSharedByEdgeGroups
 117 (
 118     const label edgeGroup1,
 119     const label edgeGroup2,
 120     DynList<label>& corners
 121 ) const
 122 {
 123     corners.clear();
 124
 125     std::pair<label, label> ep
 126     (
 127         Foam::min(edgeGroup1, edgeGroup2),
 128         Foam::max(edgeGroup1, edgeGroup2)
 129     );
 130
 131     std::map<std::pair<label, label>, labelHashSet>::const_iterator it =
 132         edgeGroupsCorners_.find(ep);
 133
 134     if( it != edgeGroupsCorners_.end() )
 135     {
 136         const labelHashSet& corn = it->second;
 137
 138         forAllConstIter(labelHashSet, corn, it)
 139             corners.append(it.key());
 140     }
 141 }
 142
 143 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 144
 145 } // End namespace Foam
 146
 147 // ************************************************************************* //