Revert "helixish: attempt at the whole thing"

This reverts commit 906edec3ef60cf3e0567c713b2a68688edb634d4.

diff --git a/findcurve.c b/findcurve.c
index c038182e685ca954002d0e0473338077c01f5af4..3c886f6fbab4c14998a9740b64996d883d1e58a4 100644 (file)
--- a/findcurve.c
+++ b/findcurve.c
@@ -122,7 +122,7 @@ int main(void) {
     gsl_siman_solve(rng,
                    startpoint,
                    cb_Efunc, cb_step, cb_metric,
-                   0, //cb_print,
+                   cb_print,
                    0,0,0, sizeof(startpoint), siman_params);
 
     printcore(startpoint);

diff --git a/helixish.py b/helixish.py
index 2c6043eefc306590146f1a05275dcbd4f47f2f19..77c311b34097692d8558208b243baa438ba46e0c 100644 (file)
--- a/helixish.py
+++ b/helixish.py
@@ -163,7 +163,7 @@ class HelixishCurve():
           delta = atan2(*(end - C)[0:2]) - atan2(*(start - C)[0:2])
           dbg('railway_ST C', C, 'start', start, 'end', end, 'delta', delta)
           if delta < 0: delta += tau
-          s = delta * abs(railway_r)
+          s = delta * railway_r
           dbg('railway_ST delta', delta, 'r', railway_r, 's', s)
           return s
 

diff --git a/moebius.py b/moebius.py
index e5a81acad9d0d458bbc823769b5cea55176e0801..f228bce4bdb34171c19b414b5e470e6a32cb52e2 100644 (file)
--- a/moebius.py
+++ b/moebius.py
@@ -81,9 +81,9 @@ class MoebiusHalf:
     m.nu      = nu
     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._beziers = [ m._bezier(theta, HelixishCurve) for theta in m._thetas ]
-    #check = int(nu/3)
-    #m._beziers[check] = m._bezier(m._thetas[check], HelixishCurve)
+    m._beziers = [ m._bezier(theta) for theta in m._thetas ]
+    check = int(nu/3)
+    m._beziers[check] = m._bezier(m._thetas[check], HelixishCurve)
   def _bezier(m, theta, constructor=DoubleCubicBezier):
     cp = [None] * 4
     cp[0] =               m.edge   .point(theta)

diff --git a/visual b/visual
index d623bfb6c14f42f5ea9c2f6097ebd22a9927610e..64d4b6cb41c603ca762eb77d40cf93ccda71bdd7 100755 (executable)
--- a/visual
+++ b/visual
@@ -18,7 +18,7 @@ from moebius import *
 from moedebug import *
 
 nv = 20
-nw = 80
+nw = 40
 
 class Visdebug(BaseVisdebug):
   def curve(vd, fn, color=color.yellow):
@@ -40,7 +40,7 @@ class Visdebug(BaseVisdebug):
     rate(1E-8)
 
 dbg_file(sys.stderr)
-#vdbg_enable(Visdebug())
+vdbg_enable(Visdebug())
 
 m = Moebius(nv, nw)