wire = 2.000;
# number of wires
-kv = 7
-kw = 4
+kv = 24
+kw = 16
# resolution
-nv = 7
-nw = 20
+nv = 80
+nw = 80
ns = 4 # around tube, should be even
nv += -nv % kv
#print('VWS',v,w,s, sigmas, file=sys.stderr)
if not w % each_w:
extents[extw].quad([ cnr for cnr in qc(v,w,sigmas,False) ])
- if False and not v % each_v and w < nw:
+ if not v % each_v and w < nw:
travers[v].rquad([ cnr for cnr in qc(v,w,sigmas,True) ])
- if False and not v % each_v:
+ 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)))