FreeFOAM The Cross-Platform CFD Toolkit
Hosted by SourceForge:
Get FreeFOAM at SourceForge.net.
            Fast, secure and Free Open Source software downloads

setToPointZone.C

Go to the documentation of this file.
00001 /*---------------------------------------------------------------------------*\
00002   =========                 |
00003   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
00004    \\    /   O peration     |
00005     \\  /    A nd           | Copyright (C) 1991-2010 OpenCFD Ltd.
00006      \\/     M anipulation  |
00007 -------------------------------------------------------------------------------
00008 License
00009     This file is part of OpenFOAM.
00010 
00011     OpenFOAM is free software: you can redistribute it and/or modify it
00012     under the terms of the GNU General Public License as published by
00013     the Free Software Foundation, either version 3 of the License, or
00014     (at your option) any later version.
00015 
00016     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
00017     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
00018     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
00019     for more details.
00020 
00021     You should have received a copy of the GNU General Public License
00022     along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
00023 
00024 \*---------------------------------------------------------------------------*/
00025 
00026 #include "setToPointZone.H"
00027 #include <OpenFOAM/polyMesh.H>
00028 #include <meshTools/pointZoneSet.H>
00029 
00030 #include <OpenFOAM/addToRunTimeSelectionTable.H>
00031 
00032 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
00033 
00034 namespace Foam
00035 {
00036 
00037 defineTypeNameAndDebug(setToPointZone, 0);
00038 
00039 addToRunTimeSelectionTable(topoSetSource, setToPointZone, word);
00040 
00041 addToRunTimeSelectionTable(topoSetSource, setToPointZone, istream);
00042 
00043 }
00044 
00045 
00046 Foam::topoSetSource::addToUsageTable Foam::setToPointZone::usage_
00047 (
00048     setToPointZone::typeName,
00049     "\n    Usage: setToPointZone <pointSet>\n\n"
00050     "    Select all points in the pointSet.\n\n"
00051 );
00052 
00053 
00054 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
00055 
00056 // Construct from components
00057 Foam::setToPointZone::setToPointZone
00058 (
00059     const polyMesh& mesh,
00060     const word& setName
00061 )
00062 :
00063     topoSetSource(mesh),
00064     setName_(setName)
00065 {}
00066 
00067 
00068 // Construct from dictionary
00069 Foam::setToPointZone::setToPointZone
00070 (
00071     const polyMesh& mesh,
00072     const dictionary& dict
00073 )
00074 :
00075     topoSetSource(mesh),
00076     setName_(dict.lookup("set"))
00077 {}
00078 
00079 
00080 // Construct from Istream
00081 Foam::setToPointZone::setToPointZone
00082 (
00083     const polyMesh& mesh,
00084     Istream& is
00085 )
00086 :
00087     topoSetSource(mesh),
00088     setName_(checkIs(is))
00089 {}
00090 
00091 
00092 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
00093 
00094 Foam::setToPointZone::~setToPointZone()
00095 {}
00096 
00097 
00098 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
00099 
00100 void Foam::setToPointZone::applyToSet
00101 (
00102     const topoSetSource::setAction action,
00103     topoSet& set
00104 ) const
00105 {
00106     if (!isA<pointZoneSet>(set))
00107     {
00108         WarningIn
00109         (
00110             "setToPointZone::applyToSet(const topoSetSource::setAction"
00111             ", topoSet"
00112         )   << "Operation only allowed on a pointZoneSet." << endl;
00113     }
00114     else
00115     {
00116         pointZoneSet& fzSet = refCast<pointZoneSet>(set);
00117 
00118         if ((action == topoSetSource::NEW) || (action == topoSetSource::ADD))
00119         {
00120             Info<< "    Adding all points from pointSet " << setName_
00121                 << " ..." << endl;
00122 
00123             // Load the sets
00124             pointSet fSet(mesh_, setName_);
00125 
00126             // Start off from copy
00127             DynamicList<label> newAddressing(fzSet.addressing());
00128 
00129             forAllConstIter(pointSet, fSet, iter)
00130             {
00131                 label pointI = iter.key();
00132 
00133                 if (!fzSet.found(pointI))
00134                 {
00135                     newAddressing.append(pointI);
00136                 }
00137             }
00138 
00139             fzSet.addressing().transfer(newAddressing);
00140             fzSet.updateSet();
00141         }
00142         else if (action == topoSetSource::DELETE)
00143         {
00144             Info<< "    Removing all points from pointSet " << setName_
00145                 << " ..." << endl;
00146 
00147             // Load the set
00148             pointSet loadedSet(mesh_, setName_);
00149 
00150             // Start off empty
00151             DynamicList<label> newAddressing(fzSet.addressing().size());
00152 
00153             forAll(fzSet.addressing(), i)
00154             {
00155                 if (!loadedSet.found(fzSet.addressing()[i]))
00156                 {
00157                     newAddressing.append(fzSet.addressing()[i]);
00158                 }
00159             }
00160             fzSet.addressing().transfer(newAddressing);
00161             fzSet.updateSet();
00162         }
00163     }
00164 }
00165 
00166 
00167 // ************************ vim: set sw=4 sts=4 et: ************************ //
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines