src/mesh/cfMesh/utilities/smoothers/geometry/meshSurfaceOptimizer/meshSurfaceOptimizerCalculateTrianglesAndAddressing.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 "demandDrivenData.H"
  29 #include "meshSurfaceOptimizer.H"
  30 #include "meshSurfaceEngine.H"
  31
  32 //#define DEBUGTriangulation
  33
  34 # ifdef DEBUGTriangulation
  35 #include "triSurf.H"
  36 #include "triSurfModifier.H"
  37 #include "helperFunctions.H"
  38 # endif
  39
  40 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  41
  42 namespace Foam
  43 {
  44
  45 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  46
  47 void meshSurfaceOptimizer::calculateTrianglesAndAddressing() const
  48 {
  49     if( triMeshPtr_ )
  50         FatalErrorIn
  51         (
  52             "void meshSurfaceOptimizer::calculateTrianglesAndAddressing() const"
  53         ) << "Addressing is already calculated!" << abort(FatalError);
  54
  55     triMeshPtr_ = new partTriMesh(*partitionerPtr_);
  56
  57     # ifdef DEBUGTriangulation
  58     const labelList& sPoints = triMeshPtr_->meshSurfacePointLabelInTriMesh();
  59
  60     forAll(sPoints, bpI)
  61     {
  62         triSurf surf;
  63
  64         const label spI = sPoints[bpI];
  65
  66         partTriMeshSimplex simplex(*triMeshPtr_, spI);
  67
  68         triSurfModifier sMod(surf);
  69
  70         pointField& pts = sMod.pointsAccess();
  71         pts.setSize(simplex.pts().size());
  72         forAll(pts, i)
  73             pts[i] = simplex.pts()[i];
  74
  75         LongList<labelledTri>& trias = sMod.facetsAccess();
  76         trias.setSize(simplex.triangles().size());
  77         forAll(trias, i)
  78         {
  79             const triFace& t = simplex.triangles()[i];
  80             trias[i] = labelledTri(t[0], t[1], t[2], 0);
  81         }
  82
  83         sMod.patchesAccess().setSize(1);
  84         sMod.patchesAccess()[0].name() = "bnd";
  85         sMod.patchesAccess()[0].geometricType() = "patch";
  86
  87         fileName sName("bndPointSimplex_");
  88         sName += help::scalarToText(bpI);
  89         sName += ".stl";
  90         surf.writeSurface(sName);
  91     }
  92     # endif
  93 }
  94
  95 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  96
  97 } // End namespace Foam
  98
  99 // ************************************************************************* //