ENH: probes: filter out points that cannot be found. Fixes #492.

This commit is contained in:
mattijs
2017-06-08 16:55:52 +01:00
parent 6186f91fea
commit ca7dbabdcf

View File

@ -129,11 +129,11 @@ void Foam::patchProbes::findElements(const fvMesh& mesh)
if (isA<emptyPolyPatch>(bm[patchi])) if (isA<emptyPolyPatch>(bm[patchi]))
{ {
WarningInFunction WarningInFunction
<< " The sample point: " << sample << " The sample point: " << sample
<< " belongs to " << patchi << " belongs to " << patchi
<< " which is an empty patch. This is not permitted. " << " which is an empty patch. This is not permitted. "
<< " This sample will not be included " << " This sample will not be included "
<< endl; << endl;
} }
else if (info.hit()) else if (info.hit())
{ {
@ -197,9 +197,12 @@ void Foam::patchProbes::findElements(const fvMesh& mesh)
elementList_ = -1; elementList_ = -1;
faceList_.setSize(nearest.size()); faceList_.setSize(nearest.size());
faceList_ = -1; faceList_ = -1;
processor_.setSize(nearest.size());
processor_ = -1;
forAll(nearest, sampleI) forAll(nearest, sampleI)
{ {
processor_[sampleI] = nearest[sampleI].second().second();
if (nearest[sampleI].second().second() == Pstream::myProcNo()) if (nearest[sampleI].second().second() == Pstream::myProcNo())
{ {
// Store the face to sample // Store the face to sample