ENH: tetOverlapVolume - no overelap when volumes are zero

This commit is contained in:
andy
2013-06-04 12:34:50 +01:00
parent 690b7e4ce4
commit 2323876fb0

View File

@ -64,6 +64,11 @@ Foam::scalar Foam::tetOverlapVolume::tetTetOverlapVol
tetPointRef::sumVolOp volInside;
tetPointRef::dummyOp outside;
if ((tetA.tet().mag() < SMALL) || (tetB.tet().mag() < SMALL))
{
return 0.0;
}
// face0
plane pl0(tetB[1], tetB[3], tetB[2]);
tetA.tet().sliceWithPlane(pl0, cutInside, outside);