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

STYLE: reuse faceCentres if they already exist

Browse Source
This commit is contained in:
Mark Olesen 2020-04-28 15:24:56 +02:00
parent cc7af66a15
commit 8cbf55acad
1 changed files with 18 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

25
src/meshTools/searchableSurfaces/triSurfaceMesh/triSurfaceMesh.C
View File

@ -540,15 +540,26 @@ void Foam::triSurfaceMesh::clearOut()
Foam::tmp<Foam::pointField> Foam::triSurfaceMesh::coordinates() const
{
auto tpts = tmp<pointField>::New(8);
auto tpts = tmp<pointField>::New();
auto& pts = tpts.ref();
// Use copy to calculate face centres so they don't get stored
pts = PrimitivePatch<triSurface::FaceType, SubList, const pointField&>
(
SubList<triSurface::FaceType>(*this, triSurface::size()),
triSurface::points()
).faceCentres();
if (triSurface::hasFaceCentres())
{
// Can reuse existing values
pts = triSurface::faceCentres();
}
else
{
typedef SubList<labelledTri> FaceListType;
// Calculate face centres from a copy to avoid incurring
// additional storage
pts = PrimitivePatch<labelledTri, SubList, const pointField&>
(
FaceListType(*this, triSurface::size()),
triSurface::points()
).faceCentres();
}
return tpts;
}