From 9c97eb52521f7fb3d50d67e5cf16ebbd1f9602fc Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Joachim=20Sch=C3=B6berl?= <joachim.schoeberl@tuwien.ac.at>
Date: Thu, 18 Aug 2016 18:01:17 +0200
Subject: [PATCH] fix for quad-faces

---
 libsrc/meshing/topology.cpp | 14 +++++++-------
 1 file changed, 7 insertions(+), 7 deletions(-)

diff --git a/libsrc/meshing/topology.cpp b/libsrc/meshing/topology.cpp
index 8e122fd3..4eac1c45 100644
--- a/libsrc/meshing/topology.cpp
+++ b/libsrc/meshing/topology.cpp
@@ -139,11 +139,12 @@ namespace netgen
                   swap (face4.I2(), face4.I4());
                 }
               
+              if (face4.I1() != v) continue;
+              
               func(face4, elnr, j, true, facedir);
                 /*
               INDEX_3 face(face4.I1(), face4.I2(), face4.I3());
               
-              if (face.I1() != v) continue;
               
               if (vert2face.Used (face))
                 {
@@ -256,10 +257,10 @@ namespace netgen
                 swap (face4.I2(), face4.I4());
               }
             
+            if (face4.I1() != v) continue;
             func(face4, elnr, 0, false, facedir);
             /*
               INDEX_3 face(face4.I1(), face4.I2(), face4.I3());
-              if (face.I1() != v) continue;
 		
               if (vert2face.Used (face))
               facenum = vert2face.Get(face);
@@ -638,7 +639,6 @@ namespace netgen
 
 
 
-
         NgProfiler::StopTimer (timer2a);
         NgProfiler::StartTimer (timer2b);
 
@@ -764,8 +764,9 @@ namespace netgen
                                  });
                 }
             });
-              /*
-                int oldnfa = face2vert.Size();
+        
+          /*
+          int oldnfa = face2vert.Size();
         int nfa = oldnfa;
         INDEX_3_CLOSED_HASHTABLE<int> vert2face(2*max_face_on_vertex+10); 
 
@@ -1010,9 +1011,8 @@ namespace netgen
 		  }
 	      }
 	  }
-*/
-
         face2vert.SetAllocSize (nfa);            
+          */
 
 	// *testout << "face2vert = " << endl << face2vert << endl;