From 34590f1b9a54145e256397eb3a8b1ada76dba090 Mon Sep 17 00:00:00 2001
From: Christopher Lackner <clackner@cerbsim.com>
Date: Tue, 26 May 2020 11:23:33 +0200
Subject: [PATCH] support already meshed edges/faces in occ mesher

---
 libsrc/occ/occgenmesh.cpp | 6 ++++--
 1 file changed, 4 insertions(+), 2 deletions(-)

diff --git a/libsrc/occ/occgenmesh.cpp b/libsrc/occ/occgenmesh.cpp
index 38902710..314d405a 100644
--- a/libsrc/occ/occgenmesh.cpp
+++ b/libsrc/occ/occgenmesh.cpp
@@ -365,6 +365,7 @@ namespace netgen
           {
             TopoDS_Face face = TopoDS::Face(exp1.Current());
             int facenr = geom.fmap.FindIndex(face);
+            if(facenr < 1) continue;
 
             if (face2solid[0][facenr-1] == 0)
               face2solid[0][facenr-1] = solidnr;
@@ -382,8 +383,7 @@ namespace netgen
 
 
     int facenr = 0;
-    int edgenr = 0;
-
+    int edgenr = mesh.GetNSeg();
 
     (*testout) << "faces = " << geom.fmap.Extent() << endl;
     int curr = 0;
@@ -449,6 +449,8 @@ namespace netgen
                     continue;
                   }
 
+                if(geom.emap.FindIndex(edge) < 1) continue;
+
                 if (geom.vmap.FindIndex(TopExp::FirstVertex (edge)) ==
                     geom.vmap.FindIndex(TopExp::LastVertex (edge)))
                   {