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

ENH: cyclicACMI: report number of blended faces. Fixes #411.

Browse Source
This commit is contained in:
mattijs
2017-02-27 11:49:19 +00:00
parent a843d84f80
commit b71c4676ce
1 changed files with 54 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

55
src/meshTools/AMIInterpolation/patches/cyclicACMI/cyclicACMIPolyPatch/cyclicACMIPolyPatch.C
View File

@ -3,7 +3,7 @@
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | \\ / O peration |
\\ / A nd | Copyright (C) 2013-2016 OpenFOAM Foundation \\ / A nd | Copyright (C) 2013-2016 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation | Copyright (C) 2017 OpenCFD Ltd.
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
This file is part of OpenFOAM. This file is part of OpenFOAM.
@ -97,6 +97,59 @@ void Foam::cyclicACMIPolyPatch::resetAMI
AMIPatchToPatchInterpolation& AMI = AMIPatchToPatchInterpolation& AMI =
const_cast<AMIPatchToPatchInterpolation&>(this->AMI()); const_cast<AMIPatchToPatchInterpolation&>(this->AMI());
// Output some stats. AMIInterpolation will have already output the
// average weights ("sum(weights) min/max/average = 1, 1, 1")
{
const scalarField& wghtsSum = AMI.srcWeightsSum();
label nUncovered = 0;
label nCovered = 0;
forAll(wghtsSum, facei)
{
scalar sum = wghtsSum[facei];
if (sum < tolerance_)
{
nUncovered++;
}
else if (sum > scalar(1)-tolerance_)
{
nCovered++;
}
}
reduce(nUncovered, sumOp<label>());
reduce(nCovered, sumOp<label>());
label nTotal = returnReduce(wghtsSum.size(), sumOp<label>());
Info<< "ACMI: Patch source uncovered/blended/covered = "
<< nUncovered << ", " << nTotal-nUncovered-nCovered
<< ", " << nCovered << endl;
}
{
const scalarField& wghtsSum = AMI.tgtWeightsSum();
label nUncovered = 0;
label nCovered = 0;
forAll(wghtsSum, facei)
{
scalar sum = wghtsSum[facei];
if (sum < tolerance_)
{
nUncovered++;
}
else if (sum > scalar(1)-tolerance_)
{
nCovered++;
}
}
reduce(nUncovered, sumOp<label>());
reduce(nCovered, sumOp<label>());
label nTotal = returnReduce(wghtsSum.size(), sumOp<label>());
Info<< "ACMI: Patch target uncovered/blended/covered = "
<< nUncovered << ", " << nTotal-nUncovered-nCovered
<< ", " << nCovered << endl;
}
srcMask_ = srcMask_ =
min(scalar(1) - tolerance_, max(tolerance_, AMI.srcWeightsSum())); min(scalar(1) - tolerance_, max(tolerance_, AMI.srcWeightsSum()));