meshscad

   1 #!/usr/bin/python
   2
   3 from __future__ import print_function
   4
   5 import signal
   6 signal.signal(signal.SIGINT, signal.SIG_DFL)
   7
   8 from moebius import *
   9 from scad import *
  10
  11 nomsize = 30;
  12 wire = 2.000;
  13
  14 # number of wires
  15 kv = 7
  16 kw = 4
  17
  18 # resolution
  19 nv = 7
  20 nw = 20
  21 ns = 4 # around tube, should be even
  22
  23 nv += -nv % kv
  24 nw += -nw % kw
  25
  26 each_v = nv / kv
  27 each_w = nw / kw
  28
  29 m = Moebius(nv, nw)
  30
  31 def dpr(v): return '%+3f %+3f %+3f' % tuple(v)
  32
  33 def points_round_wire(p, norm, acrs, sigmas):
  34   for sigma in sigmas:
  35     delta = norm * sin(sigma) + acrs * cos(sigma)
  36     r = p + wire/nomsize * delta
  37     yield r
  38
  39 def calc_sigmas(ss):
  40   return [ (s + 0.5)/ns * tau for s in ss ]
  41
  42 def make_moebius(objname):
  43   print('module %s(){' % objname)
  44   # wires:
  45   extents = [ ScadObject() for w in range(0,nw) ] # along extents
  46   travers = [ ScadObject() for v in range(0,nv) ] # along traverses
  47
  48   def qc(v, w, sigmas, is_trav):
  49     #print(' QCv,w,T',v,w,is_trav, file=sys.stderr)
  50     for ab in 0,1:
  51       p, norm, extt, trav = m.details(v + ab*(not is_trav), w + ab*is_trav)
  52       if is_trav: acrs = extt
  53       else: acrs = trav
  54       #print(' RW,ab,sx',ab,sx,
  55       #      'r=',dpr(r),
  56       #      'p=',dpr(p),
  57       #      'norm=',dpr(norm),
  58       #      'extt=',dpr(extt),
  59       #      'acrs=',dpr(acrs),
  60       #      'delta=',dpr(delta),
  61       #      's=','%4f' % sigma,
  62       #      file=sys.stderr)
  63       for r in points_round_wire(p, norm, acrs, sigmas):
  64         yield r
  65
  66   for v in range(0, nv):
  67     for w in range(0, nw+1):
  68       extw = w
  69       if extw > nw/2: extw = nw-w
  70       # each extent wire has to go round twice to meet itself
  71       # (except the middle one, if there is one, where nw/2 = nw - nw/2
  72       # and it gets to be by itself
  73       for s in range(0, ns):
  74         sigmas = calc_sigmas([s + sx for sx in (0,1)])
  75         #print('VWS',v,w,s, sigmas, file=sys.stderr)
  76         if not w % each_w:
  77           extents[extw].quad([ cnr for cnr in qc(v,w,sigmas,False) ])
  78         if False and not v % each_v and w < nw:
  79           travers[v].rquad([ cnr for cnr in qc(v,w,sigmas,True) ])
  80     if False and not v % each_v:
  81       for w in 0, nw:
  82         p, norm, extt, trav = m.details(v, w)
  83         cnrs = points_round_wire(p, norm, extt, calc_sigmas(range(0,ns)))
  84         cnrs = list(cnrs)
  85         if w: cnrs.reverse()
  86         for s in range(0, ns-1):
  87           travers[v].triangle(cnrs[s],
  88                               cnrs[s+1],
  89                               cnrs[ns-1])
  90
  91   for w in range(0, nw):
  92     print('// extent w=', w)
  93     extents[w].writeout_core(nomsize)
  94   for v in range(0, nv):
  95     print('// travers v=', v)
  96     travers[v].writeout_core(nomsize)
  97   print('}')
  98
  99 make_moebius('MoebiusMesh')
 100 print('moebiuscore_nomsize=%s;' % repr(nomsize))