src/mesh/cfMesh/utilities/boundaryLayers/triangulateNonPlanarBaseFaces/triangulateNonPlanarBaseFacesFunctions.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 "decomposeFaces.H"
  30 #include "decomposeCells.H"
  31 #include "helperFunctions.H"
  32 #include "meshSurfacePartitioner.H"
  33 #include "triangle.H"
  34
  35 # ifdef USE_OMP
  36 #include <omp.h>
  37 # endif
  38
  39 //#define DEBUGLayer
  40
  41 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  42
  43 namespace Foam
  44 {
  45
  46 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  47
  48 bool triangulateNonPlanarBaseFaces::findNonPlanarBoundaryFaces()
  49 {
  50     const pointFieldPMG& points = mesh_.points();
  51     const label nInternalFaces = mesh_.nInternalFaces();
  52
  53     meshSurfacePartitioner mPart(mesh_);
  54     const meshSurfaceEngine& mse = mPart.surfaceEngine();
  55     const labelList& faceOwner = mse.faceOwners();
  56     const faceList::subList& bFaces = mse.boundaryFaces();
  57
  58     bool hasInvalid(false);
  59
  60     # ifdef USE_OMP
  61     # pragma omp parallel for schedule(dynamic, 50)
  62     # endif
  63     forAll(bFaces, bfI)
  64     {
  65         const face& bf = bFaces[bfI];
  66
  67         //- triangle shall not be decomposed, they are flat
  68         if( bf.size() == 3 )
  69             continue;
  70
  71         //- calculate min face diagonal
  72         scalar minDist(VGREAT);
  73         const point c = bf.centre(points);
  74         forAll(bf, pI)
  75         {
  76             minDist = Foam::min(minDist, Foam::mag(c - points[bf[pI]]));
  77         }
  78
  79         forAll(bf, eI)
  80         {
  81             triangle<point, point> tri
  82             (
  83                 points[bf[eI]],
  84                 points[bf.nextLabel(eI)],
  85                 c
  86             );
  87
  88             const point triCentre = tri.centre();
  89             vector n = tri.normal();
  90             n /= (mag(n) + VSMALL);
  91
  92             forAll(bf, pI)
  93             {
  94                 const scalar d = (points[bf[pI]] - triCentre) & n;
  95
  96                 if( d > tol_ * minDist )
  97                 {
  98                     invertedCell_[faceOwner[bfI]] = true;
  99
 100                     decomposeFace_[nInternalFaces+bfI] = true;
 101                     hasInvalid = true;
 102                 }
 103             }
 104         }
 105     }
 106
 107     reduce(hasInvalid, maxOp<bool>());
 108
 109     return hasInvalid;
 110 }
 111
 112 void triangulateNonPlanarBaseFaces::decomposeBoundaryFaces()
 113 {
 114     //- decompose base faces into triangles
 115     decomposeFaces triangulator(mesh_);
 116     triangulator.decomposeMeshFaces(decomposeFace_);
 117     const VRWGraph& newFacesFromFace = triangulator.newFacesForFace();
 118
 119     //- update face subsets
 120     mesh_.updateFaceSubsets(newFacesFromFace);
 121 }
 122
 123 void triangulateNonPlanarBaseFaces::decomposeCellsIntoPyramids()
 124 {
 125     decomposeCells sc(mesh_);
 126     sc.decomposeMesh(invertedCell_);
 127 }
 128
 129 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 130
 131 } // End namespace Foam
 132
 133 // ************************************************************************* //