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: Problem with list combine gather for multiprocessor run

Browse Source
This commit is contained in:
sergio
2011-03-02 16:17:50 +00:00
parent 5b88f8e38f
commit 4e1d711858
1 changed files with 16 additions and 35 deletions
Download Patch File Download Diff File Expand all files Collapse all files

51
src/sampling/probes/patchProbes.C
View File

@ -2,7 +2,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) 2010-2010 OpenCFD Ltd. \\ / A nd | Copyright (C) 2010-2011 OpenCFD Ltd.
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -53,23 +53,15 @@ void Foam::patchProbes::findElements(const fvMesh& mesh)
{ {
const vector& sample = operator[](probeI); const vector& sample = operator[](probeI);
label faceI = meshSearchEngine.findNearestBoundaryFace(sample); label faceI = meshSearchEngine.findNearestBoundaryFace(sample);
if (faceI == -1) const point& fc = mesh.faceCentres()[faceI];
{ nearest[probeI].first() = pointIndexHit
nearest[probeI].second().first() = Foam::sqr(GREAT); (
nearest[probeI].second().second() = Pstream::myProcNo(); true,
} fc,
else faceI
{ );
const point& fc = mesh.faceCentres()[faceI]; nearest[probeI].second().first() = magSqr(fc-sample);
nearest[probeI].first() = pointIndexHit nearest[probeI].second().second() = Pstream::myProcNo();
(
true,
fc,
faceI
);
nearest[probeI].second().first() = magSqr(fc-sample);
nearest[probeI].second().second() = Pstream::myProcNo();
}
} }
@ -92,27 +84,16 @@ void Foam::patchProbes::findElements(const fvMesh& mesh)
} }
} }
// Check if all patchProbes have been found. // Check if all patchProbes have been found.
forAll(nearest, sampleI) forAll(nearest, sampleI)
{ {
label localI = nearest[sampleI].first().index(); label localI = -1;
if (nearest[sampleI].second().second() == Pstream::myProcNo())
{
localI = nearest[sampleI].first().index();
}
if (localI == -1) elementList_[sampleI] = localI;
{
if (Pstream::master())
{
WarningIn("patchProbes::findElements()")
<< "Did not find location "
<< nearest[sampleI].second().first()
<< " in any cell. Skipping location." << endl;
}
}
else
{
elementList_[sampleI] = localI;
}
} }
} }