src/meshTools/sets/topoSets/pointSet.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2009 OpenCFD Ltd.
   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 the
  13     Free Software Foundation; either version 2 of the License, or (at your
  14     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, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 \*---------------------------------------------------------------------------*/
  26
  27 #include "pointSet.H"
  28 #include "mapPolyMesh.H"
  29 #include "polyMesh.H"
  30 #include "syncTools.H"
  31
  32 #include "addToRunTimeSelectionTable.H"
  33
  34 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  35
  36 namespace Foam
  37 {
  38
  39 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  40
  41 defineTypeNameAndDebug(pointSet, 0);
  42
  43 addToRunTimeSelectionTable(topoSet, pointSet, word);
  44 addToRunTimeSelectionTable(topoSet, pointSet, size);
  45
  46
  47 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  48
  49 pointSet::pointSet(const IOobject& obj)
  50 :
  51     topoSet(obj, typeName)
  52 {}
  53
  54
  55 pointSet::pointSet
  56 (
  57     const polyMesh& mesh,
  58     const word& name,
  59     readOption r,
  60     writeOption w
  61 )
  62 :
  63     topoSet(mesh, typeName, name, r, w)
  64 {
  65     check(mesh.nPoints());
  66 }
  67
  68
  69 pointSet::pointSet
  70 (
  71     const polyMesh& mesh,
  72     const word& name,
  73     const label size,
  74     writeOption w
  75 )
  76 :
  77     topoSet(mesh, name, size, w)
  78 {}
  79
  80
  81 pointSet::pointSet
  82 (
  83     const polyMesh& mesh,
  84     const word& name,
  85     const labelHashSet& set,
  86     writeOption w
  87 )
  88 :
  89     topoSet(mesh, name, set, w)
  90 {}
  91
  92
  93 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
  94
  95 pointSet::~pointSet()
  96 {}
  97
  98
  99 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 100
 101 void pointSet::sync(const polyMesh& mesh)
 102 {
 103     // Convert to boolList
 104
 105     boolList contents(mesh.nPoints(), false);
 106
 107     forAllConstIter(pointSet, *this, iter)
 108     {
 109         contents[iter.key()] = true;
 110     }
 111     syncTools::syncPointList
 112     (
 113         mesh,
 114         contents,
 115         orEqOp<bool>(),
 116         false,          // null value
 117         false           // no separation
 118     );
 119
 120     // Convert back to labelHashSet
 121
 122     labelHashSet newContents(size());
 123
 124     forAll(contents, pointI)
 125     {
 126         if (contents[pointI])
 127         {
 128             newContents.insert(pointI);
 129         }
 130     }
 131
 132     transfer(newContents);
 133 }
 134
 135
 136 label pointSet::maxSize(const polyMesh& mesh) const
 137 {
 138     return mesh.nPoints();
 139 }
 140
 141
 142 void pointSet::updateMesh(const mapPolyMesh& morphMap)
 143 {
 144     updateLabels(morphMap.reversePointMap());
 145 }
 146
 147
 148 void pointSet::writeDebug
 149 (
 150     Ostream& os,
 151     const primitiveMesh& mesh,
 152     const label maxLen
 153 ) const
 154 {
 155     topoSet::writeDebug(os, mesh.points(), maxLen);
 156 }
 157
 158 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 159
 160 } // End namespace Foam
 161
 162 // ************************************************************************* //