From: eap <eap@opencascade.com>
Date: Mon, 17 Aug 2020 11:00:07 +0000 (+0300)
Subject: #19895 [CEA] Overlapping volume cells with viscous layer on shared faces
X-Git-Tag: V9_6_0a1~1
X-Git-Url: http://git.salome-platform.org/gitweb/?a=commitdiff_plain;h=bbf6d2faa3ad9ee967a75c65ff342bd88ab99233;p=plugins%2Fghs3dplugin.git

#19895 [CEA] Overlapping volume cells with viscous layer on shared faces
---

diff --git a/src/GHS3DPlugin/GHS3DPlugin_GHS3D.cxx b/src/GHS3DPlugin/GHS3DPlugin_GHS3D.cxx
index 64825c0..1f11791 100644
--- a/src/GHS3DPlugin/GHS3DPlugin_GHS3D.cxx
+++ b/src/GHS3DPlugin/GHS3DPlugin_GHS3D.cxx
@@ -263,9 +263,15 @@ static int checkTmpFace(const SMDS_MeshNode* node1,
         bool isDomainInPyramid = ( i3base != i3 );
         return isDomainInPyramid ? HOLE_ID : vol->getshapeId();
       }
-      else
+      else // nbNodes == 6
       {
-        return vol->getshapeId(); // triangle is a prism top
+        int i1 = vol->GetNodeIndex( node1 );
+        if (( i1 == 5 && i2 == 4 && i3 == 3 ) ||
+            ( i1 == 4 && i2 == 3 && i3 == 5 ) ||
+            ( i1 == 3 && i2 == 5 && i3 == 4 ))
+          return HOLE_ID;
+        else
+          return vol->getshapeId(); // triangle is a prism top
       }
     }
   }