patch for embedded edges in OCC surfaces

Geo/GModelIO_OCC.cpp

@@ -108,7 +108,7 @@ void OCC_Internals::addShapeToLists(TopoDS_Shape _shape)
             if(fmap.FindIndex(face) < 1){
               fmap.Add(face);
 
-              for(exp3.Init(exp2.Current(), TopAbs_WIRE); exp3.More(); exp3.Next()){
+              for(exp3.Init(exp2.Current().Oriented(TopAbs_FORWARD), TopAbs_WIRE); exp3.More(); exp3.Next()){
                 TopoDS_Wire wire = TopoDS::Wire(exp3.Current());
                 if(wmap.FindIndex(wire) < 1){
                   wmap.Add(wire);

Geo/OCCFace.cpp

@@ -42,7 +42,7 @@
 #endif
 
 OCCFace::OCCFace(GModel *m, TopoDS_Face _s, int num)
-  : GFace(m, num), s(_s)
+  : GFace(m, num), s(TopoDS::Face(_s.Oriented(TopAbs_FORWARD)))
 {
   setup();
   model()->getOCCInternals()->bind(s, num);

Geo/OCCRegion.cpp

@@ -36,14 +36,20 @@ void OCCRegion::setup()
     for(exp3.Init(shell, TopAbs_FACE); exp3.More(); exp3.Next()){
       TopoDS_Face face = TopoDS::Face(exp3.Current());
       GFace *f = model()->getOCCInternals()->getOCCFaceByNativePtr(model(),face);
-      if(f){
+      if(!f){
+        Msg::Error("Unknown face in region %d", tag());
+      }
+      else if (face.Orientation() == TopAbs_INTERNAL){
+        Msg::Info("Adding embedded face %d", f->tag());
+        embedded_faces.push_back(f);
+      }
+      else{
         l_faces.push_back(f);
         f->addRegion(this);
       }
-      else
-        Msg::Error("Unknown face in region %d", tag());
     }
   }
+
   Msg::Debug("OCC Region %d with %d faces", tag(), l_faces.size());
 }