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 = 50;
  12 wire = 3.0;
  13
  14 # number of wires
  15 kv = 24
  16 kw = 16
  17
  18 # resolution
  19 nv = 80
  20 nw = 80
  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       for s in range(0, ns):
  69         sigmas = calc_sigmas([s + sx for sx in (0,1)])
  70         #print('VWS',v,w,s, sigmas, file=sys.stderr)
  71         if not w % each_w:
  72           extents[w % nw].quad([ cnr for cnr in qc(v,w,sigmas,False) ])
  73         if not v % each_v and w < nw:
  74           travers[v].rquad([ cnr for cnr in qc(v,w,sigmas,True) ])
  75     if not v % each_v:
  76       for w in 0, nw:
  77         p, norm, extt, trav = m.details(v, w)
  78         cnrs = points_round_wire(p, norm, extt, calc_sigmas(range(0,ns)))
  79         cnrs = list(cnrs)
  80         if w: cnrs.reverse()
  81         for s in range(0, ns-1):
  82           travers[v].triangle(cnrs[s],
  83                               cnrs[s+1],
  84                               cnrs[ns-1])
  85
  86   for w in range(0, nw):
  87     print('// extent w=', w)
  88     extents[w].writeout_core(nomsize)
  89   for v in range(0, nv):
  90     print('// travers v=', v)
  91     travers[v].writeout_core(nomsize)
  92   print('}')
  93
  94 make_moebius('MoebiusMesh')
  95 print('moebiuscore_nomsize=%s;' % repr(nomsize))