- extents[w].quad([ cnr for cnr in qc(v,w,sigmas,False) ])
- if False and not v % each_v:
- travers[v].quad([ round_wire(p, norm, extt, sigmas[sx])
- for b in 0,1
- for p, norm, extt, trav in (m.details(v, w+b),)
- for sx in 0,1 ])
- for v in range(0, nv):
- print('// extent v=', v)
- extents[v].writeout_core(nomsize)
+ extents[w % nw].quad([ cnr for cnr in qc(v,w,sigmas,False) ])
+ if not v % each_v and w < nw:
+ travers[v].rquad([ cnr for cnr in qc(v,w,sigmas,True) ])
+ if not v % each_v:
+ for w in 0, nw:
+ p, norm, extt, trav = m.details(v, w)
+ cnrs = points_round_wire(p, norm, extt, calc_sigmas(range(0,ns)))
+ cnrs = list(cnrs)
+ if w: cnrs.reverse()
+ for s in range(0, ns-1):
+ travers[v].triangle(cnrs[s],
+ cnrs[s+1],
+ cnrs[ns-1])
+