CompuptNVertices also for 1D

Author: JSchoeberl

libsrc/meshing/meshclass.cpp

@@ -6952,6 +6952,14 @@ namespace netgen
                       },
                       [](auto a, auto b) { return a > b ?  a : b; },
                       numvertices);
+    numvertices =
+      ParallelReduce (LineSegments().Size(),
+                      [&](size_t nr)
+                      {
+                        return int(Max(LineSegments()[nr].Vertices()));
+                      },
+                      [](auto a, auto b) { return a > b ?  a : b; },
+                      numvertices);
     numvertices += 1-PointIndex::BASE;    
   }
 

libsrc/meshing/meshtype.hpp

@@ -1334,6 +1334,7 @@ namespace netgen
     auto PNums() const { return FlatArray<const PointIndex> (GetNP(), &pnums[0]); }
     auto PNums() { return FlatArray<PointIndex> (GetNP(), &pnums[0]); }
 
+    auto Vertices() const { return FlatArray<const PointIndex> (2, &pnums[0]); }
 
     ELEMENT_TYPE GetType() const
     {