curveopt: test a bit

Signed-off-by: Ian Jackson <ijackson@chiark.greenend.org.uk>

diff --git a/moebius.py b/moebius.py
index 9235a7755ca7764ba7192a3d96b22f6832ef533a..b3e32c5565a8eb5d9b84dfabb86e53ae63389832 100644 (file)
--- a/moebius.py
+++ b/moebius.py
@@ -7,6 +7,7 @@ from numpy import cos, sin
 from bezier import BezierSegment
 from moenp import *
 from moebez import *
+from curveopt import *
 
 import sys
 
@@ -60,6 +61,10 @@ class MoebiusHalf:
     m._thetas = [ u * tau for u in np.linspace(0, 1, nu+1) ]
     m._cp2b = BezierSegment([ (c,) for c in [0.33,0.33, 1.50]])
     m._dbeziers = [ m._dbezier(theta) for theta in m._thetas ]
+    check = int(nu/3)-1
+    m._dbezier(m._thetas[check], OptimisedCurve)
+    #m._dbezier(m._thetas[check+1], HelixishCurve)
+    #m._dbezier(m._thetas[check+2], HelixishCurve)
   def _dbezier(m, theta, dconstructor=DiscreteBezier):
     cp = [None] * 4
     cp[0] =               m.edge   .point(theta)