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])
+
for w in range(0, nw):
print('// extent w=', w)
extents[w].writeout_core(nomsize)