FLMAToSurface.C 3.87 KB
Newer Older
Franjo's avatar
Franjo committed
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145
/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | cfMesh: A library for mesh generation
   \\    /   O peration     |
    \\  /    A nd           | Author: Franjo Juretic (franjo.juretic@c-fields.com)
     \\/     M anipulation  | Copyright (C) Creative Fields, Ltd.
-------------------------------------------------------------------------------
License
    This file is part of cfMesh.

    cfMesh is free software; you can redistribute it and/or modify it
    under the terms of the GNU General Public License as published by the
    Free Software Foundation; either version 3 of the License, or (at your
    option) any later version.

    cfMesh is distributed in the hope that it will be useful, but WITHOUT
    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
    for more details.

    You should have received a copy of the GNU General Public License
    along with cfMesh.  If not, see <http://www.gnu.org/licenses/>.

Description
    Reads the AVL's surface mesh

\*---------------------------------------------------------------------------*/

#include "argList.H"
#include "triSurf.H"
#include "triSurfModifier.H"
#include "triFaceList.H"
#include "labelLongList.H"
#include "IFstream.H"

using namespace Foam;

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

int main(int argc, char *argv[])
{
    argList::noParallel();
    argList::validArgs.clear();

    argList::validArgs.append("input surface file");
    argList::validArgs.append("output surface file");
    argList args(argc, argv);

    fileName inFileName(args.args()[1]);
    fileName outFileName(args.args()[2]);

    if( inFileName.ext() != "flma" )
    {
        Info << "Cannot convert this mesh" << endl;
        return 0;
    }

    //- create the surface mesh
    triSurf ts;
    triSurfModifier tsm(ts);

    label counter;

    IFstream inFile(inFileName);

    inFile >> counter;

    //- read vertices
    pointField& points = tsm.pointsAccess();
    points.setSize(counter);
    forAll(points, pointI)
    {
        point& p = points[pointI];

        inFile >> p.x();
        inFile >> p.y();
        inFile >> p.z();
    }

    //- read facets
    inFile >> counter;
    geometricSurfacePatchList patches(1);
    patches[0].name() = "patch";
    LongList<labelledTri>& triangles = tsm.facetsAccess();
    triangles.setSize(counter);
    forAll(triangles, triI)
    {
        inFile >> counter;

        if( counter != 3 )
        {
            Info << "Facet " << triI << " is not a triangle!!" << endl;
            Warning << "Cannot convert this surface!" << endl;
            return 0;
        }

        for(label j=0;j<3;++j)
            inFile >> triangles[triI][2-j];

        triangles[triI].region() = 0;
    }

    //- read cell types
    inFile >> counter;
    forAll(triangles, triI)
        inFile >> counter;

    //- start reading selections
    inFile >> counter;
    for(label selI=0;selI<counter;++selI)
    {
        //- read selection name
        word selName;
        inFile >> selName;

        //- read selection type
        label selType;
        inFile >> selType;

        //- read selection entries
        label size;
        inFile >> size;
        labelLongList entries(size);
        for(label i=0;i<size;++i)
            inFile >> entries[i];

        //- store cell selections
        if( selType == 2 )
        {
            Info << "Adding subset " << selName << endl;
            const label setID = ts.addFacetSubset(selName);

            forAll(entries, i)
                ts.addFacetToSubset(setID, entries[i]);
        }
    }

    //- write the surface
    ts.writeSurface(outFileName);

    Info << "End\n" << endl;
    return 0;
}

// ************************************************************************* //