ENH: PatchTools: test app
This commit is contained in:
mattijs
parent
257495f5b4
commit
5b0d85c3ff
3
applications/test/PatchTools/Make/files
Normal file
3
applications/test/PatchTools/Make/files
Normal file
@ -0,0 +1,3 @@
|
||||
Test-PatchTools.C
|
||||
|
||||
EXE = $(FOAM_USER_APPBIN)/Test-PatchTools
|
||||
7
applications/test/PatchTools/Make/options
Normal file
7
applications/test/PatchTools/Make/options
Normal file
@ -0,0 +1,7 @@
|
||||
EXE_INC = \
|
||||
-I$(LIB_SRC)/finiteVolume/lnInclude \
|
||||
-I$(LIB_SRC)/meshTools/lnInclude \
|
||||
-I$(LIB_SRC)/surfMesh/lnInclude
|
||||
|
||||
EXE_LIBS = \
|
||||
-lfiniteVolume
|
||||
297
applications/test/PatchTools/Test-PatchTools.C
Normal file
297
applications/test/PatchTools/Test-PatchTools.C
Normal file
@ -0,0 +1,297 @@
|
||||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2012-2013 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM 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.
|
||||
|
||||
OpenFOAM 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 OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
Application
|
||||
testPatchTools
|
||||
|
||||
Description
|
||||
Test app for PatchTools functionality
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#include "PatchTools.H"
|
||||
#include "argList.H"
|
||||
#include "fvMesh.H"
|
||||
#include "volFields.H"
|
||||
#include "Time.H"
|
||||
#include "OBJstream.H"
|
||||
|
||||
using namespace Foam;
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
//template<class PatchType>
|
||||
//Foam::tmp<Foam::pointField>
|
||||
//areaPointNormals
|
||||
//(
|
||||
// const polyMesh& mesh,
|
||||
// const PatchType& p,
|
||||
// const labelList& meshFaces
|
||||
//)
|
||||
//{
|
||||
// // Assume patch is smaller than the globalData().coupledPatch() (?) so
|
||||
// // loop over patch meshPoints.
|
||||
//
|
||||
// const labelList& meshPoints = p.meshPoints();
|
||||
//
|
||||
// const globalMeshData& globalData = mesh.globalData();
|
||||
// const indirectPrimitivePatch& coupledPatch = globalData.coupledPatch();
|
||||
// const Map<label>& coupledPatchMP = coupledPatch.meshPointMap();
|
||||
// const mapDistribute& map = globalData.globalPointSlavesMap();
|
||||
// const globalIndexAndTransform& transforms =
|
||||
// globalData.globalTransforms();
|
||||
//
|
||||
//
|
||||
// // 1. Start off with local (area-weighted) normals
|
||||
// // (note:without calculating pointNormals
|
||||
// // to avoid them being stored)
|
||||
//
|
||||
// tmp<pointField> textrudeN(new pointField(p.nPoints(), vector::zero));
|
||||
// pointField& extrudeN = textrudeN();
|
||||
// {
|
||||
// const faceList& localFaces = p.localFaces();
|
||||
// const vectorField& faceAreas = mesh.faceAreas();
|
||||
//
|
||||
// forAll(localFaces, faceI)
|
||||
// {
|
||||
// const face& f = localFaces[faceI];
|
||||
// const vector& n = faceAreas[meshFaces[faceI]];
|
||||
// forAll(f, fp)
|
||||
// {
|
||||
// extrudeN[f[fp]] += n;
|
||||
// }
|
||||
// }
|
||||
// }
|
||||
//
|
||||
//
|
||||
// // Collect local pointFaces
|
||||
// List<List<point> > pointFaceNormals(map.constructSize());
|
||||
// {
|
||||
// const vectorField& faceAreas = mesh.faceAreas();
|
||||
//
|
||||
// forAll(meshPoints, patchPointI)
|
||||
// {
|
||||
// label meshPointI = meshPoints[patchPointI];
|
||||
// Map<label>::const_iterator fnd = coupledPatchMP.find(meshPointI);
|
||||
// if (fnd != coupledPatchMP.end())
|
||||
// {
|
||||
// label coupledPointI = fnd();
|
||||
//
|
||||
// List<point>& pNormals = pointFaceNormals[coupledPointI];
|
||||
// const labelList& pFaces = p.pointFaces()[patchPointI];
|
||||
// pNormals.setSize(pFaces.size());
|
||||
// forAll(pFaces, i)
|
||||
// {
|
||||
// pNormals[i] = faceAreas[meshFaces[pFaces[i]]];
|
||||
// }
|
||||
// }
|
||||
// }
|
||||
// }
|
||||
//
|
||||
// // Pull remote data into local slots
|
||||
// map.distribute
|
||||
// (
|
||||
// transforms,
|
||||
// pointFaceNormals,
|
||||
// listTransform()
|
||||
// );
|
||||
//
|
||||
//
|
||||
// // Combine normals
|
||||
// const labelListList& slaves = globalData.globalPointSlaves();
|
||||
// const labelListList& transformedSlaves =
|
||||
// globalData.globalPointTransformedSlaves();
|
||||
//
|
||||
//
|
||||
// pointField coupledPointNormals(map.constructSize(), vector::zero);
|
||||
//
|
||||
// forAll(meshPoints, patchPointI)
|
||||
// {
|
||||
// label meshPointI = meshPoints[patchPointI];
|
||||
// Map<label>::const_iterator fnd = coupledPatchMP.find(meshPointI);
|
||||
// if (fnd != coupledPatchMP.end())
|
||||
// {
|
||||
// label coupledPointI = fnd();
|
||||
// const labelList& slaveSlots = slaves[coupledPointI];
|
||||
// const labelList& transformedSlaveSlots =
|
||||
// transformedSlaves[coupledPointI];
|
||||
//
|
||||
// label nFaces = slaveSlots.size()+transformedSlaveSlots.size();
|
||||
// if (nFaces > 0)
|
||||
// {
|
||||
// // Combine
|
||||
// point& n = coupledPointNormals[coupledPointI];
|
||||
//
|
||||
// n += sum(pointFaceNormals[coupledPointI]);
|
||||
//
|
||||
// forAll(slaveSlots, i)
|
||||
// {
|
||||
// n += sum(pointFaceNormals[slaveSlots[i]]);
|
||||
// }
|
||||
// forAll(transformedSlaveSlots, i)
|
||||
// {
|
||||
// n += sum(pointFaceNormals[transformedSlaveSlots[i]]);
|
||||
// }
|
||||
//
|
||||
// // Put back into slave slots
|
||||
// forAll(slaveSlots, i)
|
||||
// {
|
||||
// coupledPointNormals[slaveSlots[i]] = n;
|
||||
// }
|
||||
// forAll(transformedSlaveSlots, i)
|
||||
// {
|
||||
// coupledPointNormals[transformedSlaveSlots[i]] = n;
|
||||
// }
|
||||
// }
|
||||
// }
|
||||
// }
|
||||
//
|
||||
//
|
||||
// // Send back
|
||||
// map.reverseDistribute
|
||||
// (
|
||||
// transforms,
|
||||
// coupledPointNormals.size(),
|
||||
// coupledPointNormals,
|
||||
// mapDistribute::transform()
|
||||
// );
|
||||
//
|
||||
//
|
||||
// // Override patch normals
|
||||
// forAll(meshPoints, patchPointI)
|
||||
// {
|
||||
// label meshPointI = meshPoints[patchPointI];
|
||||
// Map<label>::const_iterator fnd = coupledPatchMP.find(meshPointI);
|
||||
// if (fnd != coupledPatchMP.end())
|
||||
// {
|
||||
// label coupledPointI = fnd();
|
||||
// extrudeN[patchPointI] = coupledPointNormals[coupledPointI];
|
||||
// }
|
||||
// }
|
||||
//
|
||||
// extrudeN /= mag(extrudeN)+VSMALL;
|
||||
//
|
||||
// return textrudeN;
|
||||
//}
|
||||
|
||||
|
||||
|
||||
// Main program:
|
||||
|
||||
int main(int argc, char *argv[])
|
||||
{
|
||||
# include "addTimeOptions.H"
|
||||
argList::validArgs.append("patch");
|
||||
# include "setRootCase.H"
|
||||
# include "createTime.H"
|
||||
|
||||
# include "createMesh.H"
|
||||
|
||||
const word patchName = args[1];
|
||||
label patchI = mesh.boundaryMesh().findPatchID(patchName);
|
||||
const polyPatch& pp = mesh.boundaryMesh()[patchI];
|
||||
|
||||
const indirectPrimitivePatch& cpp = mesh.globalData().coupledPatch();
|
||||
|
||||
{
|
||||
OBJstream str(runTime.path()/"edgePatchNormals.obj");
|
||||
|
||||
labelList patchEdges;
|
||||
labelList coupledEdges;
|
||||
PackedBoolList sameEdgeOrientation;
|
||||
PatchTools::matchEdges
|
||||
(
|
||||
pp,
|
||||
cpp,
|
||||
patchEdges,
|
||||
coupledEdges,
|
||||
sameEdgeOrientation
|
||||
);
|
||||
|
||||
const pointField en
|
||||
(
|
||||
PatchTools::edgeNormals
|
||||
(
|
||||
mesh,
|
||||
pp,
|
||||
patchEdges,
|
||||
coupledEdges
|
||||
)
|
||||
);
|
||||
|
||||
forAll(en, patchEdgeI)
|
||||
{
|
||||
const edge& patchE = pp.edges()[patchEdgeI];
|
||||
//str.write(pp.localPoints()[pointI], en[pointI]);
|
||||
const point pt = patchE.centre(pp.localPoints());
|
||||
str.write(linePointRef(pt, pt + 0.1*en[patchEdgeI]));
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
return 0;
|
||||
|
||||
|
||||
// {
|
||||
// OBJstream str(runTime.path()/"unweightedPatchNormals.obj");
|
||||
//
|
||||
// const pointField pn
|
||||
// (
|
||||
// PatchTools::pointNormals
|
||||
// (
|
||||
// mesh,
|
||||
// pp,
|
||||
// identity(pp.size())+pp.start()
|
||||
// )
|
||||
// );
|
||||
// forAll(pn, pointI)
|
||||
// {
|
||||
// str.write(linePointRef(pp.localPoints()[pointI], pn[pointI]));
|
||||
// }
|
||||
// }
|
||||
// {
|
||||
// OBJstream str(runTime.path()/"areaWeightedPatchNormals.obj");
|
||||
//
|
||||
// const pointField pn
|
||||
// (
|
||||
// areaPointNormals
|
||||
// (
|
||||
// mesh,
|
||||
// pp,
|
||||
// identity(pp.size())+pp.start()
|
||||
// )
|
||||
// );
|
||||
// forAll(pn, pointI)
|
||||
// {
|
||||
// str.write(linePointRef(pp.localPoints()[pointI], pn[pointI]));
|
||||
// }
|
||||
// }
|
||||
|
||||
|
||||
Pout<< "End\n" << endl;
|
||||
|
||||
return 0;
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
Loading…
Reference in New Issue
Block a user