zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

BUG: AMI - corrected appending of neighbour faces routine

Browse Source
This commit is contained in:
andy
2013-08-08 12:21:15 +01:00
parent 307d511a98
commit e54f156e60
1 changed files with 7 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

15
src/meshTools/AMIInterpolation/AMIInterpolation/AMIMethod/AMIMethod/AMIMethod.C
View File

@ -273,9 +273,8 @@ void Foam::AMIMethod<SourcePatch, TargetPatch>::appendNbrFaces
) const ) const
{ {
const labelList& nbrFaces = patch.faceFaces()[faceI]; const labelList& nbrFaces = patch.faceFaces()[faceI];
const pointField& points = patch.points();
// filter out faces already visited from src face neighbours // filter out faces already visited from face neighbours
forAll(nbrFaces, i) forAll(nbrFaces, i)
{ {
label nbrFaceI = nbrFaces[i]; label nbrFaceI = nbrFaces[i];
@ -301,15 +300,15 @@ void Foam::AMIMethod<SourcePatch, TargetPatch>::appendNbrFaces
} }
} }
// prevent addition of face if it is not on the same plane-ish
if (valid) if (valid)
{ {
const face& f = patch[faceI]; const vector& n1 = patch.faceNormals()[faceI];
const face& nbrF = patch[nbrFaceI]; const vector& n2 = patch.faceNormals()[nbrFaceI];
const vector& n = f.normal(points);
const vector& nbrN = nbrF.normal(points);
scalar cosI = mag(nbrN & n);
if (cosI < Foam::cos(degToRad(89.0))) scalar cosI = n1 & n2;
if (cosI > Foam::cos(degToRad(89.0)))
{ {
faceIDs.append(nbrFaceI); faceIDs.append(nbrFaceI);
} }