From cbd129bf53503b9ab02538f6ba90f07fd873791f Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Joachim=20Sch=C3=B6berl?= <joachim.schoeberl@tuwien.ac.at>
Date: Mon, 9 Sep 2019 10:55:27 +0200
Subject: [PATCH] disable parallel vertex2element table, needs sorting

---
 libsrc/meshing/topology.cpp | 13 ++++++++-----
 1 file changed, 8 insertions(+), 5 deletions(-)

diff --git a/libsrc/meshing/topology.cpp b/libsrc/meshing/topology.cpp
index 8fc40937..6fd54b8b 100644
--- a/libsrc/meshing/topology.cpp
+++ b/libsrc/meshing/topology.cpp
@@ -394,14 +394,14 @@ namespace netgen
        });
     
     vert2element = TABLE<ElementIndex,PointIndex::BASE> (cnt);
-    /*
     for (ElementIndex ei = 0; ei < ne; ei++)
       {
 	const Element & el = (*mesh)[ei];
 	for (int j = 0; j < el.GetNV(); j++)
 	  vert2element.AddSave (el[j], ei);
       }
-    */
+    
+    /*
     ParallelForRange
       (tm, ne,
        [&] (size_t begin, size_t end)
@@ -413,6 +413,8 @@ namespace netgen
                vert2element.ParallelAdd (el[j], ei);
            }
        });
+       requires sorting !!!!
+    */
     
     cnt = 0;
     /*
@@ -438,14 +440,14 @@ namespace netgen
     
 
     vert2surfelement = TABLE<SurfaceElementIndex,PointIndex::BASE> (cnt);
-    /*
+
     for (SurfaceElementIndex sei = 0; sei < nse; sei++)
       {
 	const Element2d & el = (*mesh)[sei];
 	for (int j = 0; j < el.GetNV(); j++)
 	  vert2surfelement.AddSave (el[j], sei);
       }
-    */
+    /*
     ParallelForRange
       (tm, nse,
        [&] (size_t begin, size_t end)
@@ -457,7 +459,8 @@ namespace netgen
                vert2surfelement.ParallelAdd (el[j], sei);
            }
        });
-
+       requires sorting !!!
+    */
 
     cnt = 0;
     for (SegmentIndex si = 0; si < nseg; si++)