applications/utilities/surface/surfaceRefineRedGreen/surfaceRefineRedGreen.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2009 OpenCFD Ltd.
   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 the
  13     Free Software Foundation; either version 2 of the License, or (at your
  14     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, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 Application
  26     surfaceRefineRedGreen
  27
  28 Description
  29     Refine by splitting all three edges of triangle ('red' refinement).
  30     Neighbouring triangles (which are not marked for refinement get split
  31     in half ('green') refinement. (R. Verfuerth, "A review of a posteriori
  32     error estimation and adaptive mesh refinement techniques",
  33     Wiley-Teubner, 1996)
  34
  35 \*---------------------------------------------------------------------------*/
  36
  37 #include "triSurface.H"
  38 #include "triSurfaceTools.H"
  39 #include "argList.H"
  40 #include "OFstream.H"
  41
  42 using namespace Foam;
  43
  44
  45 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  46 // Main program:
  47
  48 int main(int argc, char *argv[])
  49 {
  50     argList::noParallel();
  51     argList::validArgs.clear();
  52     argList::validArgs.append("surface file");
  53     argList::validArgs.append("output surface file");
  54     argList::argList args(argc, argv);
  55
  56     fileName surfFileName(args.additionalArgs()[0]);
  57     fileName outFileName(args.additionalArgs()[1]);
  58
  59     Info<< "Reading surface from " << surfFileName << " ..." << endl;
  60
  61     triSurface surf1(surfFileName);
  62
  63     // Refine
  64     triSurface surf2 = triSurfaceTools::redGreenRefine
  65     (
  66         surf1,
  67         identity(surf1.size())  //Hack: refine all
  68     );
  69
  70     Info<< "Original surface:" << endl
  71         << "    triangles     :" << surf1.size() << endl
  72         << "    vertices(used):" << surf1.nPoints() << endl
  73         << "Refined surface:" << endl
  74         << "    triangles     :" << surf2.size() << endl
  75         << "    vertices(used):" << surf2.nPoints() << endl << endl;
  76
  77
  78     Info<< "Writing refined surface to " << outFileName << " ..." << endl;
  79
  80     surf2.write(outFileName);
  81
  82     Info << "End\n" << endl;
  83
  84     return 0;
  85 }
  86
  87
  88 // ************************************************************************* //