applications/utilities/mesh/cfMesh/surfaceFeatureEdges/surfaceFeatureEdges.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     Finds feature edges and corners of a triangulated surface
  26
  27 \*---------------------------------------------------------------------------*/
  28
  29 #include "argList.H"
  30 #include "IFstream.H"
  31 #include "fileName.H"
  32 #include "triSurf.H"
  33 #include "OFstream.H"
  34 #include "OSspecific.H"
  35 #include "demandDrivenData.H"
  36 #include <cstdlib>
  37 #include <sstream>
  38
  39 #include "triSurfaceDetectFeatureEdges.H"
  40 #include "triSurfacePatchManipulator.H"
  41
  42 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  43 // Main program:
  44 using namespace Foam;
  45
  46 int main(int argc, char *argv[])
  47 {
  48     argList::noParallel();
  49     argList::validArgs.clear();
  50     argList::validArgs.append("input surface file");
  51     argList::validArgs.append("output surface file");
  52     argList::validOptions.insert("angle", "scalar");
  53     argList args(argc, argv);
  54
  55     fileName inFileName(args.args()[1]);
  56     fileName outFileName(args.args()[2]);
  57
  58     if (outFileName == inFileName)
  59     {
  60         FatalErrorIn(args.executable())
  61             << "Output file " << outFileName
  62             << " would overwrite the input file."
  63             << exit(FatalError);
  64     }
  65
  66     scalar tol(45.0);
  67     if( args.options().found("angle") )
  68     {
  69         const scalar ang = readScalar(IStringStream(args.options()["angle"])());
  70         tol = ang;
  71     }
  72     else
  73     {
  74         Info << "Using 45 deg as default angle!" << endl;
  75     }
  76
  77     triSurf originalSurface(inFileName);
  78
  79     triSurfaceDetectFeatureEdges edgeDetector(originalSurface, tol);
  80     edgeDetector.detectFeatureEdges();
  81
  82     if( outFileName.ext() == "fms" || outFileName.ext() == "FMS" )
  83     {
  84         Info << "Writing : " << outFileName << endl;
  85         originalSurface.writeSurface(outFileName);
  86     }
  87     else
  88     {
  89         triSurfacePatchManipulator manipulator(originalSurface);
  90         const triSurf* newSurfPtr = manipulator.surfaceWithPatches();
  91
  92         Info << "Writing : " << outFileName << endl;
  93         newSurfPtr->writeSurface(outFileName);
  94
  95         deleteDemandDrivenData(newSurfPtr);
  96     }
  97
  98     Info << "End\n" << endl;
  99
 100     return 0;
 101 }
 102
 103
 104 // ************************************************************************* //