zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

ENH: tetDecomposer: continue in case of neg-tet volume

Browse Source
This commit is contained in:
mattijs
2015-10-14 14:45:01 +01:00
committed by Andrew Heather
parent f4aa81fc84
commit 9e9f690de2
1 changed files with 4 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
src/dynamicMesh/polyTopoChange/polyTopoChange/tetDecomposer.C
View File

@ -2,8 +2,8 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | \\ / O peration |
\\ / A nd | Copyright (C) 2012-2015 OpenFOAM Foundation \\ / A nd | Copyright (C) 2012-2014 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation | Copyright (C) 2015 OpenCFD Ltd.
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
This file is part of OpenFOAM. This file is part of OpenFOAM.
@ -149,7 +149,7 @@ Foam::label Foam::tetDecomposer::triIndex(const label faceI, const label fp)
const const
{ {
const face& f = mesh_.faces()[faceI]; const face& f = mesh_.faces()[faceI];
const label fp0 = mesh_.tetBasePtIs()[faceI]; const label fp0 = max(0, mesh_.tetBasePtIs()[faceI]);
// Work out triangle index on this face // Work out triangle index on this face
label thisTriI; label thisTriI;
@ -436,7 +436,7 @@ void Foam::tetDecomposer::setRefinement
} }
else else
{ {
label fp0 = mesh_.tetBasePtIs()[faceI]; label fp0 = max(mesh_.tetBasePtIs()[faceI], 0);
label fp = f.fcIndex(fp0); label fp = f.fcIndex(fp0);
for (label triI = 0; triI < f.size()-2; triI++) for (label triI = 0; triI < f.size()-2; triI++)