diff --git a/applications/utilities/mesh/generation/makeFaMesh/Make/files b/applications/utilities/mesh/generation/makeFaMesh/Make/files deleted file mode 100644 index fbd0f40cd..000000000 --- a/applications/utilities/mesh/generation/makeFaMesh/Make/files +++ /dev/null @@ -1,3 +0,0 @@ -makeFaMesh.C - -EXE = $(FOAM_APPBIN)/makeFaMesh diff --git a/applications/utilities/mesh/generation/makeFaMesh/Make/options b/applications/utilities/mesh/generation/makeFaMesh/Make/options deleted file mode 100644 index 98ba4f722..000000000 --- a/applications/utilities/mesh/generation/makeFaMesh/Make/options +++ /dev/null @@ -1,8 +0,0 @@ -EXE_INC = \ - -I$(LIB_SRC)/finiteArea/lnInclude \ - -I$(LIB_SRC)/finiteVolume/lnInclude \ - -I$(LIB_SRC)/cfdTools/general/lnInclude - -EXE_LIBS = \ - -lfiniteArea \ - -lfiniteVolume diff --git a/applications/utilities/mesh/generation/makeFaMesh/makeFaMesh.C b/applications/utilities/mesh/generation/makeFaMesh/makeFaMesh.C deleted file mode 100644 index 010c647f2..000000000 --- a/applications/utilities/mesh/generation/makeFaMesh/makeFaMesh.C +++ /dev/null @@ -1,328 +0,0 @@ -/*---------------------------------------------------------------------------*\ - ========= | - \\ / F ield | foam-extend: Open Source CFD - \\ / O peration | - \\ / A nd | For copyright notice see file Copyright - \\/ M anipulation | -------------------------------------------------------------------------------- -License - This file is part of foam-extend. - - foam-extend 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. - - foam-extend 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 foam-extend. If not, see . - -Application - makeFaMesh - -Description - A mesh generator for finite area mesh. - -\*---------------------------------------------------------------------------*/ - -#include "objectRegistry.H" -#include "Time.H" -#include "argList.H" -#include "OSspecific.H" -#include "faMesh.H" -#include "fvMesh.H" - -using namespace Foam; - -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -class faPatchData -{ -public: - word name_; - word type_; - dictionary dict_; - label ownPolyPatchID_; - label ngbPolyPatchID_; - labelList edgeLabels_; - faPatchData() - : - name_(word::null), - type_(word::null), - ownPolyPatchID_(-1), - ngbPolyPatchID_(-1) - {} -}; - - -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -int main(int argc, char *argv[]) -{ -# include "setRootCase.H" -# include "createTime.H" -# include "createMesh.H" - - // Reading faMeshDefinition dictionary - IOdictionary faMeshDefinition - ( - IOobject - ( - "faMeshDefinition", - runTime.constant(), - "faMesh", - mesh, - IOobject::MUST_READ, - IOobject::NO_WRITE - ) - ); - - wordList polyMeshPatches - ( - faMeshDefinition.lookup("polyMeshPatches") - ); - - dictionary bndDict = faMeshDefinition.subDict("boundary"); - - wordList faPatchNames = bndDict.toc(); - - List faPatches(faPatchNames.size()+1); - - forAll (faPatchNames, patchI) - { - dictionary curPatchDict = - bndDict.subDict(faPatchNames[patchI]); - - faPatches[patchI].name_ = faPatchNames[patchI]; - - faPatches[patchI].type_ = - word(curPatchDict.lookup("type")); - - word ownName = curPatchDict.lookup("ownerPolyPatch"); - - faPatches[patchI].ownPolyPatchID_ = - mesh.boundaryMesh().findPatchID(ownName); - - if ( faPatches[patchI].ownPolyPatchID_ < 0 ) - { - FatalErrorIn("makeFaMesh:") - << "neighbourPolyPatch " << ownName << " does not exist" - << exit(FatalError); - } - - word neiName = curPatchDict.lookup("neighbourPolyPatch"); - - faPatches[patchI].ngbPolyPatchID_ = - mesh.boundaryMesh().findPatchID(neiName); - - if ( faPatches[patchI].ngbPolyPatchID_ < 0 ) - { - FatalErrorIn("makeFaMesh:") - << "neighbourPolyPatch " << neiName << " does not exist" - << exit(FatalError); - } - } - - // Setting faceLabels list size - label size = 0; - - labelList patchIDs(polyMeshPatches.size(), -1); - - forAll (polyMeshPatches, patchI) - { - patchIDs[patchI] = - mesh.boundaryMesh().findPatchID(polyMeshPatches[patchI]); - - if ( patchIDs[patchI] < 0 ) - { - FatalErrorIn("makeFaMesh:") - << "Patch " << polyMeshPatches[patchI] << " does not exist" - << exit(FatalError); - } - - size += mesh.boundaryMesh()[patchIDs[patchI]].size(); - } - - labelList faceLabels(size, -1); - - sort(patchIDs); - - - // Filling of faceLabels list - label faceI = -1; - - forAll (polyMeshPatches, patchI) - { - label start = mesh.boundaryMesh()[patchIDs[patchI]].start(); - - label size = mesh.boundaryMesh()[patchIDs[patchI]].size(); - - for(label i = 0; i < size; i++) - { - faceLabels[++faceI] = start + i; - } - } - - // Creating faMesh - Info << "Create faMesh ... "; - - faMesh areaMesh - ( - mesh, - faceLabels - ); - Info << "Done" << endl; - - - // Determination of faPatch ID for each boundary edge. - // Result is in the bndEdgeFaPatchIDs list - const indirectPrimitivePatch& patch = areaMesh.patch(); - - labelList faceCells(faceLabels.size(), -1); - - forAll (faceCells, faceI) - { - label faceID = faceLabels[faceI]; - - faceCells[faceI] = mesh.faceOwner()[faceID]; - } - - labelList meshEdges = - patch.meshEdges - ( - mesh.edges(), - mesh.cellEdges(), - faceCells - ); - - const labelListList& edgeFaces = mesh.edgeFaces(); - - const label nTotalEdges = patch.nEdges(); - const label nInternalEdges = patch.nInternalEdges(); - - labelList bndEdgeFaPatchIDs(nTotalEdges - nInternalEdges, -1); - - for (label edgeI = nInternalEdges; edgeI < nTotalEdges; edgeI++) - { - label curMeshEdge = meshEdges[edgeI]; - - labelList curEdgePatchIDs(2, -1); - - label patchI = -1; - - forAll (edgeFaces[curMeshEdge], faceI) - { - label curFace = edgeFaces[curMeshEdge][faceI]; - - label curPatchID = mesh.boundaryMesh().whichPatch(curFace); - - if (curPatchID != -1) - { - curEdgePatchIDs[++patchI] = curPatchID; - } - } - - for(label pI = 0; pI < faPatches.size() - 1; pI++) - { - if - ( - ( - curEdgePatchIDs[0] == faPatches[pI].ownPolyPatchID_ - && curEdgePatchIDs[1] == faPatches[pI].ngbPolyPatchID_ - ) - || - ( - curEdgePatchIDs[1] == faPatches[pI].ownPolyPatchID_ - && curEdgePatchIDs[0] == faPatches[pI].ngbPolyPatchID_ - ) - ) - { - bndEdgeFaPatchIDs[edgeI - nInternalEdges] = pI; - break; - } - } - } - - - // Set edgeLabels for each faPatch - for(label pI=0; pI<(faPatches.size()-1); pI++) - { - SLList