src/conversion/meshReader/meshReader.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2008 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 "meshReader.H"
  28 #include "Time.H"
  29 #include "polyMesh.H"
  30 #include "faceSet.H"
  31 #include "emptyPolyPatch.H"
  32 #include "cellModeller.H"
  33 #include "demandDrivenData.H"
  34
  35 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  36
  37 const Foam::cellModel* Foam::meshReader::unknownModel = Foam::cellModeller::
  38 lookup
  39 (
  40     "unknown"
  41 );
  42
  43 const Foam::cellModel* Foam::meshReader::tetModel = Foam::cellModeller::
  44 lookup
  45 (
  46     "tet"
  47 );
  48
  49 const Foam::cellModel* Foam::meshReader::pyrModel = Foam::cellModeller::
  50 lookup
  51 (
  52     "pyr"
  53 );
  54
  55 const Foam::cellModel* Foam::meshReader::prismModel = Foam::cellModeller::
  56 lookup
  57 (
  58     "prism"
  59 );
  60
  61 const Foam::cellModel* Foam::meshReader::hexModel = Foam::cellModeller::
  62 lookup
  63 (
  64     "hex"
  65 );
  66
  67
  68 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  69
  70 void Foam::meshReader::addCellZones(polyMesh& mesh) const
  71 {
  72     cellTable_.addCellZones(mesh, cellTableId_);
  73     warnDuplicates("cellZones", mesh.cellZones().names());
  74 }
  75
  76
  77 void Foam::meshReader::addFaceZones(polyMesh& mesh) const
  78 {
  79     label nZone = monitoringSets_.size();
  80     mesh.faceZones().setSize(nZone);
  81
  82     if (!nZone)
  83     {
  84         return;
  85     }
  86
  87     nZone = 0;
  88     for
  89     (
  90         HashTable<List<label>, word, string::hash>::const_iterator
  91         iter = monitoringSets_.begin();
  92         iter != monitoringSets_.end();
  93         ++iter
  94     )
  95     {
  96         Info<< "faceZone " << nZone
  97             << " (size: " << iter().size() << ") name: "
  98             << iter.key() << endl;
  99
 100         mesh.faceZones().set
 101         (
 102             nZone,
 103             new faceZone
 104             (
 105                 iter.key(),
 106                 iter(),
 107                 List<bool>(iter().size(), false),
 108                 nZone,
 109                 mesh.faceZones()
 110             )
 111         );
 112
 113         nZone++;
 114     }
 115     mesh.faceZones().writeOpt() = IOobject::AUTO_WRITE;
 116     warnDuplicates("faceZones", mesh.faceZones().names());
 117 }
 118
 119
 120 Foam::autoPtr<Foam::polyMesh> Foam::meshReader::mesh
 121 (
 122     const objectRegistry& registry
 123 )
 124 {
 125     readGeometry();
 126
 127     Info<< "Creating a polyMesh" << endl;
 128     createPolyCells();
 129
 130     Info<< "Number of internal faces: " << nInternalFaces_ << endl;
 131
 132     createPolyBoundary();
 133     clearExtraStorage();
 134
 135     autoPtr<polyMesh> mesh
 136     (
 137         new polyMesh
 138         (
 139             IOobject
 140             (
 141                 polyMesh::defaultRegion,
 142                 "constant",
 143                 registry
 144             ),
 145             points(),
 146             meshFaces_,
 147             cellPolys_
 148         )
 149     );
 150
 151     // adding patches also checks the mesh
 152     mesh().addPatches(polyBoundaryPatches(mesh));
 153
 154     warnDuplicates("boundaries", mesh().boundaryMesh().names());
 155
 156     addCellZones(mesh());
 157     addFaceZones(mesh());
 158
 159     return mesh;
 160 }
 161
 162
 163 void Foam::meshReader::writeMesh
 164 (
 165     const polyMesh& mesh,
 166     IOstream::streamFormat fmt
 167 ) const
 168 {
 169     fileName meshDir = mesh.objectRegistry::path()/mesh.meshDir();
 170
 171     // remove some directories and files - this should be easier
 172     mesh.removeFiles(mesh.instance());
 173     if (dir(meshDir/"sets"))
 174     {
 175         rmDir(meshDir/"sets");
 176     }
 177
 178     Info<< "Writing polyMesh" << endl;
 179     mesh.writeObject
 180     (
 181         fmt,
 182         IOstream::currentVersion,
 183         IOstream::UNCOMPRESSED
 184     );
 185     writeAux(mesh);
 186 }
 187
 188
 189 void Foam::meshReader::clearExtraStorage()
 190 {
 191     cellFaces_.clear();
 192     baffleFaces_.clear();
 193     boundaryIds_.clear();
 194     baffleIds_.clear();
 195
 196     deleteDemandDrivenData(pointCellsPtr_);
 197 }
 198
 199
 200 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 201
 202 Foam::meshReader::meshReader
 203 (
 204     const fileName& fileOrPrefix,
 205     const scalar scaleFactor
 206 )
 207     :
 208     pointCellsPtr_(NULL),
 209     nInternalFaces_(0),
 210     patchStarts_(0),
 211     patchSizes_(0),
 212     interfaces_(0),
 213     baffleIds_(0),
 214     meshFaces_(0),
 215     cellPolys_(0),
 216     geometryFile_(fileOrPrefix),
 217     scaleFactor_(scaleFactor),
 218     points_(0),
 219     origCellId_(0),
 220     boundaryIds_(0),
 221     patchTypes_(0),
 222     patchNames_(0),
 223     patchPhysicalTypes_(0),
 224     cellFaces_(0),
 225     baffleFaces_(0),
 226     cellTableId_(0),
 227     cellTable_()
 228 {}
 229
 230
 231 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 232
 233 Foam::meshReader::~meshReader()
 234 {
 235     deleteDemandDrivenData(pointCellsPtr_);
 236 }
 237
 238
 239 // ************************************************************************* //