Skip to content
Snippets Groups Projects
Commit dfb652bc authored by mattijs's avatar mattijs
Browse files

ENH: findRefCell: initialise to -1 if not needed. Fixes #988.

parent c84390ab
Branches
Tags
No related merge requests found
......@@ -3,7 +3,7 @@
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
\\/ M anipulation | Copyright (C) 2018 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
......@@ -111,6 +111,10 @@ bool Foam::setRefCell
return true;
}
else
{
refCelli = -1;
}
return false;
}
......@@ -141,25 +145,14 @@ Foam::scalar Foam::getRefCellValue
FatalErrorInFunction
<< "Illegal reference cellID " << refCelli
<< ". Mesh has " << field.mesh().nCells() << ". cells."
<< exit(FatalIOError);
}
// Catch duplicates
label hasRef =
(refCelli >= 0 && refCelli < field.mesh().nCells() ? 1 : 0);
if (1 != returnReduce<label>(hasRef, sumOp<label>()))
{
FatalErrorInFunction
<< "Invalid reference cellID " << refCelli
<< ". Mesh has " << field.mesh().nCells() << ". cells."
<< exit(FatalIOError);
<< exit(FatalError);
}
scalar refCellValue = (hasRef ? field[refCelli] : 0.0);
#else
scalar refCellValue = (refCelli >= 0 ? field[refCelli] : 0.0);
#endif
scalar refCellValue = (refCelli >= 0 ? field[refCelli] : 0.0);
// Currently distributing the value to all processors. This is generally
// not needed since only the processor holding the reference cell needs
// it. Tdb.
return returnReduce(refCellValue, sumOp<scalar>());
}
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment