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: PointEdgeWave: incorrect indexing in handleCollocatedPoints

Browse Source
This commit is contained in:
mattijs
2013-01-17 14:55:13 +00:00
parent 0f05652a98
commit 3b637a0d54
1 changed files with 26 additions and 17 deletions
Download Patch File Download Diff File Expand all files Collapse all files

43
src/meshTools/algorithms/PointEdgeWave/PointEdgeWave.C
View File

@ -566,7 +566,8 @@ Foam::label Foam::PointEdgeWave<Type, TrackingData>::handleCollocatedPoints()
elems[pointI] = allPointInfo_[meshPoints[pointI]];
}
// Pull slave data onto master. No need to update transformed slots.
// Pull slave data onto master (which might or might not have any
// initialised points). No need to update transformed slots.
slavesMap.distribute(elems, false);
// Combine master data with slave data
@ -597,26 +598,34 @@ Foam::label Foam::PointEdgeWave<Type, TrackingData>::handleCollocatedPoints()
// Extract back onto mesh
forAll(meshPoints, pointI)
{
Type& elem = allPointInfo_[meshPoints[pointI]];
// Like updatePoint but bypass Type::updatePoint with its tolerance
// checking
if (!elem.valid(td_) || !elem.equal(elems[pointI], td_))
if (elems[pointI].valid(td_))
{
nEvals_++;
elem = elems[pointI];
label meshPointI = meshPoints[pointI];
// See if element now valid
if (elem.valid(td_))
{
--nUnvisitedPoints_;
}
Type& elem = allPointInfo_[meshPointI];
// Update database of changed points
if (!changedPoint_[pointI])
bool wasValid = elem.valid(td_);
// Like updatePoint but bypass Type::updatePoint with its tolerance
// checking
//if (!elem.valid(td_) || !elem.equal(elems[pointI], td_))
if (!elem.equal(elems[pointI], td_))
{
changedPoint_[pointI] = true;
changedPoints_[nChangedPoints_++] = pointI;
nEvals_++;
elem = elems[pointI];
// See if element now valid
if (!wasValid && elem.valid(td_))
{
--nUnvisitedPoints_;
}
// Update database of changed points
if (!changedPoint_[meshPointI])
{
changedPoint_[meshPointI] = true;
changedPoints_[nChangedPoints_++] = meshPointI;
}
}
}
}