From a3ccb6432c657bb94611eb26fe98b0fb72155e39 Mon Sep 17 00:00:00 2001
From: Matthias Hochsteger <mhochsteger@cerbsim.com>
Date: Tue, 1 Oct 2019 15:15:46 +0200
Subject: [PATCH] Use Mesh::CreatePoint2ElementTable() in SplitImprove()

---
 libsrc/meshing/improve3.cpp | 7 ++-----
 1 file changed, 2 insertions(+), 5 deletions(-)

diff --git a/libsrc/meshing/improve3.cpp b/libsrc/meshing/improve3.cpp
index 2bb098a7..1fd5955a 100644
--- a/libsrc/meshing/improve3.cpp
+++ b/libsrc/meshing/improve3.cpp
@@ -599,7 +599,8 @@ void MeshOptimize3d :: SplitImprove (Mesh & mesh,
   int np = mesh.GetNP();
   int ne = mesh.GetNE();
 
-  TABLE<ElementIndex,PointIndex::BASE> elementsonnode(np); 
+  auto elementsonnode = mesh.CreatePoint2ElementTable();
+
   NgArray<ElementIndex> hasbothpoints;
 
   NgBitArray origpoint(np+1), boundp(np+1);  // big enough for 0 and 1-based
@@ -645,10 +646,6 @@ void MeshOptimize3d :: SplitImprove (Mesh & mesh,
       (*testout) << "Total badness = " << bad1 << endl;
     }
 
-  for (ElementIndex ei : mesh.VolumeElements().Range())
-    for (PointIndex pi : mesh[ei].PNums())
-      elementsonnode.Add (pi, ei);
-
   mesh.MarkIllegalElements();
   if (goal == OPT_QUALITY || goal == OPT_LEGAL)
     {