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

STYLE: code clean-up

Browse Source
This commit is contained in:
andy
2013-05-08 11:36:38 +01:00
parent c51d7eb60b
commit c87f83dc46
2 changed files with 7 additions and 9 deletions
Download Patch File Download Diff File Expand all files Collapse all files

7
src/meshTools/AMIInterpolation/AMIInterpolation/AMIInterpolation.C
View File

@ -190,12 +190,13 @@ void Foam::AMIInterpolation<SourcePatch, TargetPatch>::normaliseWeights
wghtSum.setSize(wght.size()); wghtSum.setSize(wght.size());
forAll(wght, faceI) forAll(wght, faceI)
{ {
scalar s = sum(wght[faceI]); scalarList& w = wght[faceI];
scalar s = sum(w);
scalar t = s/patchAreas[faceI]; scalar t = s/patchAreas[faceI];
forAll(addr[faceI], i) forAll(w, i)
{ {
wght[faceI][i] /= s; w[i] /= s;
} }
wghtSum[faceI] = t; wghtSum[faceI] = t;

9
src/meshTools/AMIInterpolation/AMIInterpolation/AMIMethod/faceAreaWeightAMI/faceAreaWeightAMI.C
View File

@ -67,7 +67,7 @@ bool Foam::faceAreaWeightAMI<SourcePatch, TargetPatch>::processSourceFace
visitedFaces.append(tgtFaceI); visitedFaces.append(tgtFaceI);
scalar area = interArea(srcFaceI, tgtFaceI); scalar area = interArea(srcFaceI, tgtFaceI);
// store when intersection area > 0 // store when intersection fractional area > tolerance
if (area/this->srcMagSf_[srcFaceI] > faceAreaIntersect::tolerance()) if (area/this->srcMagSf_[srcFaceI] > faceAreaIntersect::tolerance())
{ {
srcAddr[srcFaceI].append(tgtFaceI); srcAddr[srcFaceI].append(tgtFaceI);
@ -227,12 +227,9 @@ Foam::scalar Foam::faceAreaWeightAMI<SourcePatch, TargetPatch>::interArea
const face& tgt = this->tgtPatch_[tgtFaceI]; const face& tgt = this->tgtPatch_[tgtFaceI];
// quick reject if either face has zero area // quick reject if either face has zero area
// Note: do not used stored face areas for target patch // Note: do not use stored face areas for target patch
const scalar tgtMag = tgt.mag(tgtPoints); const scalar tgtMag = tgt.mag(tgtPoints);
if if ((this->srcMagSf_[srcFaceI] < ROOTVSMALL) || (tgtMag < ROOTVSMALL))
(
(this->srcMagSf_[srcFaceI] < ROOTVSMALL) || (tgtMag < ROOTVSMALL)
)
{ {
return 0.0; return 0.0;
} }