src/surfMesh/surfaceFormats/tri/TRIsurfaceFormatCore.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 "TRIsurfaceFormat.H"
  28 #include "IFstream.H"
  29 #include "IOmanip.H"
  30 #include "IStringStream.H"
  31 #include "Map.H"
  32
  33 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  34
  35 Foam::fileFormats::TRIsurfaceFormatCore::TRIsurfaceFormatCore
  36 (
  37     const fileName& filename
  38 )
  39 :
  40     sorted_(true),
  41     points_(0),
  42     zoneIds_(0),
  43     sizes_(0)
  44 {
  45     read(filename);
  46 }
  47
  48
  49 // * * * * * * * * * * * * * * * * Destructors * * * * * * * * * * * * * * * //
  50
  51 Foam::fileFormats::TRIsurfaceFormatCore::~TRIsurfaceFormatCore()
  52 {}
  53
  54
  55 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  56
  57 bool Foam::fileFormats::TRIsurfaceFormatCore::read
  58 (
  59     const fileName& filename
  60 )
  61 {
  62     this->clear();
  63     sorted_ = true;
  64
  65     IFstream is(filename);
  66     if (!is.good())
  67     {
  68         FatalErrorIn
  69         (
  70             "fileFormats::TRIsurfaceFormatCore::read(const fileName&)"
  71         )
  72             << "Cannot read file " << filename
  73             << exit(FatalError);
  74     }
  75
  76     // uses similar structure as STL, just some points
  77     // the rest of the reader resembles the STL binary reader
  78     DynamicList<point> dynPoints;
  79     DynamicList<label> dynZones;
  80     DynamicList<label> dynSizes;
  81     HashTable<label>   lookup;
  82
  83     // place faces without a group in zone0
  84     label zoneI = 0;
  85     dynSizes.append(zoneI);
  86     lookup.insert("zoneI", zoneI);
  87
  88     while (is.good())
  89     {
  90         string line = this->getLineNoComment(is);
  91
  92         // handle continuations ?
  93         //          if (line[line.size()-1] == '\\')
  94         //          {
  95         //              line.substr(0, line.size()-1);
  96         //              line += this->getLineNoComment(is);
  97         //          }
  98
  99         IStringStream lineStream(line);
 100
 101         point p
 102         (
 103             readScalar(lineStream),
 104             readScalar(lineStream),
 105             readScalar(lineStream)
 106         );
 107
 108         if (!lineStream) break;
 109
 110         dynPoints.append(p);
 111         dynPoints.append
 112         (
 113             point
 114             (
 115                 readScalar(lineStream),
 116                 readScalar(lineStream),
 117                 readScalar(lineStream)
 118             )
 119         );
 120         dynPoints.append
 121         (
 122             point
 123             (
 124                 readScalar(lineStream),
 125                 readScalar(lineStream),
 126                 readScalar(lineStream)
 127             )
 128         );
 129
 130         // zone/colour in .tri file starts with 0x. Skip.
 131         // ie, instead of having 0xFF, skip 0 and leave xFF to
 132         // get read as a word and name it "zoneFF"
 133
 134         char zero;
 135         lineStream >> zero;
 136
 137         word rawName(lineStream);
 138         word name("zone" + rawName(1, rawName.size()-1));
 139
 140         HashTable<label>::const_iterator fnd = lookup.find(name);
 141         if (fnd != lookup.end())
 142         {
 143             if (zoneI != fnd())
 144             {
 145                 // group appeared out of order
 146                 sorted_ = false;
 147             }
 148             zoneI = fnd();
 149         }
 150         else
 151         {
 152             zoneI = dynSizes.size();
 153             lookup.insert(name, zoneI);
 154             dynSizes.append(0);
 155         }
 156
 157         dynZones.append(zoneI);
 158         dynSizes[zoneI]++;
 159     }
 160
 161     // skip empty groups
 162     label nZone = 0;
 163     forAll(dynSizes, zoneI)
 164     {
 165         if (dynSizes[zoneI])
 166         {
 167             if (nZone != zoneI)
 168             {
 169                 dynSizes[nZone] = dynSizes[zoneI];
 170             }
 171             nZone++;
 172         }
 173     }
 174     // truncate addressed size
 175     dynSizes.setCapacity(nZone);
 176
 177     // transfer to normal lists
 178     points_.transfer(dynPoints);
 179     zoneIds_.transfer(dynZones);
 180     sizes_.transfer(dynSizes);
 181
 182     return true;
 183 }
 184
 185
 186 // ************************************************************************* //