from moebius import *
from moedebug import *
-nv = 40
-nw = 40
+nv = 20
+nw = 80
class Visdebug(BaseVisdebug):
def curve(vd, fn, color=color.yellow):
rate(1E-8)
dbg_file(sys.stderr)
-vdbg_enable(Visdebug())
+#vdbg_enable(Visdebug())
m = Moebius(nv, nw)
#c0 = curve( color=color.red, pos = [ m.edge.point(th) for th in thetas ] )
#c1 = curve( color=color.blue, pos = [ m.midline.point(th) for th in thetas ] )
-for v in vs:
- c2 = curve( color=color.blue, pos = [ m.point(v,w) for w in ws ] )
+for vi in range(0, len(vs)):
+ v= vs[vi]
+ c2 = curve( color= (color.blue if vi == int(2*nv/3) else (0,0,0.5)),
+ pos = [ m.point(v,w) for w in ws ] )
#for w in ws:
# c3 = curve( color=color.blue, pos = [ m.point(v,w) for v in vs ] )