applications/utilities/surface/surfaceAutoPatch/surfaceAutoPatch.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
   6      \\/     M anipulation  |
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of OpenFOAM.
  10
  11     OpenFOAM is free software: you can redistribute it and/or modify it
  12     under the terms of the GNU General Public License as published by
  13     the Free Software Foundation, either version 3 of the License, or
  14     (at your option) any later version.
  15
  16     OpenFOAM 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 OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 Application
  25     surfaceAutoPatch
  26
  27 Description
  28     Patches surface according to feature angle. Like autoPatch.
  29
  30 \*---------------------------------------------------------------------------*/
  31
  32 #include "triangle.H"
  33 #include "triSurface.H"
  34 #include "argList.H"
  35 #include "surfaceFeatures.H"
  36 #include "treeBoundBox.H"
  37 #include "meshTools.H"
  38 #include "OFstream.H"
  39
  40 using namespace Foam;
  41
  42 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  43
  44 // Main program:
  45
  46 int main(int argc, char *argv[])
  47 {
  48     argList::noParallel();
  49     argList::validArgs.append("input surfaceFile");
  50     argList::validArgs.append("output surfaceFile");
  51     argList::validArgs.append("includedAngle [0..180]");
  52     argList args(argc, argv);
  53
  54     const fileName inFileName  = args[1];
  55     const fileName outFileName = args[2];
  56     const scalar includedAngle = args.argRead<scalar>(3);
  57
  58     Info<< "Surface        : " << inFileName << nl << endl;
  59
  60
  61     // Read
  62     // ~~~~
  63
  64     Info<< "Reading : " << inFileName << endl;
  65     triSurface surf(inFileName);
  66
  67     Info<< "Read surface:" << endl;
  68     surf.writeStats(Info);
  69     Info<< endl;
  70
  71
  72
  73     // Construct features from surface&featureangle
  74     // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  75
  76     Info<< "Constructing feature set from included angle " << includedAngle
  77         << endl;
  78
  79     surfaceFeatures set(surf, includedAngle);
  80
  81     Info<< nl
  82         << "Feature set:" << nl
  83         << "    feature points : " << set.featurePoints().size() << nl
  84         << "    feature edges  : " << set.featureEdges().size() << nl
  85         << "    of which" << nl
  86         << "        region edges   : " << set.nRegionEdges() << nl
  87         << "        external edges : " << set.nExternalEdges() << nl
  88         << "        internal edges : " << set.nInternalEdges() << nl
  89         << endl;
  90
  91     // Get per-edge status.
  92     boolList borderEdge(surf.nEdges(), false);
  93     forAll(set.featureEdges(), i)
  94     {
  95         borderEdge[set.featureEdges()[i]] = true;
  96     }
  97
  98     labelList faceRegion(surf.size());
  99     label nRegions = surf.markZones(borderEdge, faceRegion);
 100
 101     // Reregion triangles.
 102     forAll(surf, i)
 103     {
 104         surf[i].region() = faceRegion[i];
 105     }
 106
 107     // Create some patches
 108     surf.patches().setSize(nRegions);
 109
 110     forAll(surf.patches(), patchI)
 111     {
 112         surf.patches()[patchI].name() = "patch" + Foam::name(patchI);
 113         surf.patches()[patchI].geometricType() = "empty";
 114     }
 115
 116
 117     Info<< "Writing : " << outFileName << endl;
 118     surf.write(outFileName, true);
 119
 120     Info<< "End\n" << endl;
 121
 122     return 0;
 123 }
 124
 125
 126 // ************************************************************************* //