From 9ef202872927b4398349c69118e730d0379fad69 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Joachim=20Sch=C3=B6berl?= <joachim.schoeberl@tuwien.ac.at>
Date: Sun, 28 Feb 2016 05:32:08 +0100
Subject: [PATCH] parallel update

---
 libsrc/meshing/refine.cpp | 11 ++++++++++-
 1 file changed, 10 insertions(+), 1 deletion(-)

diff --git a/libsrc/meshing/refine.cpp b/libsrc/meshing/refine.cpp
index 77a56428..49b908ea 100644
--- a/libsrc/meshing/refine.cpp
+++ b/libsrc/meshing/refine.cpp
@@ -17,7 +17,13 @@ namespace netgen
     mesh.ComputeNVertices();
     mesh.SetNP(mesh.GetNV());
 
+    if (mesh.mlbetweennodes.Size() < mesh.GetNV())
+      {
+        mesh.mlbetweennodes.SetSize(mesh.GetNV());
+        mesh.mlbetweennodes = INDEX_2(PointIndex::BASE-1,PointIndex::BASE-1);
+      }
 
+    
     INDEX_2_HASHTABLE<PointIndex> between(mesh.GetNP() + 5);
 
 
@@ -70,7 +76,10 @@ namespace netgen
     QuickSort (parents, par_nr);
 
     for (int i = 0; i < parents.Size(); i++)
-      between.Set (parents[i], mesh.GetNV()+i+PointIndex::BASE);
+      {
+        between.Set (parents[i], mesh.GetNV()+i+PointIndex::BASE);
+        mesh.mlbetweennodes[mesh.GetNV()+i+PointIndex::BASE] = parents[i];
+      }
     
     mesh.SetNP(mesh.GetNV() + parents.Size());
     Array<bool, PointIndex::BASE> pointset(mesh.GetNP());