chiark / gitweb /
helixish: bugfixes
[moebius3.git] / meshscad
1 #!/usr/bin/python3
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 = 1.500;
13
14 # number of wires
15 kv = 24
16 kw = 16
17
18 # resolution
19 nv = 200
20 nw = 200
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 not v % each_v and w < nw:
79           travers[v].rquad([ cnr for cnr in qc(v,w,sigmas,True) ])
80     if 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))