zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

patchToPatches::intersection: Fix error part calculation for reversed targets

Browse Source
This commit is contained in:
Will Bainbridge
2023-12-21 19:26:08 +00:00
parent 431c58d3fe
commit 73a3fbe446
1 changed files with 4 additions and 9 deletions
Download Patch File Download Diff File Expand all files Collapse all files

13
src/meshTools/patchToPatch/intersection/intersectionPatchToPatch.C
View File

@ -394,14 +394,8 @@ bool Foam::patchToPatches::intersection::intersectFaces
if (tgtFaceEdgei < tgtPatch[tgtFacei].size()) if (tgtFaceEdgei < tgtPatch[tgtFacei].size())
{ {
if (reverse()) tgtFaceEdgePart_[tgtFaceEdgei] +=
{ reverse() ? ictEdgePart : -ictEdgePart;
tgtFaceEdgePart_[tgtFaceEdgei] += ictEdgePart;
}
else
{
tgtFaceEdgePart_[tgtFaceEdgei] -= ictEdgePart;
}
ictEdgeSide = -1; ictEdgeSide = -1;
} }
else else
@ -456,7 +450,8 @@ bool Foam::patchToPatches::intersection::intersectFaces
srcFaceEdgeParts_[srcFacei][srcFaceEdgei] += srcFaceEdgeParts_[srcFacei][srcFaceEdgei] +=
srcFaceEdgePart_[srcFaceEdgei]; srcFaceEdgePart_[srcFaceEdgei];
} }
srcErrorParts_[srcFacei] -= sum(tgtFaceEdgePart_); srcErrorParts_[srcFacei] +=
reverse() ? sum(tgtFaceEdgePart_) : -sum(tgtFaceEdgePart_);
srcErrorParts_[srcFacei] += errorPart; srcErrorParts_[srcFacei] += errorPart;
} }