From 7862ac9a12a97686ad164b5e266bbd113c38538a Mon Sep 17 00:00:00 2001
From: Joachim Schoeberl <joachim.schoeberl@tuwien.ac.at>
Date: Mon, 22 Mar 2010 21:41:07 +0000
Subject: [PATCH] points on 2D boundary

---
 libsrc/geom2d/genmesh2d.cpp | 1 +
 libsrc/geom2d/spline.hpp    | 6 +++++-
 2 files changed, 6 insertions(+), 1 deletion(-)

diff --git a/libsrc/geom2d/genmesh2d.cpp b/libsrc/geom2d/genmesh2d.cpp
index bec84518..b792ce39 100644
--- a/libsrc/geom2d/genmesh2d.cpp
+++ b/libsrc/geom2d/genmesh2d.cpp
@@ -39,6 +39,7 @@ namespace netgen
 
     geometry.PartitionBoundary (h, *mesh);
 
+
     // marks mesh points for hp-refinement
     for (int i = 0; i < geometry.GetNP(); i++)
       if (geometry.GetPoint(i).hpref)
diff --git a/libsrc/geom2d/spline.hpp b/libsrc/geom2d/spline.hpp
index 1487d623..0b0095e2 100644
--- a/libsrc/geom2d/spline.hpp
+++ b/libsrc/geom2d/spline.hpp
@@ -345,8 +345,12 @@ void SplineSeg<D> :: Partition (double h, double elto0,
       while (j < curvepoints.Size() && (l >= curvepoints[j] || i == n))
 	{
 	  frac = (curvepoints[j]-lold) / (l-lold);
-	  mark = pold + frac * (p-pold);
 	  edgelength = i*dt + (frac-1)*dt;
+	  // mark = pold + frac * (p-pold);
+	  mark = GetPoint (edgelength);
+	  
+	  // cout << "mark = " << mark << " =?= " << GetPoint (edgelength) << endl;
+
 	  {
 	    PointIndex pi1 = -1, pi2 = -1;