exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Added check for faces with insufficient points

Browse source
This commit is contained in:
Hrvoje Jasak 2013-12-16 14:29:51 +00:00
parent 52f4078d64
commit ba3d7f689b
1 changed files with 8 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

10
src/foam/meshes/meshShapes/face/face.C
View file

@ -479,6 +479,12 @@ Foam::point Foam::face::centre(const pointField& meshPoints) const
{
// Calculate the centre by breaking the face into triangles and
// area-weighted averaging their centres
if (size() < 3)
{
FatalErrorIn("point face::centre(const pointField& meshPoints) const")
<< "Face with fewer than 3 points detected: " << *this
<< abort(FatalError);
}
// If the face is a triangle, do a direct calculation
if (size() == 3)
@ -495,7 +501,7 @@ Foam::point Foam::face::centre(const pointField& meshPoints) const
label nPoints = size();
point centrePoint = point::zero;
for (register label pI=0; pI<nPoints; pI++)
for (register label pI = 0; pI < nPoints; pI++)
{
centrePoint += meshPoints[operator[](pI)];
}
@ -504,7 +510,7 @@ Foam::point Foam::face::centre(const pointField& meshPoints) const
scalar sumA = 0;
vector sumAc = vector::zero;
for (register label pI=0; pI<nPoints; pI++)
for (register label pI = 0; pI < nPoints; pI++)
{
const point& nextPoint = meshPoints[operator[]((pI + 1) % nPoints)];