src/mesh/cfMesh/meshLibrary/utilities/triSurfaceTools/triSurfaceExtrude2DEdges/triSurfaceExtrude2DEdges.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 "triSurfaceExtrude2DEdges.H"
  29 #include "triSurfModifier.H"
  30 #include "boundBox.H"
  31
  32 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  33
  34 namespace Foam
  35 {
  36
  37 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  38
  39 triSurfaceExtrude2DEdges::triSurfaceExtrude2DEdges(const triSurf& surface)
  40 :
  41     surf_(surface)
  42 {}
  43
  44 triSurfaceExtrude2DEdges::~triSurfaceExtrude2DEdges()
  45 {}
  46
  47 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  48
  49 void triSurfaceExtrude2DEdges::extrudeSurface(triSurf& newSurf) const
  50 {
  51     triSurfModifier sMod(newSurf);
  52
  53     //- set patches
  54     geometricSurfacePatchList& patches = sMod.patchesAccess();
  55
  56     patches.setSize(1);
  57     patches[0].name() = "patch0";
  58     patches[0].geometricType() = "patch";
  59
  60     //- check if the edges are in the x-y plane
  61     const pointField& sPoints = surf_.points();
  62     const boundBox bb(sPoints);
  63
  64     if( Foam::mag(bb.max().z() - bb.min().z()) > SMALL )
  65         FatalErrorIn
  66         (
  67             "void triSurfaceExtrude2DEdges::extrudeSurface(triSurf&) const"
  68         ) << "Cannot extrude edges which are not in the x-y plane!"
  69           << exit(FatalError);
  70
  71     //- copy points
  72     pointField& pts = sMod.pointsAccess();
  73     pts.setSize(2 * sPoints.size());
  74
  75     const label nOffset = sPoints.size();
  76     const scalar zOffset = 0.1 * bb.mag();
  77
  78     forAll(sPoints, pI)
  79     {
  80         pts[pI] = pts[pI+nOffset] = sPoints[pI];
  81         pts[pI+sPoints.size()].z() += zOffset;
  82     }
  83
  84     //- create triangles from feature edges
  85     LongList<labelledTri>& triangles = sMod.facetsAccess();
  86     const edgeLongList& edges = surf_.featureEdges();
  87
  88     triangles.setSize(2 * edges.size());
  89     forAll(edges, eI)
  90     {
  91         const edge& e = edges[eI];
  92         const label tI = 2 * eI;
  93         triangles[tI] = labelledTri(e[0], e[1], e[1]+nOffset, 0);
  94         triangles[tI + 1] = labelledTri(e[0], e[1]+nOffset, e[0]+nOffset, 0);
  95     }
  96 }
  97
  98 const triSurf* triSurfaceExtrude2DEdges::extrudeSurface() const
  99 {
 100     triSurf* sPtr = new triSurf();
 101
 102     extrudeSurface(*sPtr);
 103
 104     return sPtr;
 105 }
 106
 107 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 108
 109 } // End namespace Foam
 110
 111 // ************************************************************************* //