src/mesh/cfMesh/cartesianMesh/cartesianMeshExtractor/cartesianMeshExtractorDecomposeSplitHexes.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 "demandDrivenData.H"
  30 #include "decomposeFaces.H"
  31 #include "decomposeCells.H"
  32 #include "hexMatcher.H"
  33
  34 //#define DEBUGMesh
  35
  36 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  37
  38 namespace Foam
  39 {
  40
  41 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  42
  43 void cartesianMeshExtractor::decomposeSplitHexesIntoTetsAndPyramids()
  44 {
  45     if( !decomposeSplitHexes_ ) return;
  46
  47     Info << "Decomposing split-hex cells" << endl;
  48
  49     const faceListPMG& faces = mesh_.faces();
  50
  51     //- decompose faces which have more than 4 vertices
  52     boolList decompose(faces.size(), false);
  53
  54     label nDecomposed(0);
  55     forAll(faces, faceI)
  56     {
  57         if( faces[faceI].size() > 4 )
  58         {
  59             ++nDecomposed;
  60
  61             decompose[faceI] = true;
  62         }
  63     }
  64
  65     reduce(nDecomposed, sumOp<label>());
  66
  67     Info << "Decomposing " << nDecomposed
  68         << " faces with more than 4 vertices" << endl;
  69
  70     if( nDecomposed != 0 )
  71     {
  72         //- decompose marked faces into triangles
  73         decomposeFaces(mesh_).decomposeMeshFaces(decompose);
  74     }
  75
  76     //- decompose cells with 24 faces
  77     const cellListPMG& cells = mesh_.cells();
  78     decompose.setSize(cells.size());
  79     decompose = false;
  80
  81     hexMatcher hex;
  82     forAll(cells, cellI)
  83     {
  84         if( !hex.matchShape(true, faces, mesh_.owner(), cellI, cells[cellI]) )
  85         {
  86             ++nDecomposed;
  87             decompose[cellI] = true;
  88         }
  89     }
  90
  91     reduce(nDecomposed, sumOp<label>());
  92
  93     Info << "Decomposing " << nDecomposed
  94         << " cells into tetrahedra and pyramids" << endl;
  95
  96     if( nDecomposed )
  97     {
  98         //- decompose marked cells into tets and pyramids
  99         decomposeCells dc(mesh_);
 100         dc.decomposeMesh(decompose);
 101     }
 102
 103     Info << "Finished decomposing split-hex cells" << endl;
 104 }
 105
 106 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 107
 108 } // End namespace Foam
 109
 110 // ************************************************************************* //