songtianlun/openfoam
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

BUG: AMIInterpolation: reset cached data. Fixes #3334

Browse Source
This commit is contained in:
mattijs 2025-03-17 13:59:44 +00:00
parent 5cc36dc5b7
commit a77aaa7582
1 changed files with 40 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -6,7 +6,7 @@
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2011-2017 OpenFOAM Foundation
Copyright (C) 2015-2024 OpenCFD Ltd.
Copyright (C) 2015-2025 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@ -923,6 +923,36 @@ bool Foam::AMIInterpolation::calculate
addProfiling(ami, "AMIInterpolation::calculate");
// Clear storage (only needed if src/tgt become zero size)
{
if (srcMagSf_.size())
{
srcMagSf_.resize_nocopy(srcPatch.size());
}
srcAddress_.resize_nocopy(srcPatch.size());
srcWeights_.resize_nocopy(srcPatch.size());
srcWeightsSum_.resize_nocopy(srcPatch.size());
if (srcCentroids_.size())
{
srcCentroids_.resize_nocopy(srcPatch.size());
}
if (tgtPatch.size())
{
tgtMagSf_.resize_nocopy(tgtPatch.size());
}
tgtAddress_.resize_nocopy(tgtPatch.size());
tgtWeights_.resize_nocopy(tgtPatch.size());
tgtWeightsSum_.resize_nocopy(tgtPatch.size());
if (tgtCentroids_.size())
{
tgtCentroids_.resize_nocopy(tgtPatch.size());
}
}
if (surfPtr)
{
srcPatchPts_ = srcPatch.points();
@ -948,20 +978,13 @@ bool Foam::AMIInterpolation::calculate
}
// Note: use original communicator for statistics
label srcTotalSize = returnReduce
const label srcTotalSize = returnReduce
(
srcPatch.size(),
sumOp<label>(),
UPstream::msgType(),
comm_
);
label tgtTotalSize = returnReduce
(
tgtPatch.size(),
sumOp<label>(),
UPstream::msgType(),
comm_
);
if (srcTotalSize == 0)
{
@ -974,6 +997,14 @@ bool Foam::AMIInterpolation::calculate
return false;
}
const label tgtTotalSize = returnReduce
(
tgtPatch.size(),
sumOp<label>(),
UPstream::msgType(),
comm_
);
// Calculate:
// - which processors have faces
// - allocates a communicator (geomComm_) for those