src/conversion/ensight/part/ensightPart.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | foam-extend: Open Source CFD
   4    \\    /   O peration     | Version:     3.2
   5     \\  /    A nd           | Web:         http://www.foam-extend.org
   6      \\/     M anipulation  | For copyright notice see file Copyright
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of foam-extend.
  10
  11     foam-extend 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     foam-extend is distributed in the hope that it will be useful, but
  17     WITHOUT ANY WARRANTY; without even the implied warranty of
  18     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  19     General Public License for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 \*----------------------------------------------------------------------------*/
  25
  26 #include "ensightPart.H"
  27 #include "addToRunTimeSelectionTable.H"
  28 #include "dictionary.H"
  29 #include "ListOps.H"
  30
  31 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  32 namespace Foam
  33 {
  34    defineTypeNameAndDebug(ensightPart, 0);
  35    defineTemplateTypeNameAndDebug(IOPtrList<ensightPart>, 0);
  36    defineRunTimeSelectionTable(ensightPart, istream);
  37 }
  38
  39 Foam::List<Foam::word> Foam::ensightPart::elemTypes_(0);
  40
  41 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  42
  43 bool Foam::ensightPart::isFieldDefined(const List<scalar>& field) const
  44 {
  45     forAll(elemLists_, elemI)
  46     {
  47         const labelList& idList = elemLists_[elemI];
  48
  49         forAll(idList, i)
  50         {
  51             label id = idList[i];
  52
  53 #           ifndef Intel
  54             if (id >= field.size() || std::isnan(field[id]))
  55 #           else
  56             if (id >= field.size())
  57 #           endif
  58             {
  59                 return false;
  60             }
  61         }
  62     }
  63     return true;
  64 }
  65
  66
  67 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  68
  69 Foam::ensightPart::ensightPart
  70 ()
  71 :
  72     number_(0),
  73     name_(""),
  74     elemLists_(0),
  75     offset_(0),
  76     size_(0),
  77     isCellData_(true),
  78     matId_(0),
  79     meshPtr_(0)
  80 {}
  81
  82
  83 Foam::ensightPart::ensightPart
  84 (
  85     label partNumber,
  86     const string& partDescription
  87 )
  88 :
  89     number_(partNumber),
  90     name_(partDescription),
  91     elemLists_(0),
  92     offset_(0),
  93     size_(0),
  94     isCellData_(true),
  95     matId_(0),
  96     meshPtr_(0)
  97 {}
  98
  99
 100 Foam::ensightPart::ensightPart
 101 (
 102     label partNumber,
 103     const string& partDescription,
 104     const polyMesh& pMesh
 105 )
 106 :
 107     number_(partNumber),
 108     name_(partDescription),
 109     elemLists_(0),
 110     offset_(0),
 111     size_(0),
 112     isCellData_(true),
 113     matId_(0),
 114     meshPtr_(&pMesh)
 115 {}
 116
 117
 118 Foam::ensightPart::ensightPart(const ensightPart& part)
 119 :
 120     number_(part.number_),
 121     name_(part.name_),
 122     elemLists_(part.elemLists_),
 123     offset_(part.offset_),
 124     size_(part.size_),
 125     isCellData_(part.isCellData_),
 126     matId_(part.matId_),
 127     meshPtr_(part.meshPtr_)
 128 {}
 129
 130
 131 // * * * * * * * * * * * * * * * * Selectors * * * * * * * * * * * * * * * * //
 132
 133 Foam::autoPtr<Foam::ensightPart> Foam::ensightPart::New(Istream& is)
 134 {
 135     word partType(is);
 136
 137     istreamConstructorTable::iterator cstrIter =
 138         istreamConstructorTablePtr_->find(partType);
 139
 140     if (cstrIter == istreamConstructorTablePtr_->end())
 141     {
 142         FatalIOErrorIn
 143         (
 144             "ensightPart::New(Istream&)",
 145             is
 146         )   << "unknown ensightPart type " << partType << endl << endl
 147             << "Valid ensightPart types are :" << endl
 148             << istreamConstructorTablePtr_->sortedToc()
 149             << exit(FatalIOError);
 150     }
 151
 152     return autoPtr<ensightPart>(cstrIter()(is));
 153 }
 154
 155
 156 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 157
 158 Foam::ensightPart::~ensightPart()
 159 {}
 160
 161
 162 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 163
 164 void Foam::ensightPart::reconstruct(Istream& is)
 165 {
 166     dictionary dict(is);
 167     dict.lookup("id") >> number_;
 168     dict.lookup("name") >> name_;
 169     dict.readIfPresent("offset", offset_);
 170
 171     // populate elemLists_
 172     elemLists_.setSize(elementTypes().size());
 173
 174     forAll(elementTypes(), elemI)
 175     {
 176         word key(elementTypes()[elemI]);
 177
 178         elemLists_[elemI].clear();
 179         dict.readIfPresent(key, elemLists_[elemI]);
 180
 181         size_ += elemLists_[elemI].size();
 182     }
 183
 184     is.check("ensightPart::reconstruct(Istream&)");
 185 }
 186
 187
 188 void Foam::ensightPart::renumber(labelList const& origId)
 189 {
 190     // transform to global values first
 191     if (offset_)
 192     {
 193         forAll(elemLists_, elemI)
 194         {
 195             labelList& idList = elemLists_[elemI];
 196             forAll(idList, i)
 197             {
 198                 idList[i] += offset_;
 199             }
 200         }
 201
 202         offset_ = 0;
 203     }
 204
 205     if (origId.size())
 206     {
 207         forAll(elemLists_, elemI)
 208         {
 209             inplaceRenumber(origId, elemLists_[elemI]);
 210         }
 211     }
 212 }
 213
 214
 215 // ************************************************************************* //