src/mesh/cfMesh/meshLibrary/cartesianMesh/cartesianMeshExtractor/cartesianMeshExtractor.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | cfMesh: A library for mesh generation
   4    \\    /   O peration     |
   5     \\  /    A nd           | Author: Franjo Juretic (franjo.juretic@c-fields.com)
   6      \\/     M anipulation  | Copyright (C) Creative Fields, Ltd.
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of cfMesh.
  10
  11     cfMesh 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     cfMesh 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 cfMesh.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 Description
  25
  26 \*---------------------------------------------------------------------------*/
  27
  28 #include "cartesianMeshExtractor.H"
  29 #include "meshOctree.H"
  30
  31 // #define DEBUGSearch
  32
  33 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  34
  35 namespace Foam
  36 {
  37
  38 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  39
  40 void cartesianMeshExtractor::clearOut()
  41 {
  42     deleteDemandDrivenData(leafCellLabelPtr_);
  43 }
  44
  45 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  46
  47 // Construct from octree and mesh data
  48 cartesianMeshExtractor::cartesianMeshExtractor
  49 (
  50     meshOctree& octree,
  51     const IOdictionary& meshDict,
  52     polyMeshGen& mesh
  53 )
  54 :
  55     octreeCheck_(octree, meshDict, false),
  56     mesh_(mesh),
  57     decomposeSplitHexes_(false),
  58     leafCellLabelPtr_(new labelList(octree.numberOfLeaves(), -1))
  59 {
  60 }
  61
  62 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
  63
  64 cartesianMeshExtractor::~cartesianMeshExtractor()
  65 {
  66     clearOut();
  67 }
  68
  69 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  70
  71 void cartesianMeshExtractor::decomposeSplitHexes()
  72 {
  73     decomposeSplitHexes_ = true;
  74 }
  75
  76 void cartesianMeshExtractor::createMesh()
  77 {
  78     Info << "Extracting polyMesh" << endl;
  79
  80     //- create points and pointLeaves addressing
  81     createPointsAndAddressing();
  82
  83     //- create the mesh
  84     createPolyMesh();
  85
  86     //- decompose split-hex cells into tetrahedra and pyramids
  87     decomposeSplitHexesIntoTetsAndPyramids();
  88
  89     //- remove unused vertices
  90     polyMeshGenModifier(mesh_).removeUnusedVertices();
  91
  92     Info << "Mesh has :" << nl
  93         << mesh_.points().size() << " vertices " << nl
  94         << mesh_.faces().size() << " faces" << nl
  95         << mesh_.cells().size() << " cells" << endl;
  96
  97     if( Pstream::parRun() )
  98     {
  99         label nCells = mesh_.cells().size();
 100         reduce(nCells, sumOp<label>());
 101         Info << "Total number of cells " << nCells << endl;
 102     }
 103     if( mesh_.cells().size() == 0 )
 104     {
 105         FatalErrorIn
 106         (
 107             "void cartesianMeshExtractor::createMesh()"
 108         ) << "There are no cells in the mesh!"
 109         << nl << "The reasons for this can be fwofold:"
 110         << nl << "1. Inadequate mesh resolution."
 111         << nl << "2. You maxCellSize is a multiplier of the domain length."
 112         << " This can be reolved by reducing the maxCellSize by a fraction."
 113         << "i.e. 2.49999 instead of 2.5." << exit(FatalError);
 114     }
 115
 116     Info << "Finished extracting polyMesh" << endl;
 117 }
 118
 119 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 120
 121 } // End namespace Foam
 122
 123 // ************************************************************************* //