src/mesh/cfMesh/utilities/boundaryLayers/triangulateNonPlanarBaseFaces/triangulateNonPlanarBaseFaces.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 "triangulateNonPlanarBaseFaces.H"
  29 #include "dictionary.H"
  30 #include "demandDrivenData.H"
  31
  32 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  33
  34 namespace Foam
  35 {
  36
  37 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  38
  39 triangulateNonPlanarBaseFaces::triangulateNonPlanarBaseFaces
  40 (
  41     polyMeshGen& mesh
  42 )
  43 :
  44     mesh_(mesh),
  45     invertedCell_(mesh_.cells().size(), false),
  46     decomposeFace_(mesh_.faces().size(), false),
  47     tol_(0.5)
  48 {}
  49
  50 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
  51
  52 triangulateNonPlanarBaseFaces::~triangulateNonPlanarBaseFaces()
  53 {}
  54
  55 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  56
  57 void triangulateNonPlanarBaseFaces::setRelativeTolerance(const scalar tol)
  58 {
  59     tol_ = tol;
  60 }
  61
  62 void triangulateNonPlanarBaseFaces::triangulateLayers()
  63 {
  64     if( findNonPlanarBoundaryFaces() )
  65     {
  66         Info << "Decomposing twisted boundary faces" << endl;
  67
  68         decomposeBoundaryFaces();
  69
  70         decomposeCellsIntoPyramids();
  71     }
  72     else
  73     {
  74         Info << "All boundary faces are flat" << endl;
  75     }
  76 }
  77
  78 void triangulateNonPlanarBaseFaces::readSettings
  79 (
  80     const dictionary& meshDict,
  81     triangulateNonPlanarBaseFaces& triangulator
  82 )
  83 {
  84     if( meshDict.found("boundaryLayers") )
  85     {
  86         const dictionary& layersDict = meshDict.subDict("boundaryLayers");
  87
  88         if( layersDict.found("optimisationParameters") )
  89         {
  90             const dictionary& optLayerDict =
  91                 layersDict.subDict("optimisationParameters");
  92
  93             if( optLayerDict.found("relFlatnessTol") )
  94             {
  95                 const scalar relTol =
  96                     readScalar(optLayerDict.lookup("relFlatnessTol"));
  97
  98                 triangulator.setRelativeTolerance(relTol);
  99             }
 100         }
 101     }
 102 }
 103
 104 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 105
 106 } // End namespace Foam
 107
 108 // ************************************************************************* //